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20.  Abstract  (continued) 

The  model  is  applied  to  simulate  the  freeze -dehydration  of  beef  meat  by 
microvmve  dielectric  heating  at  24^0  MHz.  The  results  of  the  calculations 
show  that  the  tricrowave  freeze-drying  process  is  mainly  controlled  by 
heat  transfer.  The  study  of  the  effects  of  the  process  constraints  and 
variables  upon  the  drying  time  and/or  the  freeze-dryer  output  indicates 
that  an  optimal  operation  (highest  product  yield)  would  be  for  an  opera- 
tion near  corona  and  overheating  (or  melting)  conditions.  Drying  times 
as  short  as  1^  hours  should  be  possible,  if  corona  is  to  occur  beyond 
a peak  strength  of  255  V/cm. 

E3q)erlmental  drying  curves  have  been  obtained  for  the  freeze -drying  of 
beef  meat  with  microwave  heating  at  2450  MHz.  The  experimental  drying  curves 
are  in  good  agreejnent  with  those  simulated  by  the  model.  The  small  dif- 
ferences may  be  attributed  to  an  uncertainty  on  the  numerical  values  used 
for  the  physical  properties  in  the  model  or  non-ideal  eacperimental  condi- 
tions. An  experimental  value  of  the  heat  transfer  coefficient  of 
1.9  X 10“4  cal/sec/cm^/°C  (at  the  surface  of  the  sanple  exposed  to  the 
vacuum)  is  also  obtained. 
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Dehydration  has  been  used  for  millenia  to  extend  storage  life  of  food 
products.  If  properly  done  and  followed  by  appropriate  packaging,  a wide 
variety  of  food  products  can  be  maintained  without  significant  change  at 
ambient  temperatures  for  very  long  periods.  Wei^t  savings  may  be  as  hi^ 
as  90^.  Accordingly,  dehydration  has  always  played  an  important  role  in  the 
preparation  of  military  rations. 

Efforts  to  improve  the  process  have  culminated  in  freeze-diying  or 
dehydration  by  sublimation.  The  results  have  been  dramatic  and  with  many 
products  consumer  acceptance  approaches  that  for  fresh  or  frozen  counterparts, 
'flie  cost  of  freeze-dr„  'ng  is  hi^  due  to  the  equipment  required  and  low 
through-put  rate.  Present  procedures  are  alow  because  as  the  product  dries 
from  the  surface  inward  the  energy  for  sublimation  must  pass  by  conduction 
through  dry  product  which  is  an  effective  thermal  insulator.  Transfer  of 
microwave  energy  is  nearly  unaffected  by  the  dry  layer.  Thus,  theoretically, 
microwave  freeze-drying  should  permit  drying  much  more  rapidly. 

Investigations  of  microwaVe  freeze-dr.ving  have  indicated  many  problems. 

This  study  covers  theoretical  and  ejq)erimental  approaches  toward  their  solution. 
Included  is  a mathematical  simulation  of  the  drying  process  which  it  is 
hoped  will  constitute  a valuable  tool  for  future  process  design,  optimization 
and  con'’  'ol. 

This  study  was  undertaken  under  a grant  from  US  Army  Natick  Laboratories 
No.  DAAG37-71-G-0001,  Project  Number  1T161101A91A07,  Ix>*House  Laboratory 
Independent  Research.  This  manuscript  was  previously  submitted  to 

Worcester  Pol’^iechnic  Institute  in  partial  fulfillment  oi  requirements  for 
the  PhD  degree  in  Chemical  Engineering.  The  project  Officer  for  this 
grant  was  Dr.  Clarence  K.  Wadsworth. 
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- Chapter  I 


INTRODUCTION 


For  many  years,  dehydrated  food  products  have  received  a lot  of 
attention  as  they  greatly  simplify  storage  and  conservation  problems. 

In  effect,  dehydrated  food  products  can  be  kept  at  room  temperature 
over  long  periods  of  time  if  properly  packaged.  On  the  other  hand, 
the  loss  of  the  moisture  content  represents  a saving  in  weight  of 
approximately  75%. 

Various  means  of  dehydration  have  been  developed  to  suit  the 
needs  of  the  industry.  Among  these,  the  freeze -dr ying  process  is  of 
special  interest  as  it  is  most  suited  for  heat  sensitive  materials 
such  as  biologicals  and  food  products. 

In  this  technique,  which  was  developed  in  the  early  1950 's  the 
material  to  be  dried  is  first  deep  frozen  and  then  dehydrated  by 
sublimation  of  its  frozen  moisture  content.  As  the  dehydration  pro- 
ceeds an  outer  dried  layer  forms  while  the  frozen  core  retreats.  It 
has  been  observed  experimentally  that  a sharp  sublimation  zone  separates 
distinctly  the  nearly  completely  dried  layer,  and  the  fully  hydrated 
and  frozen  core.  The  sublimation  is  maintained  by  keeping  the  sub- 
limation zone  temperature  below  the  melting  point  of  the  frozen  material. 
A low  partial  pressure  of  water  vapor  in  the  surrounding  atmosphere 
provides  the  driving  force  necessary  for  the  dehydration  to  take 
place.  Water  vapor  partial  pressures  of  the  order  of  .1  mmHg  are 
usually  maintained  by  means  of  a low  temperature  condenser  (-40  to  -50°C) . 
The  dehydration  is  generally  conducted  under  vacuum  (.1  to  1 mmHg) 
in  order  to  decrease  the  resistance  to  the  water  vapor  flow  through 
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tht  porous  dried  layer,  although,  a sublination  under  atmospheric 
conditions  is  physically  possible  (Hohner  (1970)). 

Although  it  has  proven  very  successful  and  is  widely  used, 
conventional  freese-drying  processes  use  an  external,  surface 
heating  (i.e.,  hot  plates,  infra  reo  lamps)  to  supply  the  necessary 
enthalpy  of  sublimation  to  the  product.  The  poor  thermal  conductivity 
of  food  products  results  in  long  drying  times  due  to  upper -limits  on 
the  temperature  gradients  allowable  in  the  product  in  order  to  avoid 
melting  and/or  thermal  degradation  of  the  dried  material.  Overheating 
of  the  dried  product  is  particularly  undesirable  with  respect  to  the 
reversibility  of  the  f reese-arying  process  as  poor  rehydration 
characteristics  may  develop  (Kuprianoff  (1961)). 

Microwave  dielectric  heating  has  shown  a great  potential  in 
accelerating  f reeee-drying  (Hoover  et  al.  (1966)).  Its  feasibility 
was  demonstrated  on  an  experimental  scale  in  the  late  1950* s by 
Jackson  et  al.  (1957)  and  Copson  (1958).  Hoover  et  al.  (1966) 
have  shown  that  use  of  microwave  energy  to  replace  conventional 
(external)  heating  may  reduce  drying  times  by  factors  of  3 - 13 
which  would  allow  the  use  of  a continuous  or  at  least  semi -continuous 
process.  The  shorter  drying  times  are  due  to  the  fact  that  microwave 
energy  is  dissipated  throughout  the  bulk  of  the  material  being  dried, 
particularly  in  the  frozen  core.  Thus,  the  problem  of  supplying 
the  enthalpy  of  sublimation  by  conduction  from  the  surface  can  be 
overcome. 

However,  the  microwave  freeze-drying  process  has  remained  in  the 
experimental  stage.  The  lack  of  a systematic  theoretical  study  of 
the  process  as  well  as  problems  with  the  relatively  new  microwave 
heating  technology  seem  to  have  hampered  the  development  of  this 
new  technique.  Among  the  problems  involved  in  the  design  of  a 
microwave  freeze-dryer,  the  most  critical  ones  are  the  corona  dis- 
charge occuring  at  low  pressures,  the  difficulty  of  obtaining  a 
uniform  heating  of  the  frozen  product,  and  the  difficulty  of  coupling 
efficiently  microwave  energy  with  the  material  being  dried. 

The  question  of  corona  discharge  in  the  microwave  applicator  is 
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of  prim*  importanc*  at  ita  formation  drastically  d*cr*at*t  th*  power 
•fficiancy  and  may  reault  in  a«v*r*  burnt  of  th*  product.  Thus,  th* 
electric  field  in  the  microwave  applicator  mutt  be  kept  below  its 
break-down  value.  This  conatitutea  an  additional  constraint  imposed 
on  th*  proceta.  Previous  studies  have  shown  that  th*  electric  field 
strength  at  which  breakdown  (corona  discharge)  occurs  is  a function 
of  the  pressure  In  the  vacuum,  the  design  of  the  microwave  generator 
and  the  geometry  of  the  microwave  applicator  (MacDonald  (1966), 

Gould  & Kenyon  (1971)).  Experimental  results  Ly  Gould  S.  Kenyon  on 
the  phenomenon  of  gaseous  breakdown  with  reference  to  f reeze-drying 
indicate  breakdown  field  strengths  of  200  to  300  V.^cm  RMS  in  the 
pressure  range  2 to  .1  mmHg.  In  this  pressure  range,  the  highest 
strengths  are  reached  for  the  lowest  pressure  valies. 

Research  has  been  undertaken  at  the  U.S.  Army,  Natick  Laboratories, 
in  an  effort  to  develop  the  microwave  freeze-drying  technique. 

As  part  of  this  work,  the  present  investigation  is  intended  to 
provide  a mathematical  simulation  of  the  freeze-drying  process  using 
microwave  dielectric  heating.  It  will  constitute  a useful  tool 
in  the  future  for  process  design,  optimization  and  control. 

The  objectives  of  this  work  can  be  stated  as  follows: 

1.  Develop  a mathematical  mbdel  to  describe  the  process. 

2.  Study  the  response  of  the  mathematical  model  to  various 
freeze-drying  conditions  in  order  to  test  its  validity 
and  provide  design  guidelines.  Of  prime  concern  will  be 
the  effects  of  the  process  constraints  and  variables  on 
the  drying  time. 

3.  Obtain  experimental  drying  curves  of  beef  samples  and 
compare  them  with  predicted  drying  curves  as  a verification 
of  the  simulation. 


- Chapter  II 


LITERATURE  SURVEY 


During  the  last  twenty-five  years,  there  have  been  numerous 
publications  on  the  various  asjjetts  of  conventional  freeze  drying. 
However,  very  little  information  could  be  found  on  microwave  freeze- 
drying. It  should  be  pointed  out  that  there  still  is  a need  for 
a systematic  theoretical  study  of  the  conventional  f reeze-drying 
process . 

1 . The  microwave  freeze-dr  vino  process 

Most  of  the  work  done  on  this  subject  falls  in  a few  isolated 
experimental  studies  which  were  directed  toward  the  development  of 
this  new  technique.  Jackson  et  al.  (1957),  Copson  (1962)  successfully 
accelerated  freeze-drying  with  microwave  energy  on  an  experimental 
scale.  They  found  that  microwave  heating  should  allow  drastic 
reductions  of  the  drying  time,  although  engineering  prob’’  ,ms,  particu- 
larly glow  discharge  may  complicate  its  application.  Harper  (1961), 
Kan  & Yeaton  (1961)  obtained  experimental  data  on  the  dielectric 
properties  of  various  food  products  between  microwave  frequencies 
of  500  and  3000  MHz.  Hoover  et  al.  (1966)  obtained  comparative 
drying  curves  for  radiant  and  microwave  (at  632  and  915  MHz)  freeze- 
drying of  various  food  products.  They  demonstrated  that  use  of 
microwave  freeze-drying  may  reduce  drying  times  by  factors  of  3 to  13 
as  compared  to  radiant  (surface)  heating.  They  found  that  with 
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mlcrowavt  hMtlng  th«  drying  tln«i  wert  tiientially  Indeptndtnt  of 
tht  thickneii  of  th«  lampl*.  It  should  bt  notsd  that  th«  mlcrowava 
power  Input  level  was  Increased  in  their  experiments,  In  the  same 
ratio  as  that  of  the  thickness.  They  reported  variations  of  the 
glow  power  vs  pressure  for  the  microwave  applicator  which  they  used. 

In  a succeeding  paper,  they  used  engineering  data  derived  from  their 
experimental  results  to  discuss  the  feasibility  of  the  microwave 
freeze -drying  process. 

General  analyses  about  the  microwave  freeze-drying  process 
have  been  reported  in  the  recent  years.  Copson  (1962)  published 
a book  in  which  he  described  in  detail  the  theoretical  and  engineering 
prob]  'ms  arising  in  the  application  of  microwave  heating  to  freeze- 
drying. Parker  (1968)  made  an  extensive  use  of  the  data  of  Hoover 
et  al . , to  illustrate  the  possibilities  of  microwave  freeze-drying, 
as  one  of  the  multiple  applications  of  microwaves  in  Power  Engineering 
Other  reviews  have  been  done  by  Decareau  (1962),  Meryman  (1964), 
and  Harper  et  al . (1962). 

Based  on  Hoover’s  experimental  data  for  microwave  freeze-drying 
at  915  MHz,  Hammond  (1967)  made  an  economic  evaluation  of  microwave 
vs  radiant  heating  and  concluded  that  both  methods  were  economically 
competitive.  However,  he  also  pointed  out  that  his  analysis  was 
very  sensitive  to  the  assumptions  which  he  made. 

More  recently,  Gould  and  Kenyon  (1971)  studied  experimentally 
the  problem  of  corona  breakdown  in  microwave  freeze-drying  at  2450  MHz 
They  also  derived  an  equation  which  shows  the  effect  of  a dielectric 
load  on  the  gas  discharge  in  a microwave  cavity. 

2.  Modeling  the  f reeze-drvinq  process 

There  has  not  been,  yet,  any  comprehensive  analysis  of  the 
freeze-drying  process  using  microwave  heating.  Previous  analyses 
have  been  limited  to  simplified  approaches  using  a quasi -state 
assumption  and  were  aimed  at  preliminary  design  applications  (Copson 
(1962)). 

On  the  other  hand,  there  have  been  many  publications  concerned 
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with  a theoretical  description  of  the  conventional  freeze-drying 
process  using  external  heatir*;.  The  reader  is  referred  for  more 
information  to  the  excellent  review  of  King  (1971),  and  an  extensive 
collection  of  articles  on  Lyophilisat ion  edited  by  Rey  (1964). 

Although  theoretical  studies  of  conventional  freeze  dr ving  are 
not  applicable  to  microwave  freeze-drying  (this  point  will  be 
discussed  below),  it  is  of  interest  to  review  some  of  the  most 
recent  and  complete  studies  as  they  emphasize  the  need  for  a general, 
transient  analysis  of  the  freeze-drying  process. 

These  studies  have  ge)erally  employed  the  simplifying  assumption 
of  a pseudo-steady  state.  This  assumption  is  not  applicable  to 
microwave  freeze-drying  because  of  the  fast  temperature  transients 
occu.’'ing  in  the  sample.  However,  this  assumption  is  generally  quite 
good  in  the  radiant  heating  case  as  shown  by  Sandall  et  al.  (1967). 
Hill  and  Sunderland  '1971)  have  presented  a meoretical  analysis 
considering  the  simultaneous  energy  anu  ma.s  transfer  during  dehy- 
dration. Although  they  considered  the  time  variation  of  the  interface 
temperature,  they  have  neglected  the  trar  -,ient  behavior  of  temperature 
and  concentration  both  in  the  ice-core  aid  in  the  dried  layer. 

Fox  and  Thompson  (19^2)  considered  the  transient  transport  equations 

in  the  dried  layer  for  two  limiting  cases:  quasi-steady  state  with 

0 

heat  leakage  in  the  ice-core  and  transient  solutions  with  a uniform 
ice-core  temperature.  Hohner  (1970)  considered  the  transient 
behavior  of  both  mass  and  energy  transport  equation  with  a numerical 
technique. 

Thus,  there  is  seemingly  a need  for  a general  analysis  of  the 
freeze-drying  process  which  would  account  for  all  transient  effects 
in  both  the  frozen  and  the  dried  regions.  Such  a study  would  permit 
the  simulation  of  the  radiant  freeze-drying  process  and  thus,  allow 
a systematic  theoretical  study. 

The  objective  of  this  work  is  to  develop  a general  study  using 
microwave  dielectric  heating.  Microwave  freeze -drying  is  indeed  a 
more  general  case.  In  effect,  it  includes  a heat  generation  term 
In  the  energy  equation,  whereas,  conventional  freeze-drying,  using 
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•urfac*  hMtlng  do«t  not.  A nodal  of  tho  mlcxowava  f rtesa-drylng 
pxocoiis  would  b«  xaadlly  applieabla  to  tha  eonvantlonal  ona  by 
aattlng  tha  volunatxlc  hMt  ganaxatlon  taxm  to  aaxo  and  adjusting 
the  boundaxy  conditions  of  the  equations  governing  the  process  to 
account  for  the  surface  heating. 

3 . Physical  properties  of  beef  meat  vitn  reference  to  microwave 
f reeze-drvino 

3 . 1 Thermal  conductivities 

The  thermal  conductivity  of  dried  beef  .meal  nas  been  measured 
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by  Harper  (1961,  1962)  at  various  pressures  in  the  range  10  - lu  rotnHg. 

The  measurements  were  made  using  heat  flow  transducers,  located  on 
both  sides  of  the  sample.  He  also  derived  an  equation  to  predict 
the  thermal  conductivity  of  the  porous  material  which  showed  very 
good  agreement  with  his  experimental  data.  Other  thermal  conductivity 
values  have  been  reported  for  dried  beef.  However,  they  are  generally 
the  results  of  calculations  from  the  drying  rate  data  and  may  lack 
accuracy  as  they  depend  on  the  assumptions  employed  in  the  calculations 
(Massey  and  Sunderland  (1967)). 

Quite  a few  data  (see  Chapter  6,  Sec.  2)  have  been  reported 
on  the  thermal  conductivity  of  frozen  beef  meat  as  it  is  of  prime 
importance  in  refrigeration  applications.  The  various  literature 
sources  have  been  listed  in  Table  2.1  with  the  composition  of  the 
beef  samples  used  in  the  measurements. 

Table  2.1  - Sources  of  thermal  conductivity 
data  for  frozen  beef 


Sour  re 

% moisture 

% fat 

Lentz  (1961) 

75 

.9 

Miller  & Sunderland 

69.5 

(1963) 

Hill  et  al.  (1967) 

76.5 

2.35 

Hill  et  al.  (1967) 

78.7 

1.4 
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These  data  have  generAlly  been  obtained  by  measuring  the  boat  flux 
through  a slab  of  the  samole  placed  between  a cold  and  a hot  plate. 
Accuracies  better  than  3%  are  quoted  by  the  authors. 

3 . 2 Dif f usivities 

The  experimental  data  on  the  effectiv'e  diffusivity  of  beef 
meat  are  very  sparse.  They  are  practically  limited  to  measurements 
by  Harper  (1962)  and  Dyer  & Sunderland  (1966).  Harper  (1962)  has 
measured  the  diffusivity  of  water  vap>or  in  dried  beef  at  38 °C 
under  1 atm.  He  mentirned  that  the  measurements  were  not  easy 
to  reproduce  and  that  there  was  an  uncertainty  of  10  to  1536  on  the 
results.  The  measurements  consisted  of  recording  the  gain  in  weight 
of  sodium  perchlorate  placed  above  the  sample  through  which  water 
vapor  was  allowed  to  diffuse  under  atmospheric  pressure. 

Dyer  & Sunderland  (1966)  have  obtained  diffusivity  data  in  the 
range  1.44  to  2.29  mmHg , at  an  unspecified  temperature.  These 
data  were  compared  to  the  results  of  a theory  which  they  derived, 
to  predict  effective  diff usivities  in  the  transition  flow  regime. 

The  agreement  is  of  the  order  of  15  - 20%.  They  measured  the  v;ater 
vapor  flow  by  means  of  humidity  transducers  located  on  both  sides 
of  the  sample  v'thin  an  evacuated  chamber. 

3.3  Dielectric  properties 

There  has  been  data  reported  in  the  literature  on  the  dielectric 
properties  of  beef  meat,  and  other  food  products.  A review  of  these 
data  has  been  published  by  Nelson  (1971).  However,  most  of  these 
data  have  been  obtained  for  microwave  heating  applications  and  are 
usually  limited  to  the  fresh  product  at  temperatures  above  freezing. 
Indeed,  no  data  could  be  found  for  frozen  and  dried  beef  meat  at 
2450  MHz.  The  most  readily  applicable  data  were  those  reported 
by  Kan  & Yeaton  (1961)  at  3000  MHz.  These  data  cover  the  range 
•59.4*C  to  4.4*C  for  frozen  beef  (in  5.6'C  increments)  and  -51,1*C 
to  82.2*C  (in  11.2'’C  increments).  The  measurements  were  made  using 
the  shorted  line  method  with  coaxial  standing  wave  detect  'rs. 
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Results  of  measurements  at  1000  MHz  and  150  MHz  are  also  reported. 

Additional  data  points  have  been  obtained  by  Harper  et  al. 
(1962)  for  frozen  beef  at  -4*  and  -10*C  for  a frequency  of  2000  MHz. 
Bengtsson  I*  Risman  (1971)  also  reported  measurements  for  frozen 
raw  beef  at  2.8  MHz,  by  the  cavity  perturbation  technique,  at  two 
temperatures:  -20‘’C  and  -10“C. 

3.4  Specific  heat  of  frozen  and  dried  beef 

Experimental  data  for  frozen  and  dried  beef  have  been  reported 
by  Awbery  & Griffith  (1933)  and  Rey  (1964). 


I 
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PART  I 


miHEMATICAL  SIMULATION  OF  A FREEZE -DRYING  PROCESS  USING 

MICROWAVE  ENERGY 


Chapter  III:  The  Freeze-Drying  Process 
Chapter  IV:  The  Mathematical  Model 
Chapter  V:  The  Method  of  Solution 


Chapter  III 


THE  FREEZE -DRYING  PROCESS 


Although  this  investigation  has  been  conducted  with  reference 
to  the  f reere-dehydration  of  beef  meat,  it  can  easily  be  extended 
to  the  f reere-drying  of  other  products,  The  simulation  can  also 
be  applied  to  separation  processes  involving  sublimation.  Dielectric 
heating  is  considered  in  this  investigation  but,  it  should  be 
noted  that  it  does  not  limit  the  applicability  of  the  analysis  to 
the  specific  area  of  dielectric  heating.  On  the  contrary,  it 
deals  with  a more  general  case  than  the  conventional  freeze-dehydration 
process  using  external  heating.  In  effect,  the  present  analysis 
includes  a volumetric  heat  generation  term  in  the  governing  energy 
balance  equations.  The  present  analysis  can  be  used  to  simulate 
the  conventional  (surface  heating)  f reezo-drying  process  by  choosing 
the  appropriate  boundary  conditions. 

However,  in  extreme  cases  where  the  process  parameters  are 
drastically  different  from  those  assumed  here,  a revision  of  some 
of  the  assumptions  may  be  necessary.  For  example,  for  much  greater 
pore  sizes  and/or  higher  operating  pressures,  the  description  of 
the  water  vapor  flow  through  the  dried  layer  by  Fick*s  equation 
may  not  be  adequate  as  hydrodynamic  flow  predominates. 

1 . General  description 

In  the  freeze-drying  technique,  the  material  to  be  dried  is 
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first  dstp  frossn  (around  -35*C)  ivnd  than  dahydrattd  by  dirsct 
sublimation  of  th«  frosan  moistura  absorbad  in  its  pora  structure. 

As  tha  dahydration  proceeds,  the  frozen  core  retreats  inward  while  . 
a porous,  dried  outer  layer  forms.  A sublimation  zone  separates 
distinctly  the  frozen  core  and  the  dried  layer,  as  observed  ex- 
perimentally by  this  investigator  and  other  previous  studies 
(Clark  (19b8),  Hatcher  (1964),  Hardin  (1965)). 

For  the  sublimation  to  take  place,  the  temperatures  in  the 
sublimation  zone  must  be  kept  below  the  triple  iint  of  the  absorbed 
moisture  (around  -l.S*C  for  beef  meat  ,1 . The  water  vapor  from  the 
phase  change  at  the  frozen  front  flows  through  the  porous  dried 
layer  to  the  surrounding  atmosphere  where  it  is  condensed  by  a 
vapor  trap  (moisture  sink).  Low  operating  pressures  are  usually 
used  (.05  to  ImmHg)  in  order  to  insure  a small  resistance  to  mass 
transfer.  In  tha  case  under  investigation,  the  heat  necessary 
for  sublimation  is  supplied  to  the  material  by  means  of  dielectric 
heating  at  a microwave  frequency  of  2450  MHz. 

2.  Sharpness  of  the  sublimation  zone 

The  assumption  of  a uniformly  retreating  ice  front  with  a 
sharp  boundary  (URIF  Model),  has  been  commonly  used  by  previous 
investigators,  and  is  generally  in  good  agreement  with  experiments. 
However,  some  questions  have  been  raised  regarding  the  degree  of 
sharpness  of  the  front.  Bralsford  (1967),  while  using  the  URIF 
Model,  suggests  that  liquid  diffusion  is  likely  to  broaden  the 
sublimation  zone,  particularly  as  the  temperature  approaches  the 
melting  point  of  the  frozen  material.  Moisture  measurements  by 
Brajnikov  et  al.  (1969),  and  Hatcher  (1964),  indicate  the  pos- 
sibility for  a diffuse  sublimation  zone  as  thick  as  1 to  2 mm. 

But,  as  suggested  by  King  (1971),  it  could  reflect  the  limit  of 
resolution  of  the  moisture  measurement  method  which  they  used. 

Othex  evidences  for  a diffuse  front  have  been  reported  and  are  dis- 
cussed by  King  (1971). 

On  the  other  hand,  visual  examination  of  partially  dried  meat 
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samples  (Hardin  (19f5))  and  measurement  of  temperatui  ! histories 
during  the  drying  process  by  Hatcher  (1964)  have  shown  a sharp 
transition  zone.  It  should  be  noted  that  Hatcher'';  results  seem 
to  show  some  broadening  of  the  transition  zone  i.i  the  latter  stage 
of  drying,  although  they  could  be  interpreted  by  a non -planar  but 
still -sharp  front.  An  additional  fact  in  support  of  the  snarp 
frozen  front  is  that,  as  already  mentioned,  models  of  conventional 
freeze-drying  which  have  used  the  URIF  Model  of  an  infinitely 
sharp  front  have  been  in  good  agreement  with  experimental  results 
(Bralsford  (1967),  Sandall  et  al.  (1967)). 

Thus  the  hypothesis  of  a sharp  sublimation  front  seems  to  be 
a reasonable  assumption  for  the  description  of  the  freeze -drying 
mechanism,  despite  the  fact  that  some  broadening  of  the  front  may 
occur  as  the  temperature  approaches  the  melting  point. 

3.  Heat  transfer 

An  important  characterization  of  the  freeze-drying  process  is 
its  high  energy  needs.  As  compared  to  other  processes,  such  as 
drying  by  vaporization,  f reeze-drying  requires  a greater  amount  cf 
energy.  In  effect,  the  heat  of  sublimation  corresponds  to  the 
sum  of  the  heat  of  vaporization  (595cal/g  at  O'C)  and  the  heat  of 
fusion  (80cal/g)  and  is  675  cal/g  for  ice. 

Thus,  one  may  expect  that  the  rate  of  heat  transfer  to  the 
sublimation  front  will  be  a determining  factor  for  the  total  drying 
time.  Indeed,  vacuum  f reeze-drying  processes  are  heat  transfer 
limited  as  shown  by  studies  with  conventional,  surface  heating 
and,  microwave  dielectric  heating  (Ma  & Peltre  (1973)).  This 
is  due  to  tha  fact  that  the  maximum  temperature  is  to  be  kept  below 
60*C  to  avoid  thermal  degradation  in  the  dried  layer,  and  below 
-1.5*C  in  the  frozen  region  to  eliminate  melting  or  vaporization. 

This  results  in  an  upper  limit  for  allowable  temperature  gradients. 

The  sublimation  front  acts  as  a heat  sink  toward  which  the 
heat  flows  from  both  the  frozen  and  dried  regions.  The  heat  transfer 
in  both  regions  is  essentially  due  to  conduction.  The  use  of 


effective  thermal  conductivities  and  effective  specific  heats  enables 
one  to  treat  the  porous  medium  as  a continuum.  A small  but  appre- 
ciable amount  of  heat  is  lost  to  the  sublimation  in  the  dried  layer 
due  to  convection  as  shown  by  Dyer  & Sunderland  (1968). 

Regarding  the  very  low  thermal  conductivity  of  the  solid 
matrix,  the  gas  filling  the  pores,  in  the  dried  layer,  is  expected 
to  have  a significant  contribution  to  the  overall  thermal  conductivity 
of  the  dried  layer.  Indeed,  experimental  results  (Harper  (1962), 
Sandall  et  al.  (1967))  show  the  effective  thermal  conductivity  in 
the  dried  layer  to  be  a function  of  pressure.  On  the  other  hand, 
the  effective  thermal  conductivity  in  the  frozen  region  varies 
significantly  with  the  temperature  (Lentz  (1961)).  These  effects 
have  been  taken  into  account  in  the  derivation  of  the  model. 


4.  The  microwave  power  conversion 

With  the  microwave  dielectric  heating,  energy  is  generated 
throughout  the  volume  of  the  material  due  to  the  interaction  of 
tho  electric  field  with  the  free  and  bound  charges  of  the  dielectric. 
The  magnetic  component  of  the  electromagnetic  field  does  not 
generate  any  energy  because  its  direction  is  normal  to  the  velocity 
of  the  charges.  From  these  considerations,  and  assuming  a linear 
and  conductive  dielectric,  a formula  can  be  derived  to  relate 
the  time-averaged  power  generated  in  a unit  volume  of  dielectric 
to  the  electric  field  peak  strength  in  the  dielec  ric,  E^, 


u = 


*^e  p2 

T^d 


(3.1) 


It  can  be  further  transformed  by  expressing  the  electric  conductiv- 
ity in  terms  of  the  loss  factor  of  the  dielectric  €"  (also  known 
as  the  imaginary  part  of  the  complex  dielectric  constant). 

u = TT  f € €"  (3.2) 
o d ' ' 


where  f is  the  frequency  of  the  wave. 

If  a polarization  of  the  electric  field  normal  to  the  surface 
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oi  the  dielectric  material  exists,  then  the  peak  strength  in  the 
dielectric  near  the  outer  surface  can  be  related  to  that  in  the 
surrounding  vacuum  by, 

•2  '*2  U 

E = E/  ^ + e (3.3) 

d 

On  the  other  hand,  if  the  field  is  tangent  to  the  surface, 
it  is  continuous.  Then,  one  has, 

= E (3.4) 

Equation  (3.3)  is  used  as  a dominant  polarization  of  the 
field  normal  to  the  sample  surface  is  expected  in  the  experiments 
(see  Chapter  VIII,  Sec.  2.3).  However,  it  should  be  noted  that 
the  true  situation  must  lie  between  these  two  extreme  cases, 
depending  on  the  amount  and  the  direction  of  polarization  of  the 
electric  field. 

If  one  neglects  the  attenuation  effect,  the  strength  of  the 
electric  field  is  constant  throughout  the  sample  and  equal  to 
the  value  at  the  surface  (for  a uniform  field  distribution  in 
the  applicator).  Penetration  depths (E  decreased  to  1/e  of  its 
value  at  the  surface)  at  2450  MHz,  for  the  operating  conditions 
encountered,  aie  greater  than  3.5  cm,  as  opposed  to  the  maximum 
half  sample  thickness  used  of  1.27  cm. 

With  these  assumptions,  the  power  conversion  formula  reduces 
to, 

03  = K E^  (3.5) 

where  the  dissipation  coefficient  K is  defined  as, 

K = TT  f€^€"/(€'^  + €"^)  (3.6) 

The  effect  of  the  temperature  on  the  dissipation  coefficient 
for  both  the  frozen  and  the  dried  regions  is  of  special  importance. 
According  to  data  by  Kan  & Yeat'^n  (1961)  as  the  temperature 
increases,  K can  increase  by  a factor  of  4 for  the  frozen  meat 
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b«tw««n  -40*C  and  0”C,  and  by  a factor  of  2 for  tba  drlad  noat 
batwaan  0*C  and  50*C.  Tba  incraaaa  in  K valua  with  incraaiing 
temparature  bas  an  acceleration  effect  on  the  power  absorption 
by  the  dielectric.  This  may  lead  to  a situation  known  as  a 
"rutiaway"  condition. 

Finally,  it  is  desired  to  avoid  any  melting  of  the  frozen 
region  in  the  course  of  the  dehydration  because  of  the  very  high 
loss  factor  of  the  wet  meat  which  may  trigger  the  runaway  situation 
previously  mentioned.  Therefore,  the  constraint  on  the  temperature 
in  the  frozen  core  will  be  that  no  melting  occurs  at  any  time. 

5.  Mass  transfer 

Although  operation  under  atmospheric  pressure  is  possible 
(Hohner  (1970)),  freeze-drying  is  generally  conducted  u ider  vacuum 
(.05  to  1 mmHg)  in  order  to  decrease  the  resistance  to  the  flow 
of  water  vapor  through  the  dried  layer.  This  type  of  operation 
is  referred  to  as  vacuum  freeze -drying.  The  higher  mass  transfer 
rate  reached  under  vacuum  yields  lower  temperatures  of  the  frozen 
core.  This  permits  a higher  heat  input  to  the  sample  and  leads 
to  shorter  drying  times.  The  vacuum  operation  is  practically 
necessary  for  microwave  f reeze-drying  as  short  drying  times  are 
normally  required  for  the  effective  use  of  microwave  energy. 

At  the  operating  pressure  range  mentioned  above,  the  mean 
free  path  of  the  molecules  inside  the  pores  of  the  dried  layer  is 
of  the  order  of  magnitude  of  the  mean  pore  diameter  (100  (i  as 
observed  by  Harper  (1962)  for  beef  meat).  Thus  the  flow  regime 
corresponds  to  the  transition  region  between  molecular  flow  and 
slip  flow.  Table  3.1  shows  the  values  of  the  Knudsen  number  at 
various  pressures  for  dried  beef  meat. 

The  variations  of  the  thermal  conductivity  of  the  dried 
layer  with  pressure  (cf.  Fig.  6.4)  supports  the  fact  that  the 
flow  is  in  the  transition  regime.  Fig.  6.4  shows  that  the  effective 
thermal  conductivity  is  constant  at  pressures  below  .01  mmHg 
and  corresponds  to  that  of  the  solid  matrix.  As  the  pressure 
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Increases  the  effective  thermal  conductivity  increases.  When 
the  pressure  reaches  about  10  mmHg,  the  effective  thermal  conductiv- 
ity starts  to  level  off  and  becomes  indepeiident  of  the  pressure 
again.  This  latter  stage  corresponds  to  the  beginning  of  the 
continuous  flow  regime  for  which  the  kinetic  theory  of  gases 
predicts  that  the  thermal  conductivity  does  not  depend  on  the 
pressu.re . 


Table  3.1  - Variations  of  Knudsen  Number  as  a 
function  of  pressure  (Temperature:  38 “C) 

Mean  free  Knudsen 


Pressure  path  Number® 

mmHg  ^ 

760  .065 

10  5 .05 

1 50  .5 

.1  500  5 

.05  1000  10 


a - calculated  assuming  a pore  diameter  of  100  n and  that 
the  mean  free  path  varies  inversely  with  the  pressure 


Rather  complex  equations  have  been  proposed  to  describe 
the  gas  flow  in  the  transition  region  (Wakao  et  al.  (1965)). 
However,  previous  studies  particularly  by  Sandall  et  al.  (1967) 
and  Dyer  & Sunderland  (1966)  show  that  it  is  possible  to  describe 
the  flow  by  Pick's  Equation  using  an  effective  diffusivity. 

Values  of  the  effective  diffusivity  have  been  determined  experi- 
mentally (Harper  (1962),  Dyer  & Sunderland  (1966),  see  Chapter 
VI,  Sec.  4).  Thus,  Pick's  Equation  is  used  in  this  study. 

6 . Moisture  sorption  in  the  dried  layer 

As  the  water  vapor  from  the  sublimation  flows  through  the 
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dried  layer  a finite  amount  of  moisture  is  expected  to  be  adsorbed 
by  the  dried  meat.  Sorption  isotherms  have  been  obtained  for 
beef  meat  by  Saravacos  & Stinchfield  (1965).  However,  calculations 
by  Sandall  (1966)  for  typical  conditions  of  freeze-drying  of 
turkey  breast  have  shown  that  the  average  sorbed  moisture  content 
is  approximately  3%  of  the  initial  moisture  content  or  1.5% 
moisture  by  weight.  These  calculations  were  based  on  adsorption 
data  of  King  et  al.  (1968)  and  assumed  an  external  heating 
through  the  dried  layer.  This  surprisingly  low  number  has  been 
attributed  to  the  fact  that  rather  high  temperatures  and  low 
pressures  exist  in  the  pores  of  the  dried  layer  A , ■'rger  amount 
of  adsorption  may  be  expectod  for  microwave  freeze-drying  as  the 
temperature  through  the  dried  layer  is  somewhat  lower.  However, 
our  own  experimental  observations  show  that  the  total  moisture 
remaining  in  the  dried  layer  once  the  ice  core  has  been  removed 
is  less  than  5%  of  the  initial  total  weight  of  the  sample. 

On  the  basis  of  these  observations  and  in  order  to  simplify 
the  model  it  seems  reasonable  to  neglect  the  adsorption  effect 
as  the  drying  is  generally  limited  to  the  removal  of  the  ice. 
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- ChaptBx  IV 


THE  MATHEMATICAL  MODEL 


A one-dimensional  model  of  the  freeze-drying  process  using 
microwave  dielectric  heating  is  derived  with  the  following  assumptions 
which  were  discussed  in  detail  in  Chapter  III : 

i - The  sublimation  zone  has  a zero  thickness  and  retreats 
uniformly. 

ii  - An  effective  diffusivity  can  be  used  to  describe  the  vapor 
flow  through  the  dried  layer  by  Pick's  Equation, 
iii  - The  p>ower  generated  by  the  microwaves  in  a unit  volume 
can  be  calculated  by  Eq&  (3.5)  and  (3.6)  from  the  field 
strength  in  the  dielectric  (if  one  assumes  a polarization 
of  the  electric  field  normal  to  the  slab  surface), 
iv  - The  moisture  sorption  effect  in  the  dried  layer  can  be 
neglected. 

The  material  to  be  freeze-dried  is  assumed  to  have  the  geometry 
of  an  infinite  slab  as  shown  in  Fig.  4.1.  The  left  edge  of  the 
slab  is  insulated  (or  center  line  of  a symmetric  slab)  while  the 
right  edge  is  exposed  to  a vacuum  at  temperature  T and  water  vapor 
concentration  Cp. 


1 . The  transport  equations 

Application  of  the  principles  of  conservation  of  mass  and 
energy  in  the  dried  layer  gives 


aw  _ 

3x  ^ 3t 


(4.1) 


at  “d 


(4.2) 
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Figure  4.1  - MndeJ  with  expected  forms  of  the  concentration 

and  temperature  profiles  for  the  microwave  freeze 
drying  process. 


wh«r«  the  concentration  C is  that  of  the  water  vapor  in  the  gas 
phase  inside  the  pores.  All  other  values  are  effective  quantities. 
In  !;he  frozen  region,  one  has 


9Tp 

at  “f 


(4.3) 


Assume  that  Fick's  Equation  can  be  used  to  describe  the  mass  trans- 
fer in  the  dried  layer.  Then 

W » - D |§  (4.4) 


Assume  that  the  heat  transfer  in  the  dried  layer  is  governed  by 
conduction  and  convection.  Then 


!!s 

dx 


+ C W (T_ 
pw  D 


T .c) 

ref ' 


(4.5) 


where  is  the  reference  temperature  for  enthalpy.  Assume  that 

heat  transfer  in  the  frozen  region  is  due  to  conduction.  Then 


®F 


3T, 


F 3x 


(4.6) 


Substitution  of  Equations  (4.4),  (4.5),  and  (4.6)  into  Equations 
(4.1),  (4.2),  and  (4.3)  gives 


i 92 

22 

- 2 22  = 0 
D 3t 

3x^ 

D 3x 

3x 

C W.  ^ 

DW  1 

i !l2- 

3x" 

3x 

3x 

“d 

' ‘'d 

tiL 

3kp 

Hz  1. 

3x^ 

3x 

3x  ■ Op 

3t 

’ ^F 

The  variation  of  the  mass  flux  W with  location  has  been 
neglected  in  the  convection  term  of  the  heat  transfer  equation  (4.8). 
It  has  been  shown  that  the  convective  term  is  small  compared  to 
the  other  terms  in  the  energy  equation  (Dyer  & Sunderland  (1968)). 

It  also  has  been  shown  that  the  mass  flux  does  not  vary  significantly 
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through  the  dried  layer  (Ma  i Peltre  (1973)). 


2.  The  boundary  conditions 

At  the  insulated  surface  (or  center  line)  we  have 


X = 0 


ax 


= 0 


for  t \ 0 


At  the  outer  surface  of  the  dried  layer,  we  have 

^9T, 


X = L 


■ *'d( ax)i 


for  t S.  0 


(4.10) 


(4.11) 


Due  to  the  low  ambient  pressure,  the  resistance  to  the  mass  trans 
fer  at  the  outer  surface  can  be  neglected.  Thus 

for  t"\  0 


X = L 


C(L,t)  = C, 


(4.12) 


At  the  interface  (ice-front)  a thermodynamic  equilibrium  between 
ice  and  vapor  is  assumed.  Thus,  the  concentration  of  the  water 
vapor  can  be  related  to  the  ice  temperature  by  the  equilibrium 
relationship,  or  mathematically  by 

X = X(t)  C.  = f,,(T.)  for  t 0 (4.13) 

1 Cl 

From  Pick's  Equation,  we  have 

(4.14) 

An  energy  balance  at  the  interface  also  gives 

« = X(t)  fort^O  (4.15) 

There  should  not  be  any  discontinuity  in  temperature  at  the  inter 
face.  Thus 


X = X(t)  T.  = Tjj  = T. 


for  t \ 0 


(4.16) 


Finally,  the  .ate  of  movement  of  the  interface  position  can 
be  related  to  the  t \te  of  sublimation  by 


- dX 


(4.17) 
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3.  Initio; 

As  ths  thlcknsss  of  th«  dried  layer  is  initially  aero,  it  is 
not  possible  to  define  the  initial  temperature  and  concentration 
profiles  in  the  dried  layer.  In  order  to  facilitate  the  solution, 
an  arbitrary,  but  small  initial  thickness  8^  of  the  dried  layer  is 

chosen.  An  initial  thickness  of  .1%  of  the  total  slab  thickness  L* 
is  used  in  the  numerical  solution.  Thus,  in  the  frozen  region,  we 
have 

t = o'*’  ■ ^o  for  0 / X / L-8^  (4.18) 


In  the  dried  layer,  the  initial  temperature  profile  is  assumed  linear 
and  is  given  by 


+ 

t = 0 


for  L-8  Z X Z L 
o - - 


(4.19) 


where  the  surface  temperature  is  given  by  the  boundary  condition 
at  the  outer  surface,  Eq.  (4.11). 

The  initial  concentration  profile,  which  is  also  assumed  to 
be  linear,  is  determined  from  the  initial  rate  of  sublimation  as 
no  equilibrium  is  established  yet  between  the  water  vapor  and  ice 
at  the  interface. 

w. 

t = o'"  C(x,o'")  = (L  - X)  for  L-8^  Z x / L 

(4.20) 

The  initial  rate  of  sublimation  W,  is  given  by  the  energy  balance 
at  the  interface,  Eq.  (4.15). 


4.  Fixed  boundary  transformation  and  normalization 

Equations  (4.7),  (4.8),  and  (4.9)  are  a set  of  partial  differ- 
ential equations  which  are  to  be  solved  with  a moving  boundary. 

In  order  to  facilitate  the  solution  of  these  equations,  Landau's 
transformation  (Bankoff  (1964))  is  used  to  convert  the  problem 
of  a moving  boundary  to  that  of  a fixed  boundary.  The  following 
transformations 
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z 


L - X 
L - X(t) 


(4.21) 


y ® 


X 

X(t) 


(4.22) 


are  used  in  conjunction  with  the  following  normalized  variables 


and 


S a 


(4.23) 

(4.24) 

(4.25) 

(4.26) 

(4.27) 

(4.28) 

(4.29) 


Substitution  of  the  transformation  and  normalization  equations 
into  Equations  (4.7),  (4.8),  and  (4.9)  gives 


Oz' 


D*  3z  " D*  dtr 


Ml  Xm  C.2  gr 

Sz  ■ gr  ■ ° 


(4.30) 


2 

g e. 
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r 1 

/ c p a , . 
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f pw  ,1  \ 

i] 

(c*  gz 

^ PdS“d  <^d  I 

PpS 

k°  (Tj^  - T^) 


(4.31) 
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(4. 


3y 


3y  a*  V d«j  3y  a*  ar  a*  . o 


° l‘ 


dcj  3y 


a*  o 

F'  R o' 


The  boundary  and  initial  conditions  become 
30„ 


y = 0 


2 = 0 


3y 


= 0 


for  0 


= 1 


'D  h L (1-S)  32 

^ = 0 for  t!'^  0 

At  the  interface  where  y = 2 = 1,  one  has 
Tf  = fi(®i)  for  tr\0 


for  0 


ffo 

s ay  PpCp  1 -s  32 


L m w. 

s 1 

‘'f  (^R-^o) 


D F i 


for  *t\  0 


(4.33) 

(4.34) 

(4.35) 

(4.36) 
f or  tr^O  (4.37) 

(4.38) 


and  the  movement  of  the  interface  is  given  by 


dr 


W.  L 
1 

P ffao 


where  the  rate  of  sublimation  can  be  represented  by 


w = 
i L(l-S) 


f— ^ 

\32/z=1 


The  initial  conditions  become 

©p  = 0 for  0 Z y / 1 


«D  = U-  *)  Qs.o 


(4.39) 


(4.40) 

(4.41) 

(4.42) 


for  0 Z 2 Z 1 

where  © is  given  by  the  boundary  condition  at  the  outer  surface, 
d I o ^ 

_ h L (1-S) 

®S,o  - ilTT-Sj.lig  for  I = 0 (4.43) 


Finally,  one  has 


r * 


L(1-S)W, 


i^_o 


■>o<Sp-S> 


for  0 Z * ^ 1 (4.44) 


- Chaptsr  V - 


THE  NUMERICAL  METHOD 


Under  the  assumptions  made  earlier,  the  set  of  partial 
differential  Equations  (4.30),  (4.31),  and  (4.32)  with  the  initial 
and  boundary  conditions,  Equations  (4.33)  through  (4.44),  constitutes 
a one-dimensional  mathematical  model  of  the  microwave  freeze- 
drying process.  Solutions  of  these  equations  will  provide  one 
with  a means  to  predict  the  time  variations  of  the  dependent 
variables  such  as  temperature  and  conce.*tration  whithin  the 
sample  in  the  course  of  the  dehydration  and  obtain  some  information 
on  drying  rate  and  drying  time,  necessary  for  design  purposes. 

The  partial  differential  equations  to  be  solved  form  a set 
of  three  parabolic  equations  with  phase  change  at  a moving  boundary. 
Even  in  the  simpler  case  of  one  equation,  it  is  known  matrematically 
as  a very  cumbersome  problem  and  a general  solution  has  been 
found  only  in  special  instances  (Bankoff  (1964)),  such  as  that  of 
a serai-infinite  slab,  often  referred  to  as  the  Von  Neumann  problem. 

No  general  solution  has  been  found  yet  for  the  finite  slab  with 
different  properties  on  both  sides  of  the  moving  boundary.  In 
the  case  under  consideration,  there  are  three  equations  which  are 
coupled  through  their  boundary  conditions  at  the  moving  boundary 
(interface).  Moreover,  all  physical  properties  of  the  sample 
are  considered  as  a function  of  the  dependent  variables  (temperature 
and  water  vap>or  concentration).  Thus,  they  are  considered  to 
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vary  with  both  time  and  location.  Consaquantly,  an  analytical 
or  even  a semi -analytical  solution  are  discarded  as  most  unlikely, 
and  a numerical  technique  is  used  to  obtain  the  solutions. 

1 . The  finite-difference  method  of  approximation 

The  fixed  boundary  domain  of  Chapter  IV,  Sec.  4,  is  divided 
into  intervals  as  shown  in  Fig.  5,1,  The  space-axis  is  divided 
into  M intervals  of  equal  length  Ly  in  the  frozen  region  and  K 
intervals  of  length  Lz  in  the  dried  layer.  The  time  axis  is  divided 
in  intervals  of  duration  tX. 


INTERFACE 


VACUUM 

t^n+1 


Figure  5.1  - Numerical  grid  in  the  fixed  boundary  domain 


The  size  of  a'C  is  chosen  so  that  the  numerical  solution 
is  stable  and  accurate  with  reasonable  computation  time  (see 
discussion,  Sec.  4).  The  numerical  grid  thus  obtained  is  used 
to  derive  finite-difference  equations  from  the  partial  differential 
equations . 

For  this  purpose  the  partial  differential  equations  are  approx- 
imated by  the  Crank -Nicol son  method  which  offers  the  best  approx- 
imation for  tne  time -derivatives  and  thus  the  smallest  truncation 
error  (Smith  (1965)). 

The  general  technique  used  in  the  derivation  of  the  finite- 
difference  equations  is  described  here  and  the  reader  is  referred 
to  Appendix  A for  the  complete  equations  whl.>'  are  rather  voluminous 
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due  to  the  complexity  of  the  mathematical  system. 

If  one  uses  U as  a dummy  variable,  then  according  to  the 
Crank-Nlcolson  method,  the  following  approximations  are  made: 
with 


“k.„  = 

one  has 


(5.1) 


(5.2) 


I 


(Vl,nn'"k-l,n.l  ^ Vl,n'"k- 


.^2 

(• 


,ax‘^/k,n+% 


2 Ax 
2U. 


2 Ax 


(5.3) 


k.n-i-l'^  \-L.n+l  + \-t-l.n~  \-l.nj 


Ax 


Ax 

(5.4) 


(-) 

\9wk,  n+% 


\.n^l  -\.n 

Ar 


(5.5) 


where  Xj^  stands  for  or  yj^.  This  procedure  is  used  to  derive 
all  the  finite-difference  equations.  Typically,  an  equation  such 
as 


— - + a(x.'C)  — + b(x,r)  ^ - c(x,'C) 

9x 


(5.6) 


is  approximated  by 


\-t-l.n-t-l~^\,n-*-l'*'  \-l,n^-l  + a /^k-t-1  .ni-l~  \-l,n+l  + b.  ..\.n-t-l~  ^k.n 
o 2Ax 


Ax 


ATT 


+^k+l.n  ^'^k.n^  ^k-l.n  + a,  ^k+l.n  ^k-l.n  + b,  ^k.n+l"  ^k.n  = c.  . , + c. 

2^ " Ar 


(5.7) 

It  should  be  noted  that  the  means  of  the  products  are  evaluated  as 
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^“•''^n+Jj  " ^ ^“n+l''n+l 


n n‘ 


(S.8) 


2 . Method  of  solutioi 

Although  it  is  clear  that  the  coefficients  in  Equation  (5.6) 
are  time  dependent  and  therefore  unknown  at  time  f , let  us 
suppose  for  simplicity  that  they  are  indeed  known,  "xheir 
determination  will  be  discussed  in  the  next  section.  Equation 
(5.7)  can  then  be  rearranged  by  grouping  all  known  terms  on  the 
right  hand  side  of  the  equal  sign,  to  give  the  general  form 


""k-l  ^ <^k  \ ^ U, 


k '^k+l  “ ^k 


k = 2 K 


(5.9) 


Eq.  (5.9)  actually  represents  a set  of  K-l  linear  equations  at 
K-1  grid  points  between  the  boundaries:  x^,  x^ x,^.  Two 

more  relations  are  obtained  from  the  boundary  conditions  which 
are  written  in  a general  manner: 

at  the  outer  edge  (center -line  or  external  surface) 


"a  - \ 


and  at  the  interface 

P U + U an 

K ^K+1 


at  X » 


*1 


at  X = X 


“K+l 


(5.10) 


(5.11) 


Details  on  the  derivation  of  these  relations  as  well  as  their 
actual  complete  form  are  given  in  Appendix  A. 

Equations  (5.9),  (S.IO).  and  (5.11)  form  a set  of  K+l  linear 
equations  which  are  to  be  solved  with  K+l  unknowns  Uj^,  with 
k **  1,...,K+1.  The  set  is  solved  in  a sequential  manner  by 
elimination  (Gauss's  elimination  technique  - Smith  (1965)).  A 
forward  elimination  reduces  Equation  (5.9)  to  the  form 

A 


u , + S u 

k '^k  k+l 


k = 1 K 


(5.12) 


where  the  S's  and  A's  are  kncwn.  Equations  (5.11)  and  (5.13) 
are  solved  simultaneously  for  The  U^'s  are  then  determined 

by  Equation  (5.12)  through  back  subsitution. 
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Details  of  the  algorithm  are  given  In  Appendix  B. 


3 . The  general  numerical  scheme  of  solution 

Before  the  equations  can  be  solved  to  obtain  the  temperature 
and  concentration  profiles  at  method  outlined  in 

the  previous  section,  the  current  position  of  the  interface  and 
values  of  the  properties,  (i.e.,  thermal  diffusivity,  diffusivity 
and  dissipation  coefficient)  must  be  evaluated.  As  they  are  functions 
of  the  dependent  variables  (temperature  and  concentration),  their 
values  have  to  be  assumed  on  the  basis  of  their  previous  history 
and  iterated  if  necessary. 

The  interface  motion  is  related  theoretically  to  the  rate 
of  sublimation,  or  the  rate  at  which  energy  is  consumed  at  the 
interface,  by  Equation  (4.39) 


dr 


L W. 

1 

o 

P crap 


This  equation  is  used  to  calculate  the  position  of  the  interface 

%.*!• 


s = s - at:  L w. 

n+1  n i,n+H 

o ’ ^ 

P tfap 


(5.14) 


The  rate  of  sublimation 

rate  at  “Cr  and  earlier, 
n 


at  is  not  known.  One  only  knows 

n+1  ' 

Thus,  is  approximated  by 


the 


W.  , 

i,n+^ 


W.  + (w.  - W.  ) = W.  , + 

i,n-^  i,n  i,n-l 


(AW. ) (5.15) 

1 n 


From  the  profiles  calculated  at  C"  t ^ value  of  W.  , is  obtained 

n+1’  i,n+^ 

by  an  energy  balance  at  the  interface  (see  Appendix  A for  details 
of  the  calculations)  and  used  to  recalculate  the  interface  position 
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This  process  is  repeated  until  the  variation  of 
is  only  on  the  7th  digit.  A fixed  number  of  two  iterations  was 
found  to  be  5$ufficient  in  most  cases. 

The  proper tie.s.  i.e.,  diffusivity,  thermal  diffusivity  and 
diss  ipati  on  coefficient,  vary  w'tlj  both  time  and  location.  They 
are  functions  of  the  dependent  variables  (temperature  and  concen- 
tration). They  are  calculated  from  their  current  value  through 
functional  relationships  derived  from  literature  data  (see 
Chapver  VI).  Consequently,  as  their  exact  value  cannot  be  known 
unless  the  dependent  variables  are  estimated,  their  \'alue  at 
is  approximated  by  that  and  iterated  if  need  be. 

In  fact,  the  physical  properties  are  not  strong  functions  of  the 
temperature  and  concentration.  Thus,  they  are  not  expected  to 
vary  signif  icantly  over  one  t icne -increment , considering  the  small 
si*e  of  the  time-increment  used  (At  ^ 7.2  sec  as  compared  to  a 
total  drying  time  of  2 to  10  hours  (7,200  to  36,000  sec)).  This 
approximation  proved  to  be  satisfactory  for  the  thermal  diffusivity 
and  the  dissipation  coefficient. 

However,  the  approximation  is  not  sufficient  for  the  di;tfus- 

ivity,  as  small  fluctuations  on  D*  are  magnified  in  the  term 
3D* 

in  Equation  (4.30),  which  results  in  high  fluctuations  on 
the  concentration.  It  is  necessary  to  obtain  a better  accuracy 
of  D*  to  eliminate  the  concentration  fluctuations.  This  was 
done  by  inserting  the  calculation  of  the  diffusivity  inside 
the  iteration  loop  for  the  interface  position,  so  that  the 
value  of  D*  is  updated  once  the  temperature  and  concentration 
profiles  have  been  calculated  at 

This  problem  does  not  occur  with  the  thermal  diffusivity 

for  the  time-increment  used,  although  the  energy  equations 

det*  9(1* 

comprise  a similar  term  D or  F.  This  can  be  explained  by  the 

fact  that  a*  and  a*,  havi  a mulh^  smaller  value  than  that  of  D*. 

D r 

Fig.  5.2  gives  the  flow  chart  of  the  general  numerical  scheme. 

The  boundary  conditions  at  the  interface  are  such  that, 
once  the  interface  position  and  the  properties  have  been  estimated. 
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a second  iteration  loop  is  needed  to  calculate  the  profiles. 

In  effect,  as  the  heat  and  mass  transfer  equations  are  coupled 
at  the  interface,  no  explicit  form  of  the  boundary  condition  - 
such  as  Eq.  (5.11)  - can  be  obtained  for  these  equations  ((4.30), 
(4.31),  and  (4.32)).  Thus,  the  temperature  at  the  interface 
is  assumed  on  the  basis  of  its  previous  history, 


T.  ^ = T.  + (T. 

i,n+l  i,n  1, 


T.  , ) = T. 

i,n-l  1, 


(AT.)  , 


(5.16) 


and  then  iterated.  From  the  assumed  interface  temperature,  the 
concentration  at  the  interface  can  also  be  calculated  through 
the  equilibrium  relationship.  Thus,  the  temperature  and  con- 
centration profiles  can  be  calculated.  A theoretical 

value  is  then  computed  from  a full  energy  balance  at  the  interface 
and  compared  to  the  assumed  value.  Unless  an  agreement  on  the 
6th  digit  is  reached,  a new  value  is  assumed  for  the  interface 
temperature  T^  according  to 

+ (1  ->^)  Tr^  (5.17) 

and  the  iteration  proceeds.  The  coefficient  ^ in  Equation  (5.17) 
is  used  to  force  and  accelerate  the  convergence.  It  is  constantly 
updated,  by  the  scheme  described  and  justified  in  Appendix  B, 

With  this  scheme,  convergence  is  generally  attained  in  2 to  4 
iterations.  The  general  flow  chart  of  the  algorithm  is  described 
in  Fig.  5.2. 


4.  Stability,  accuracy  of  the  solution,  and  simulation  speed 

The  stability  and  the  accuracy  of  the  solution,  as  well  as 
the  simulation  speed  are  essentially  affected  by  the  size  of  the 
space  and  time  increments.  In  simple  cases,  a Fourier  analysis 
of  the  mathematical  system  can  provide  a stability  criterion  in 
terms  of  the  ratio  . It  often  results  in  a maximum  value  for 
this  ratio,  beyond  which  the  solution  becomes  unstable  (i.e., 
the  numerical  solution  oscillates  about  the  true  solution). 
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Unf ortunataly,  due  to  the  high  degree  of  complexity  of  the 
nathematical  ayeten  encountered  here,  such  an  analysis  is  expected 
to  be  extremely  cunberaone.  Thus,  variations  of  are  used 
to  study  the  behavior  of  the  solution. 

The  space  increment  is  usually  fixed  and  its  value  decided 
on  the  ground  of  the  accuracy  desired  in  the  determination  of 
the  profiles  of  the  dependent  variables.  When  the  profiles  are 
smooth  lurves  with  no  sharp  slope  changes,  20  space  intervals  are 
sufficient  to  give  a good  description  of  the  curves.  Inis  is  done 
here.  Both  the  frozen  region  and  the  dried  layer  are  divided 
into  20  equal  intervals  (Az  = Ay  = 0.05).  As  the  time-variation 
of  the  variables  is  small,  one  may  think  that  a large  value  of 
the  time  increment  At  can  be  used  without  an  appreciable  loss  of 
accuracy.  For  instance,  for  a total  drying  time  of  6 hours 
(216,000  sec),  one  might  be  able  to  use  a At  as  high  as  3.6  min 
(216  sec)  with  a reasonable  accuracy,  and  thus  a high  simulation 
speed.  However,  the  stability  requirement  is  such  that  much  smaller 
values  of  At  had  to  be  used,  especially  at  the  beginning  of  the 
simulation.  Table  5.1  shows  the  results  of  the  stability  study 
for  different  At,  as  a function  of  time  for  a typical  calculation. 


Table  5.1  - stability  of  the  numerical  solution 
Time  increment  Value  of  the  time  increment  At 

number  0.072  sec  0.72  sec  7.2  sec 

0 - 100  stable  unstable  unstable 

100  - 500  stable  stable  unstable 

500  - 00  stable  stable  stable 


Hence,  a At  of  .072  sec  is  used  for  the  first  100  time  increments. 
A At  of  .72  sec  is  used  for  the  next  400  time  increments  and  a At 
of  7.2  sec  was  used  for  the  rest  of  the  computation.  The  corres- 
ponding, normalized  time  increment  A'C  is  computed  from  the  norm- 
alization equation  (4.23).  The  total  number  of  time  increments 
necessary  to  reach  the  end  point  (S  * 0.001)  varied  from  1,500 
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to  8,000  depandlng  on  th«  total  drying  tlmt.  Th«  average  total 
number  of  iterations  per  time  increment  is  of  the  order  of 
6 to  10.  The  computer  time  required  is  of  the  order  of  3 to  5 min 
of  CPU  time  per  hour  of  total  drying  time  on  a PDP-10  digital 
computer . 
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- PART  II 


APPLICATION  TO  THE  FREEZE-DRYING  OR  BEEF  MEAT  AT  2450  MHz 


Chapter  VI:  Numerical  Data  Used  in  the  Calculations 


Chapter  VII:  Theoretical  Results  and  Discussion 
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- Chapter  VI  - 


NUMERICAL  DATA  USED  IN  THE  CALCULATIONS 


This  chapter  is  concerned  with  the  presentation  of  the 
numerical  values  of  the  physical  constants  and  the  functional 
relationships  of  the  non-constant  properties  used  in  the  mathe- 
mathical  simulation  of  the  f ree2e-drying  of  a one-dimensional 
beef  meat  slab  using  microwave  dielectric  heating  at  2450  MHz. 

The  slab  is  cut  normal  to  the  fibers. 

1 . The  constants 

The  numerical  values  used  for  the  physical  constants  are 

summarized  in  Table  6,1.  Some  were  assumed  because  of  the  lack 

of  experimental  data  in  the  literature.  This  was  the  case  for 

the  heat  transfer  coefficient.  Although  vacuum  freeze-drying 

has  been  on  the  industrial  scene  for  many  years,  no  value  of  the 

heat  transfer  coefficient  could  be  found.  This  is  apparently 

due  to  the  general  feeling  that  it  is  not  a significant  parameter 

and,  should  be  determined  according  to  a particular  design.  The 

value  given  in  Table  6.1  was  estimated  from  the  initial  slope 

of  a preliminary  experimental  drying  curve.  Introduction  of  a 

continuous  recording  of  the  sample  weight  during  drying  permitted 

a more  accurate  estimation  of  the  heat  transfer  coefficient 

from  the  start-up  portion  of  the  drying  curve  (see  Chapter  IX). 

-4  2 

It  was  then  found  that  an  h value  of  the  order  of  2x10  cal/cm  /sec/®C 
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Table  6.1  - Numerical  values  of  the  constants  used 
in  the  simulation. 


Constant 

Numerical  value 
used 

Origin 

h 

5x10  ^cal/cm^/ 
sec/*C 

assumed 

P 

.92  g/cm^ 

CRC  Hdbk 

"d 

.32  g/cm^ 

measured 

. 96  g/cm^ 

measured 

C 

pw 

.5  cal/g/»C 

assumed 

s 

.36  cal/g/»C 

Awberry  & Griffiths  (1933) 

s 

.43  cal/g/*C 

Rey  (1964) 

AH 

s 

675  cal/g 

Hohner  (1970) 

after  Threlkeld  (1962) 
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would  be  more  suitable  for  the  freeze-dryer  used  in  the  experiments 
This  value  is  used  for  the  comparison  of  the  experimental  drying 
curves  with  those  predicted. 

It  was  also  assumed  that  the  specific  heat  of  the  water 
vapor  at  the  usual  freeze-drying  operating  conditions  (-50®C  to 
60°C;  .05  to  1 mmHg)  could  be  approximated  by  that  at  100®C  and 
1 atm.  (exactly  .4836  cal/g/®C  after  Lange's  Handbook  of  Chemistry) 
It  may  be  somewhat  different,  due  to  the  difference  in  temperature 
and  pressure  as  well  as  the  relative  importance  of  the  wall  effects 
in  the  transition  flow  regime. 

The  average  porosity  (void  fraction)  is  computed  from  the 
initial  moisture  content  and  the  densities  of  pure  ice  and  frozen 
meat , 


^ \ ^ (6-1) 

0 

The  shrinkage  effect  occuring  on  the  dried  product  is  neglected 
and  the  same  porosity  value  is  used  for  the  frozen  region  and 
the  dried  layer.  It  is  taken  as  that  of  the  frozen  meat. 

2 . Effective  thermal  conductivity  of  frozen  beef 

Numerical  data  from  various  literature  sources  were  gathered 
for  heat  flow  parallel  to  the  fibers.  An  attempt  was  made  to 
fit  the  data  with  a functional  relationship  of  the  following 
form , 

kp  = a + h,^/T^-~F  (6.2) 

This  form  is  chosen  because  a plot  of  kp  vs  temperature  resembles 
a parabola  with  axis  parallel  to  the  horizontal  temperature 
axis,  and  tangent  at  the  vertex  parallel  to  the  vertical  axis  at 
about  □‘’C.  A rather  good  fit,  as  shown  in  Fig.  6.1  was  obtained. 
The  parameter  T^  is  adjusted  in  order  to  reach  the  best  fit, 
and  could  correspond  to  to  the  average  fusion  temperature  of  the 
frozen  meat. 


Some  further  interesting  correlationc  were  found.  A plot 
of  a in  Equation  (6.2)  vs  fat  content  shows  that  a varies  linearly 
with  the  fat  content  while  it  does  not  seem  to  vary  significantly 
with  the  water  content.  On  the  contrary,  Fig.  6.3,  b does  not 
seem  to  depend  on  the  fat  content  but  varies  linearly  with  the 
water  content.  It  is  found  that,  approximately 


a 

b 


(2.95  - 25  x^)  X lO"^, 
[02  + 1.8  (x^  - .65) 


cal/sec/cm/°K  (6.3) 

X 10  cal/sec/cm/( ®K) ^ ^ 

(6.4) 


and  on  the  average, 
= 271. 5*K 


(6.5) 


Equation  (6.2),  along  wi.h  (6.3),  (6.4),  and  (6.5),  is  used  to 
determine  kp  in  the  simulation,  with  the  temperature  T expressed 
in  ‘’K.  a and  b are  determined  from  the  fat  and  water  content  of 
a given  choice  of  meat  by  Equations  (6.3)  and  (6.4). 


3 . Effective  thermal  conductivity  of  dried  beef 

The  experimental  data,  obtained  at  various  pressures  by 
Harper  (1962),  are  used  to  determine  a functional  relationship 
between  the  thermal  conductivity  and  the  total  pressure.  For 
this  purpose  a 4th  degree  polynomial  is  least  square  fit  to 
Harper's  data  in  terms  of  the  natural  logarithm  of  the  pressure 
expressed  in  mmHg  (see  Fig.  6.4).  The  following  expression  was 
obtained 

■-  10"^(1.126  + 0.1351  Log  P + .01565(Log  P)^ 

(6.6) 

- 3.187  X 10'^(Log  P)^  - 5.483  x lo"^(Log  P)"^ 

where  P is  in  mmHg  and  k^^  in  cal/sec/ciii/'C. 

Additional  data  have  been  obtained  by  other  investigators, 
particularly  by  Massey  & Sunderland  (1967).  However,  they  are 
mostly  the  results  of  calculations  for  typical  treeze-dr ying 
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Figure  6.1  - Thermal  conductivity  of  frozen  beef  vs 
for  various  m.oxsture  and  fat  contents. 


Symbol 

Literature  Source 

% moisture 

% fat 

▲ 

Lentz  (1961) 

75 

.9 

• 

Miller  & Sunderland 
(1963) 

69.5 

# 

Hill  et  al.  (1967) 

76.5 

2 . 35 

■ 

Hill  et  al.  (1967) 

78.7 

1.4 
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Thermal  Conductivity,  Cal/ sec/cm /®C 


vs  moi'^ture  content  of  frozen  lieef . 


Figure  6 . L’  - Variations  of  parameter  a in  Ecj . (6.2) 

vs  fat  content  of  frozen  beef. 
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•xparlments,  batsd  on  a quasi -steady  state  model. 


4.  g»tgtivf 

A simplified  two-parameter  model,  as  proposed  by  Sandall  et  al. 
(1967),  is  used  to  describe  the  variations  of  the  effective  dif- 
fusivlty  D with  the  pressure  in  the  transition  flow  regime, 


D = 


(6.7) 


’■ '.:^re  P is  the  pressure  in  atm,  D®  the  bulb  diffusivity  at  1 atm, 
and  the  !Inudsen  diffusivity. 

Although  the  description  of  the  water  vapor  mass  transfer 
in  vacuum  f reeae-drying  processes  by  Pick's  Equation  seems  widely 
accepted,  few  experimental  data  could  be  found  on  the  diffusivity 
for  dried  beef  in  the  transition  flow  regime.  Sandall  et  al . 
(1967)  diffusivity  measurements  for  D®  and  in  the  case  of 
turkey  breast  were  used  in  che  theoretical  study.  They  were 
corrected  to  the  expected  average  temperature  of  the  dried  layer 
(20*C),  assuming  that  D®  varies  as  T and  as 
Further  temperature  effects  were  neglected  in  the  model  as  the 
relative  variation  on  D is  less  than  20%.  So  that  the  following 
relationship  was  finally  used  to  calculate  the  diffusivity  as  a 
function  of  pressure. 


D = 


78.5 
3.4  + P 


(6.8) 


where  P is  in  mmHg  and  D in  cm  /sec.  This  was  done  as  an  approx- 
imation, due  to  the  lack  of  experimental  data  on  for  beef 
meat.  It  seems  to  be  a reasonable  approximation  at  low  pressures 
if  one  considers  the  rather  good  agreement  between  the  diffusivity 
calculated  at  1.44  mmHg,  from  Sandall *s  O'  and  values  (for 
turkey  breast)  by  Eq.  (4.7),  and  that  re^rted  by  Dyer  & Sunderland 
(1966)  for  beef  meat  at  the  same  pressure.  These  data  are  shown 


in  Table  6.2. 
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5.  Pijiip»tiQn  co«iiici«)t  of  izoMM  mnd  driad  ba«f 

Th«  dlalactrlc  properties  of  frosen  end  dried  beef  are  expected 
to  bo  a function  of  the  fat,  moisture  and  salt  contents  as  well 
as  the  temperature  and  the  frequency  of  the  microwave  field. 
Unfortunately,  only  the  temperature  effect  for  frozen  and  dried 
beef  is  available  at  3000  Ml^z  as  reported  by  Kan  and  Yeaton 
(1961).  Kan  and  Yeaton 's  data  were  assumed  to  hold  at  2450  MHz 
and  used  to  determine  the  dissipation  coefficient  K defined  by 
Equation  (3.6).  Figures  6.5  and  6.6  show  the  plots  of  K^^  and 

K_  vs  temperature.  The  analytical  expressions  reached  for  these 

^ 3 2 

polynomials  are  (K  is  in  cal/sec/cm  /(V/cm)  ; T in  “K) : 

Kjj  * 10'S(14.39  - .1577  T ♦ 5.577  x 10"^T^ 

(6.9) 

-5.924  X 10“^T^) 

Kj,  « 10"y(294.7  - 3.628  T ♦ 1.464  x 10‘^T^ 

(6.10) 

-1.922  X 10"V) 

6.  Ice-water  vapor  equilibrium  relationship 

The  equilibrium  data  of  pure  ice,  were  used  in  this  work 
although  a depletion  of  the  equilibrium  vapor  pressure  of  as 
much  as  20%  may  occur  in  frozen  beef  as  reported  by  Dyer  et  al . 

(1966).  The  equilibrium  vapor  pressure  at  the  ice  front  was 
related  to  the  concentration  of  water  vapor  in  the  pore  using 
the  perfect  gas  law.  Pig.  6.7  shows  the  plot  of  the  equilibrium 
water  vapor  concentration  and  pressure  in  the  pore  vs  ice  temperature. 
The  following  analytical  expression  is  used  in  the  simulation 
to  depict  the  ice -vapor  equilibrium, 

= exp( -63.23  + 0.2969  T^  - 4.038  x lo’SJ)  (6.11) 

where  C.  is  the  equilibrium  concentration  of  water  vapor  at  the 
^ 3 

ice-front  inside  the  pore,  in  g/cm  , and  is  the  ice  equilibrium 
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X 10 


^ Least  Square  Fit  Model 


Concentration,  g/cnn?  ^ 
Pressure,  mm  Hg  ■ 


Temperoture,  ®C 


temperature, *K. 
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- Chapter  VII  - 


THEORETICAL  RESULTS 


The  mathematical  movlel  derived  earlier  has  been  applied  to 
the  freeze-drying  of  beef  meat  with  microwave  energy  at  2450  MHz. 
The  sample  is  assumed  to  have  the  geometry  of  a one-dimensional 
symmetric  slab  with  drying  on  both  sides.  Using  the  previous 
numerical  data  (Chapter  VI),  calculations  were  carried  out  in  order 
to  study  the  effects  of  the  physical  constraints  of  the  process 
upon  the  total  drying  time  as  a function  of  the  process  variables. 
The  results  of  these  calculations  are  presented  and  discussed 
below. 

It  should  be  noted  that  the  following  results  are  obtained 
by  employing  a typical  set  of  numerical  data.  Care  must  be  exerted 
to  extend  the  present  results  to  other  "practical"  systems  as  it 
is  important  to  have  more  specific  and  comprehensive  experimental 
data.  These  results  are  primarily  intended  as  a verification 
of  the  response  of  the  mathematical  model  under  typical  operating 
conditions  to  serve  as  design  guidelines. 

1 . Typical  simulation  outcome 

It  is  possible  to  investigate  quantitatively  the  variations 
of  any  characteristic  properties  through  simulation.  Of  special 
interest  are  the  temperature  and  water  vapjor  concentration  profiles 
within  the  sample.  The  variations  of  the  sample  weight  (or  the 
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lc«-<front  poiltlon)  and  th«  charactarlitlc  ttmp«ratur«i  with  timt 
ar«  alto  of  intaratt  at  thay  raflact  the  ganr^ral  fraaaa-drylng 
parformanca  for  glvan  '^paratlng  condltlont. 

Of  Importanca,  too  It  tha  quaatlon  of  anargy  convartlon 
afflclancy  for  sublimation,  which  can  be  Illustrated  by  the  com- 
parative histories  of  the  microwave  power  absorbed  in  the  sample 
and  the  power  consumed  by  the  sublimation. 

The  operating  conditions  used  in  the  calculations  are  shown 
in  Table  7.1. 

Fig.  7.1  shows  typical  concentration  and  temperature  profiles 
during  the  freeze-drying  process. 

The  concentration  profiles  in  the  dried  layer  exhibit  a 
rather  small  curvature.  In  a previous  study  (Ma  & Peltre  (1973)) 
of  the  same  process  using  a simpler  model  in  which  the  pioperties 
were  assumed  uniform  and  constant,  the  concentration  profiles 
were  found  to  be  essentially  straight  lines.  Thus,  this  small 
curvature  may  be  attributed  to  the  variation  of  the  diffusivity 
with  location,  as  the  time  derivative  in  the  mass  balance 
equation  is  negligible  compared  to  the  space  derivatives,  due 
to  the  slow  motion  of  the  interface  and  high  value  of  the  diffus- 
ivity. 

On  the  other  hand,  the  temperature  profiles  in  both  the 
frozen  and  dried  regions  are  parabolic  in  shape  which  is  expected 
from  a quasi -steady  state  analysis.  The  locations  of  the  temperature 
maxima  in  both  the  frozen  and  dried  layers  have  specific  signi- 
ficance as  they  represent  the  maximum  temperature  in  these  regions 
at  any  time.  In  fact,  the  temperature  in  the  frozen  region  always 
should  be  kept  below  the  melting  point  of  the  frozen  material 
(around  -1.5®C)  while  an  upper  limit  of  60°C  is  often  quoted  for 
the  dried  layer  in  order  to  avoid  the  thermal  degradation  of  the 
dried  product.  As  shown  in  Fig.  7.1,  the  maximum  temperature 
is  at  the  insulated  wall  (center  line)  in  the  frozen  region 
while  in  the  dried  layer  it  is  located  at  the  external  surface 
during  the  major  portion  of  the  dehydration  process.  However, 
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10  X Concentration,  G/CM 


the  vertex  of  the  parabola  shifts  from  right  to  left  as  time  elapses. 
Near  the  end  of  the  process,  the  maximum  temperature  is  located 
inside  the  dried  layer.  Thus,  hot  spots  inside  the  dried  layer 
may  be  createa. 

Fig.  7.2  shows  the  maximum  temperatures  in  the  frozen  and 
dried  regions,  and  the  interface  temperature  and  position  as  a 
function  of  time. 

The  temperatures  drop  to  a minimum  at  the  beginning  of  the 
process,  during  the  start-up  stage  as  the  ice-front  temperature 
drops  to  the  "wet  bulb"  temperature  of  the  chamber.  This  may 
be  attributed  to  the  high  mass  transfer  rate  resulting  from  the 
small  initial  thickness  of  the  dried  layer,  as  indicated  on 
Fig.  7.3  by  the  high  initial  value  of  the  power  consumed  by 
sublimation.  The  power  consumed  by  sublimation  at  the  interface 
is  given  by, 

Q.  = AH  W.  (7.1) 

Table  7.1  - Assumed  operating  conditions 

E 125  V/cm 

Pp  . 29  tnmHg 

pw  .075  mmHg 

T 20*C 

K 

T -15°C 

o 

L 1.27  cm 

x .73 

w 

x^  . 009 

The  high  sublimation  rate  causes  the  loss  in  sensible  heat  of 
the  frozen  core. 

As  the  "wet  bulb"  temper?».ture  is  different  from  the  frost 
temperature  in  the  chamber,  a small  drying  rate  exists  due  to 
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heat  supplied  from  the  surrounding  atnosphere  to  the  sample. 

As  the  thickness  of  the  dried  layer  i.icreascs  due  to  the  subli- 
mation, the  mass  transfer  resistance  through  the  dried  layer 
increases.  Thus,  the  ice-front  pressure  and  temperature  increase. 
The  slow  deevease  of  the  power  consumed  by  the  sublimation 
(supplied  from  the  surrounding)  at  the  interface  in  Fig.  7.3 
reflects  the  augmentation  of  the  resistance  to  the  heat  transfer 
as  the  thickness  of  the  dried  layer  increases. 

As  the  microwave  power  is  turned  on  at  time  t = 30  min, 
and  as  the  mass  transfer  resistance  is  increased  due  to  the 
relatively  fast  receding  of  the  ice-front,  the  temperatures  start 
to  rise  to  a maximum  and  then  decrease  again.  As  the  majority 
of  the  microwave  power  is  absorbed  by  the  frozen  region,  it  will 
absorb  le  srgy  when  its  volume  is  reduced  as  shown  in  Fig. 

7.3.  Ths  s the  power  input  decreases  the  drying  rate  slows 
down.  The  .^ecrease  of  the  ice-front  temperature  indicates  that 
the  continuous  increase  of  the  resistance  to  the  mass  transfer 
(due  to  the  growth  of  the  dried  layer)  is  less  significant  than 
the  slow  down  of  the  drying  rate  (as  the  power  input  decreases). 

It  indicates  that  the  process  is  essentially  heat  transfer  con- 
trolled. As  the  ice-front  temperature  decreases  the  entire 
frozen  core  cools  down. 

The  maxima  of  the  total  microwave  power  absorbed  (D.^)  and, 

in  the  frozen  region  only  (O.^,),  (Fig.  7.3)  reflects  the  fact 

r 

that  the  dissipation  coefficient  increases  with  the  temperature 
in  the  frozen  and  dried  regions.  iT  ^ and  Ap  in  Fig.  7.3  have 
been  calculated  by, 

.Qp(t)  = 

and  similarly  for  /Ip,-  The  w's  are  defined  by  Eq.  (3.5). 

The  strong  similarity  of  the  drying  rate  curv  ith  that 


^X(t) 


<0p(x,t)dx 


(7.3) 
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of  the  power  Input,  suggests  again  that  the  heat  transfer  is  the 
controlling  mechanism.  However,  the  resistance  to  the  mass  transfer 
and  the  driving  force  are  significant  as  they  cause  the  actual 
sample  temperature  variations. 

The  re!.ative  positions  of  the  power  input  and  the  drying 
rate  curve  (Fig.  7.3)  require  additional  comments.  As  the  micro- 
wave  power  has  just  been  turned  on,  the  power  input  is  slightly 
less  than  the  actual  power  consumed  by  the  sublimation.  This 
is  due  to  the  extra  energy  coming  from  the  heat  supplied  from 
the  surrounding  atmosphere  to  the  sample,  as  its  surface  temperature 
is  well  below  the  chamber  temperature.  As  more  microwave  power 
is  absorbed  in  the  dried  layer  the  external  surface  temperature 
continuously  increases  (Fig.  7.2)  and  exceeds  the  chamber  temper- 
ature. Thus,  the  heat  flux  at  the  surface,  which  is  initially 
from  the  surrounding  to  the  sample,  continuously  decreases,  and 
finally  reverses  itself  as  Tg  becomes  greater  than  T^^.  This 
effect  is  shown  on  Fig,  7.1  by  the  shape  of  the  temperature 
profiles  in  the  dried  layer  and  by  the  time  history  of  Tg  (Fig.  7.2). 
It  accounts  for  the  fact  that  crosses  over  in  Fig.  7.3. 

2 . Process  variables 

For  a given  sample  composition,  the  microwave  freeze  drying 

process  can  be  characterized  by  the  following  process  variables: 

the  vacuum  chamber  total  pressure  Pj^,  the  vacuum  chamber  partial 

w 

pressure  of  water  vapor  P^,  the  ambient  •♦■emperature  1^^,  the 
initial  temperature  T^,  the  sample  size  or  thickness  for  the 
one -dimensional  case,  and  the  electric  field  peak  strength  in 
the  vacuum  E.  The  operating  conditions  are  determined  by  a given 
set  (E,  P Pp,  T T , L)  of  the  process  variables  which  will 
be  referred  to  as  operation  point.  It  should  be  noted  that  the 
solution  of  the  mathematical  model  does  not  require  that  the 
process  variables  be  constants  anjd  they  could  be  programmed  as 
a function  of  time.  However,  this  possibility  is  not  considered 
in  this  application. 
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Also,  the  charobtr  tmparaturs  Is  raally  mantionad  hara  for 
rafaranca,  but  is  takan  constant  and  aqual  to  20*C  as  it  is  not 
axpactad  to  be  a significant  parameter.  Moreover,  its  measurement 
is  questionable  under  vacuum  since  heat  conduction  is  limited 
and  mostly  radiant  heat  transfer  predominates.  Indeed,  one  could 
think  of  the  heat  transfer  coefficient  h as  a lumped  parameter 
which  contains  all  possible  forms  of  heat  exchanges  between  the 
sample  and  the  surrounding  vacuum  at  (20®C  here). 

Similarly,  the  initial  temperature  of  the  frozen  sample 
is  a difficult  {parameter  to  control  due  to  the  various  manipulations 
which  precede  the  actual  drying.  Fortunately,  it  is  not  expected 
to  be  a significant  parameter  either,  as  dielectric  heating  is 
a rapid  process.  In  effect,  a very  small  amount  of  time  is  necessary 
to  bring  the  material  up  to  the  same  state  as  the  one  which  would 
have  been  obtained  with  a different  initial  temperature.  Cal- 
culations were  made  for  a one-inch  thick  slab  at  two  differe.. . 
initial  temperatures:  -1S*C  and  -30“C  with  an  electric  field 
E - 125  V/cm  and  pressures  * Pr  ~ .075  mmHg.  The  drying 
curves  were  found  to  be  practically  identical.  The  difference 
on  the  normalized  interface  position  5 is  less  than  .002  for 
most  of  the  drying.  The  temperature  profiles  in  both  the  frozen 
and  dried  regions  a.re  essentially  the  same  except  diuring  the 
start-up  stage  of  the  process.  But  the  temperature  difference 
is  only  significant  in  a short  period  of  time  at  the  beginning 
of  the  process  and  reduces  to  less  than  1®C.  As  the  microwave 
power  is  turned  on,  the  original  temperature  difference  becomes 
unnoticeable.  Thus,  a constant  initial  temperature  of  -15°C  is 
assumed. 

3 . Feasible  domain  of  operation 

The  operating  range  for  the  process  variables  is  limited 
to  a finite  domain  which  will  be  referred  to  as  feasible  domain 
of  operation.  This  is  due  to  the  fact  that  exist  upper 

limits  on  the  electric  field  strength  in  the  vacuum  and  the 
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temperature  maxima  in  the  frozen  and  dried  portion  of  the  material 
being  dried.  These  physical  constraints  which  have  been  discussed 
earlier  can  be  expressed  mathematically  by, 


E L Eg(Pj^)  (7.4) 


^F.max.(E-^R*  ^R’  ^ Wting  (^-5) 

Vmax.(^*PR>  Pr-  ^ ^herm.  degrad. 

Constraints  (7.4),  (7.5),  and  (7.6)  are  to  be  satisfied  at  any 

time  in  order  to  ensure  a satisfactory  operation  without  corona 
discharge,  melting  of  the  frozen  material  and  thermal  degradation 
of  the  dried  product. 

It  has  been  demonstrated  (Goulii  and  Kenyon  (t971))  that  the 

break-down  electric  field  strength  in  an  industrial  microwave 

applicator  for  given  gases  is  a function  of  the  vacuum  pressure 

and  the  design  of  the  microwave  applicator.  It  is  particularly 

sensitive  to  the  geometry  of  the  applicator.  Thus,  although 

Constraint  (7.4)  has  a great  practical  c-'oortance,  it  cannot  be 

incorporated  in  this  theoretical  application  as  it  must  be 

determined  for  each  particular  design. 

The  feasible  operation  domain  is  then  determined  by  Constraints 

(7.5),  and  (7.6).  In  the  four -dimensional  space  of  the  process 

w 

variables  (E,  P^^,  L)  this  domain  corresponds  to  the  innermost 


volume  enclosed  by  the  coordinate  frame  planes  and 
surfaces  defined  by  the  tight  constraints. 

the  two  limiting 

(E. 

F,max.  ' ’ 

P,,  P;^,  L)  = -1.5»C 

(7.7) 

T (F 

D,max.  ' ’ 

P^,  P]J,  L)  = 60”C 

(7.8) 

This  is  illustrated  in  Fig.  7.4.  The  limiting 
defined  by  the  tight  constraint 

surface 

"F,max.(^- 

Pj^,  Pp,  1.27)  = -1.5”C 

(7.9) 
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has  baen  calculatad  for  a one-inch  thick  symmetric  slab  (L  • or 
1.27  cm).  The  maximum  electric  field  strength,  allowable  before 
melting  of  the  frozen  material  occurs,  was  determined  at  various 
pressures.  The  accuracy  of  the  electric  field  is  5 V/cm.  The 
feasible  domain  with  respect  to  Constraint  (7.5)  is  the  inner 
volume  comprised  between  the  limiting  surface,  the  first  bisector 
plane  Pj^  * Pp  (since  Pp  ^ Pp)  sind  the  planes  E = 0 and  P^  = 0. 
Thus,  if  the  operating  conditions  correspond  to  a point  inside 
this  volume,  no  melting  will  occur  at  any  time  during  the  process. 
On  the  contrary,  if  an  operation  point  is  taken  outside  this 
domain  it  will  result  in  melting  of  the  frozen  material  in  the 
course  of  drying.  However,  Constraint  (7.5)  is  not  necessarily 
satisfied  and  the  upper  limit  of  60*C  may  well  be  exceeded. 

Figs.  7.5  and  7.6  show  the  orthogonal  projections  of  the 
limiting  surface  in  the  plane  Pp  * 0,  and  P^  ~ 0.  Fig.  7.6 
shows  the  variations  of  the  maximum  allowable  electric  field 
strength  , (such  that  no  melting  occurs)  with  the  water 

vapor  partial  pressure  at  various  total  pressures.  Fig.  7.5 
shows  the  variations  of  E with  the  total  pressure  at  various 
partial  pressures.  Fig.  7.5  indicates  that  melting  occurs  at 
a lower  E value  as  the  pressure  increases.  This  is  to  be  expected 
as  the  diffuslvity  in  the  dried  layer  decreases  with  increasing 
total  pressure,  resulting  in  a higher  mass  transfer  resistance 
of  the  dried  layer.  This  gives  a higher  ice-front  pressure  and 
temperature  for  a similar  drying  rate  (same  E) . As  the  pressure 
approaches  100  mmHg,  the  electric  field  has  to  be  decreased  drasti- 
cally in  order  to  ensure  drying  without  melting.  Similarly, 

Fig.  7.6  shows  that  at  a given  pressure,  an  increase  of  the 
partial  pressur«  ' water  vapor  in  the  vacuum  requires  that  the 
electric  field  be  decreased.  Thus,  the  mass  transfer  driving 
force  is  decreased,  which  results  in  a higher  pressure  and  temper- 
ature at  the  ice-front  in  order  to  maintain  the  same  drying 
rate.  However,  the  plateau  in  the  curves  indicates  that  at 
higher  pressures  Pp  this  effect  becomes  insignificant 
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jntll  the  partial  orassure  Pp  reaches  a high  enough  value,  then 

the  electric  field  must  be  decreased  drastically  again  in  order 

to  keep  the  sample  frosen.  The  triple  point,  which  corresponds 

w 

approximately  to  Pp  = 4.6  mmHg,  constitutes  an  absolute  limit 

boyond  which  it  is  not  possible  to  keep  the  sampl'.-  from  meltiiig. 

w 

Another  absolute  limit  is  constituted  by  the  curve  Pp  = Pp  on 

both  Figs.  7.5  an"  7.6.  It  reflects  the  physical  constraint 
w 

Pp  ^ Pp  since  a partial  pressure  cannot  be  greater  than  the  total 
pressure . 

Considering  the  fact  that  the  drying  time  is  expected  to 

decrease  substantially  as  the  electric  field  is  increased,  Fig. 

7.4  shows  that  an  optimal  operation  would  be  to  use  pressures 
w 

Pp  and  Pp  as  low  as  possible,  as  a higher  field  strength  could 

be  applied  without  melting.  Figs.  7.5  and  7.6  also  show  that  as 
w 

Pp  or  Pp  decrease  and  approach  .2  mmHg,  little  gain  is  made  o“. 
the  maximum  allowable  field  strength.  This  is  due  to  the  fad 
that  at  low  pressures  the  vapor  flow  regime  in  the  dried  layer 
is  mostly  Knudsen  diffusion  and  thus  the  diffusivity  becomes 
independent  of  the  pressure.  This  is  why  a lower  lirit  of 
about  50  ^Hg  is  considered  for  vacuum  freeze-drying.  Lower  pressures 
would  represent  cost  increase  for  the  vacuum  system  but  no  signi- 
ficant gain  in  keeping  the  sample  frozen.  A likfily  point  of 
operation  is  at  pressure,*  Pp  / .2  mmHg  and  Pp  close  to  75  ^Hg, 
using  a low  temperature  condenser  (vapor  trap;  below  -40“C). 

The  entire  feasible  operation  domain  would  be  determined  by 
a similar  study  conducted  at  various  thicknesses,  including  the 
last  Constraints  (7.6)  in  order  to  ensure  also  that  no  thermal 
degradation  of  the  dried  product  will  occur  due  to  excessive 
heating  of  the  dried  layer. 

Figs.  7.7  and  7.8,  show  typical  variations  of  the  maximum 

temperatures  at  any  time  in  the  frozen  and  dried  layers  with  the 

electric  field  strength  at  various  pressures,  for  a vapor  pressure 
w 

Pp  = .075  mmHg.  The  temperature  maximum  at  any  time  in  the  frozen 
region  is  found  to  vary  linearly  with  the  square  of  the  electric 
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- Variation-'  of  !ho  paxirur  1 o?:,  porat  uro  in  the  frozen 
reqion  fat  T^y  tiro)  wit^"  the  electric  field,  at 
variou'^  total  pressure-.  I = 1.2“  cr  ; i’^''  = . )'5  , 


so: 


field  strength  in  the  vacuum.  The  operation  domain  Is  limited  to 

the  region  enclosed  within  the  lines  (1),  (2),  and  (3)  on  both 

figures.  Line  (1)  reflects  the  no  melting  constraint  (T_  / 

r , ma  X , ® 

-1.5*C)  while  line  (2)  corresponds  to  the  phyf.ical  constraint 
w 

(Pp  ^ I-in®  (3)  corresponds  to  the  constraint  of  no  thermal 

degrr.dation  (T_  / 60“C).  It  5 s als:o  found  that  while  the 

temperature  maximum  at  any  time  in  the  frozen  region  is  relatively 
sensitive  to  variations  of  the  total  pressure,  particularly  at 
higher  pressures,  the  temperature  maximum  in  the  dried  layer  does 
not  vary  significantly  with  it.  It  has  been  pointed  out  earlier 
that  the  temperatures  in  the  froaen  region  are  expected  to 
increase  as  the  resistance  to  the  mass  transfer  through  the  dried 
layer  is  increased  when  the  total  pressure  becomes  higher.  The 
decrease  of  the  temperature  maximum  in  the  dried  layer  as  the 
total  pressure  is  increased  can  be  attributed  to  the  increase  of 
the  effective  thermal  conductivity  (Pig.  6.4).  Thus,  for  the 
same  heat  transfer,  the  temperature  in  the  dried  layer  is  expected 
to  be  somewhat  lower.  Both  temperature  maxima  are  found  to  be 
ctreug  functions  of  the  electric  field.  It  can  be  also  seen  from 
Fig.  7.7  and  7.8  that  for  pressures  of  .075  and  .29  mmHg,  the 
process  will  be  actually  limited  by  the  overheating  of  the  dried 
layer  as  the  temperature  exceeds  60®C  before  melting  of  the  frozen 
core . 

Figs.  7.9  and  7.10  show  typical  variations  of  the  maximum 
temperature  reached  in  the  frozen  region  (Tp  ) and  in  the 
dried  layer  (T^^  ),  with  respect  to  the  sample  size  (slab 

half  thickness  L).  The  calculations  have  been  carried  out  for 
Pp  = Pp  = .075  mmHg,  a reduction  of  a 1/4"  thick  slab  by  1/64" 
would  result  approximately  in  a 40®C  drop  of  the  temperature  maximum. 
Although  one  may  not  have  expected  such  a strong  ...  jonder.cy  of 
the  temperature  maxima  upon  the  slab  thickness,  it  is  comprehensible 
for  the  temperature  and  pressure  of  the  ice-front  to  decrease 
as  the  slab  thickness  is  reduced  because  of  the  lower  resistance 
to  the  mass  transfer  in  the  dried  layer,  which  results  in  a general 
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Tempo  rature 


cooling  of  th«  sample.  The  vtliarpnats  of  the  dependency  at  high 
electric  field  strengt' ■ may  be  attributed  to  the  fact  that  the 
energy  generated  in  the  sample  for  a given  thickness  varies  as 
Thus,  the  quantitative  variations  oi  ‘he  temperature  maxima  at 
various  field  strengths  will  constitute  a very  critical  information 
for  the  designer  when  deciding  on  an  optimal  slnb  thickness. 

L.iwer  temperatures  will  . 'w  one  to  increase-  the  microwave 
power  input  and  thus  shorten  the  drying  time . This  point  will 
be  discussed  in  the  next  section. 

Figs.  7. <3  and  7.10  show  that  the  60*C  limit  is  exceeded  at 
the  assumed  operating  pressure?  before  melting  occurs  for  all 
values  of  L considerea.  However,  it  must  be  noted  that  these 
findings  depend  on  the  assumptions  which  we  made  for  the  numerical 
value  of  certain  properties  (particularly  the  diff usivity) . 

It  was  found  (Part  III)  that  the  effective  diff usivity  might  have 
been  overestimated.  The  true  temperatures  in  the  frozen  region 
may,  thus,  be  substantially  higher  than  those  indicated  here. 

4.  Effects  of  the  process  variables  on  the  drvino  time 

Calculations  have  been  carried  out  at  various  total  pressures, 
electric  field  strengths  and  sample  size  in  order  to  study  the 
variations  of  the  total  drying  time  (including  the  start-up 
period)  with  these  parameters.  The  partial  pressure  of  the  water 
vapor  in  the  vacuum  is  a less  important  parameter.  In  effect, 
the  drying  time  is  not  expected  to  vary  significantly  over  the 
small  feasible  range  of  the  partial  pressure  of  water  vapor. 

As  one  may  have  expected,  the  drying  time  is  a strong  function 

of  the  electric  field  strength.  The  drying  time  has  been  plotted 

in  Fig.  7.11  vs  the  square  of  the  electric  field  (the  abscissa 

of  Fig.  7.11  is  a quadratic  scale)  at  various  pressures  and  a 

vapor  partial  pressure  of  .075  mmHg.  From  the  small  curvature 

of  the  curves  it  can  be  seen  that  the  drying  time  varies  approx- 
2 

imately  with  E , within  a relatively  small  range.  In  othpr  words, 
the  drying  time  is  approximately  proportional  to  the  microwave 
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power  input. 

It  is  found  also  that  an  increase  of  the  total  pressure  has 
an  appreciable  effect,  although  small,  in  decreasing  the  drying 
time,  as  shown  on  Fig.  7.12.  This  is  due  to  the  increase  of  the 
dissipation  coefficient  (mainly  in  the  frozen  region)  as  the 
temperature  increases  in  the  sampl«?  with  the  resistance  to  the 
mass  transfer.  This  effect  has  been  discussed  earlier.  On  the 
c her  hand,  the  increase  of  the  ;esi stance  to  the  mass  transfer 
should  result  in  a slow  down  of  the  mass  transfer  rate.  The  fact 
that  the  drying  time  actually  decreases  rather  than  increases  is 
another  indication  that  the  freeze-drying  process  using  dielectric 
heating  is  heat  transfer  controlled.  However,  if  the  pressure 
is  further  increased  the  upper  curve  on  Fig.  7.12  would  be  ex- 
pected to  go  to  a minimum  and  start  to  rise  as  the  process  becomes 
mass  transfer  controlled.  However,  for  this  to  happen  a field 
strength  low  enough  would  have  to  be  used,  otherwise  melting  of 
the  frozen  core  would  occur  first.  The  lower  curve  on  Fig.  7.12 
illustrates  the  effect  of  the  pressure  increase  in  raising  the 
frozen  core  temperature  due  to  the  higher  resistance  to  the  mass 
tr..nsfer  in  the  dried  layer. 

The  different  calculation  points  used  to  determine  the 
limiting  surface  in  Sec.  3 of  the  melting-free  domain  of  operation, 
have  also  been  plotted  in  Fig.  7.11.  It  is  interesting  to  note 
that,  although  they  were  obtained  at  different  pressures  and  vapor 
partial  pressures  of  water,  they  all  fall  on  the  same  smooth  rurve. 
This  suggests  that  the  drying  time  depends  only  on  the  electric 
field  strength.  This  supports  the  argument  presented  in  the 
previous  peiragraph.  The  pressures  affect  the  drying  time  beef  use 
of  their  effect  on  the  temperatures,  which  in  turn  affect  the 
dissipation  coefficient.  Thus,  if  the  temperatures  are  approx- 
imately the  same,  as  it  is  the  case  along  the  limiting  surface, 
the  effects  of  the  pressures  on  the  drying  time  are  nil. 

The  variations  of  the  total  drying  time  with  L in  Fig.  7.13 
seem  to  be  in  contradiction  with  Hoover  et  al.  (1966)  results. 
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Drying  Time,  Hrs. 


rature 


who  have  obtained  experimental  drying  times  essentially  independent 
of  the  thickness.  The  results  shown  in  Fig.  7.13  may  be  supported 
by  an  approximate  calculation  of  the  total  drying  time  t2  as 
shown  below.  The  total  drying  time  can  be  obtained  by  rearranging 
Equation  (4.  7)  to  give 

dt  = -p  0 ^ (7.10) 

i 

The  drying  rate  can  be  related  to  the  microwave  power  input, 
by  neglecting  the  variations  of  the  sensible  heat  and  the  losses 
at  the  outer  surface.  A combination  of  Equations  (7.1)  and  (3.5) 
gives 

“i  =1^  ["f  ^ * *'-D  I*-  - X)] 


where  the  variations  of  Kp  and  with  loc.ation  have  been  neglected. 
Thus, 


dt 


P ^ ^^s  dX 

p2  KpX  + Kj^(L-X) 


(7.12) 


Integration  of  Equation  (7.12)  between  t = 0 and  t = ox 

X = L and  X = 0,  gives  (neglecting  the  variations  of  K and  K 

r U 

with  time) , 


p cr  AH 

^2  " h ^ ^2 

(Kp-  Kj,) 

Equation  (7.13)  shows  that 

reduced,  the  overall  drying  time 

product,  remaining  to  dry  at 

decreases.^  If  the  thickness  is 
Xi 

eventually  t as  „ goes  to  zero 
i L 

the  end  of  the  start-up  period), 
only  radiant  heating,  for  which 
with  the  samp’’ e thickness  L. 


1 + 


L y 


(7.13) 


when  the  sample  thickness  is 
t^  decreases  as  the  fraction  of 
the  end  of  the  start-up  period 
further  reduced,  t^  reaches 
(all  frozen  material  removed  at 
The  drying  process  then  becomes 
the  drying  time  goes  to  zero 
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Half  Thickness,  In. 

Efforts  of  tho  sample  thickness  on  the  total 
dryino  time  for  various  field  strengths. 


The  above  variations  of  with  L are  in  agreement  and 
confirm  those  predicted  by  the  model  as  shown  in  Fig.  7.13. 

It  should  be  added  that  the  dissipation  coefficients  and 
Kp  are  expected  to  affect  the  drying  time  as  the  thickness  is 
changed.  In  effect,  they  vary  with  the  temperature  which  is 
a function  of  the  thickness  L (see  Figs.  7.9  and  7.10).  Similarly, 
although  the  mass  transfer  is  not  the  controlling  mechanism, 
a variation  in  L results  in  a different  resistance  to  the  mass 
transfer  ,Mch  in  turn  may,  directly  or  indirectly,  vary  the 
drying  time. 

Thfi  maximum  of  the  curves  at  higher  thicknesses  in  Fig.  7,13 
can  be  attributed  to  the  temperature  inc,'  "jase  (mainly  in  the  frozen 
core)  due  to  the  increase  of  tr.e  mass  transfer  resistance  when 
the  thickness  is  augmented,  as  discussed  previously.  The  higher 
temperatures  in  the  sample  increase  the  dissipation  coefficient 
in  both  the  frozen  and  dried  regions,  which  causes  more  microwave 
power  to  be  absorbed  by  the  sample,  and  thus,  shortens  the  drying 
time.  Fig.  7.13  seems  also  to  indicate  that  smaller  slab  thick- 
nesses would  enable  one  to  use  higher  electric  field  strengths 
to  shorten  drastically  the  drying  time  '..ithout  melting. 

Fig.  7.14  shows  a plot  of  the  amount  of  frozen  product 
dried  per  hour  (per  unit  area)  vs  the  sample  half  thickness  L 
at  various  electric  field  strengths.  It  can  be  seen  from  this 
figure  that  as  the  slab  thickness  is  reduced  the  output  of  the 
microwave  freeze-dryer  continuously  decreases.  Thus,  for  a given 
electric  field  strength  the  optimal  thickness  yielding  the 
highest  output  of  dried  product  is  the  largest  thickness  without 
occurrence  of  melting  of  the  frozen  material  or  overheating  of 
the  dried  product.  Fig.  7.14  shows  that,  under  the  assumptions 
made  here,  that  overheating  of  the  dried  product  would  occur 
earlier  than  melting  of  the  frozen  core.  However,  it  should  be 
mentioned  again,  and  as  it  will  be  discussed  in  the  experimental 
part  of  this  study,  that  the  numerical  value  of  the  diffusivity 
may  have  been  overestimated  in  this  case.  This  would  m^an  that 
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the  actual  temperatures  in  the  sample  are  indeed  somewhat  higher 
than  those  computed  here  and  thus  melting  may  occur  earlier  than 
predicted. 

On  the  other  hand,  the  negative  slope  of  the  "overheating" 
limiting  curve  in  Fig.  7.14,  indicates  that  us  ' of  smaller  thick- 
nesses would  increase  the  freeze-dryer  output  as  hight  • field 
strengths  are  made  possible.  However,  as  the  electric  field 
strength  is  increased,  a corona  discharge  may  develop  in  the 
microwave  applicator,  reducing  drastically  the  power  conversion 
efficiency  and  possibly  inducing  burns  of  the  product.  Thus, 
the  optimal  thickness  yielding  the  maximum  freeze-dryer  output 
is  that  which  would  correspond  to  an  operation  near  corona  and 
overheating  occurrence  (or  melting).  For  instance,  if  corona 
was  to  develop  slightly  above  255  V/cm  the  optimal  1 hickness 
would  correspond  to  the  optimal  operation  point  shown  in  Fig,  7,14. 

It  is  important  to  note  that  the  previous  discussion  had  assumed 

w 

an  operation  under  low  pressures  (Pj^  = Pj^  = ,075  mmHg)  in  order 
to  ensure  the  lowest  possible  temperatures  in  the  frozen  region 
and  the  dried  layer. 

Under  these  conditicns  and  the  assumptions  made,  if  corona 
were  to  occur  beyond  2o5  7/cm,  drying  times  as  short  as  an  hour 
and  one  half  (including  the  30  min  start-up  period),  should  be 
possible  for  a symmetric  slab.  This  would  result  in  a maximum 
freeze-drver  output  of  approximately  0.4  grams  of  frozen  product 
per  hour,  per  unit  area. 
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- PART  III 


CORRELATION  WITH  EXPERIMENTAL  DRYING  CURVES 


Chapter  VIII;  Experimental  Set-Up  and  Procedure 

Chapter  IX:  Experimental  Results  - Correlation  with  Simulation 
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- Chapter  VIII  - 


EXPERIMENTAL  SET-UP  AND  PROCEDURE 


1 . Experimental  set -up 

1.1  The  microwave  freeze-drver 

A schematic  of  the  experimental  system  is  given  in  Fig.  6.1. 
The  microwave  generator  (A)  is  a converted  microwave  oven  (Speed 
Range  110)  which  uses  a magnetron  Amperex  DX-206  to  deliver  1 KW 
of  adjustable  power  at  2450  MHz  ^ 25  MHz.  The  microwave  applicator 
cavity  (B)  is  located  inside  a vacuum  chamber  (C)  and  contains 
the  sample  to  be  dried  (D).  The  microwave  power  is  transriitted 
to  the  cavity  (B)  by  a rectangular  wave  guide  WR-264,  through 
various  control,  monitoring  and  tuning  devices.  The  power  from 
the  generator  is  sent  by  a four-port  circulator  (E)  into  the  line 
section  where  an  H-tuner  (F)  varies  the  admittance  to  a dry 
resistor  (G)  (matched  load  of  less  than  1%  reflection).  The 
fraction  of  power  reflected  at  the  H-tuner  (determined  by  the 
plunger  location)  travels  back  to  the  circulator  where  it  is 
diverted  into  the  termination  line  section  (microwave  applicator 
and  freeze-dryer ) . It  encounters  a bidirectional  coupler  (H) 
which  separates  the  power  traveling  forward  from  that  traveling 
reverse  (reflected  from  the  termination).  The  reverse  power 
travels  back  to  the  circulator  where  it  is  sent  to  a second  dry 
resistor  (I)  in  order  to  prevent  harmful  reflections  to  the 
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uro  H.l  - Schematic  nf  the  evperimental  set-up 


magnetron.  The  forward  and  reverse  powers  in  the  bidirectional 
coupler  are  sampled  by  two  small  antenna  probes.  The  bidirectional 
coupler  is  connected  to  a twist  (J)  which  rotates  90*  the  original 
polarization  of  the  electric  field.  The  electric  field  enters 
the  cavity  (B)  with  a vertical  polarization  (parallel  to  the 
figure  plane).  A 3-Stub  tuner  (K)  is  used  in  conjunction  with  an 
E-H  tuner  (L)  to  tune  the  cavity  (B).  The  microwave  power  enters 
the  cavity  (B)  through  a wave  guide  vacuum  transition  (M)  using 
a fiber  glass  window. 

The  multimode  rectangular  cavity  (B)  has  dimensions  39  x 39 
X 51  cm  and  is  made  of  perforated  aluminum  sheet  in  order  to  retain 
the  microwavj  power,  but  let  the  water  vapor  from  the  dehydration 
flow  through  to  the  main  chamber  (C)  and  finally  the  vapor  trap. 

An  oversized  condenser  is  used  for  the  vapor  trap.  The  sample  to 
be  freeze-dried  (D)  is  hung  inside  the  cavity  (B)  with  its  main 
surface  horizontal  and  is  located  in  front  of  the  wave  guide 
mouth  as  shown  in  Fig.  0.1.  It  is  suspended  to  a cantilever  beam, 
used  in  the  weight  measuring  system,  by  means  of  teflon  hooks 
and  fiber  glass  strings. 

Finally,  n-propyl  alcohol  is  circulated  at  a constant  flow 
rate  in  a glass  loop  inside  the  cavity.  It  was  intended  to  serve  a 
double  purpose.  First,  in  order  to  obtain  a fairly  uniform  field 
in  the  multimode  cavity,  a reasonable  load  is  necessary.  By 
increasing  the  load  one  decreases  the  Q of  the  cavity  and  broadens 
its  frequency  band  at  resonance.  Thus,  more  modes  (eigen  frequencies) 
can  be  admitted  inside  the  cavity  which  superimpose  and  average  out 
to  yield  a relatively  uniform  field.  The  increase  of  the  load  was 
found  necessary,  due  to  the  small  size  of  the  sample  used  (about 

3 

60  to  90  cm  ) as  compared  to  the  large  volume  of  the  cavity  (about 
4 3 

7.5  X 10  cm  ) . n-propyl  alcohol  was  chosen  as  the  circulating 
fluid  because  its  dielectric  constant  (€'=  3."')  is  comparable  to 
that  of  frozen  beef  meat  (€'  = .3.6).  Second,  it  was  intended 
as  a means  of  monitoring  the  electric  field  inside  the  cavity, 
as  it  will  be  discussed  below. 
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1.2  Th«  datfc  aLCQuiiition  «v»t«n 

A ep«ci»l  lyittm  was  dsslgnsd  by  this  Invsstlgstor  to  allow 
continuous  rscording  of  ths  sample  weight  in  the  course  of  the 
dehydration.  The  sample  is  suspended  by  means  of  a fiber  glass 
string  to  the  free  end  of  a cantilever  beam  through  a knife  edge. 

The  strain  induced  at  the  fixed  end  by  the  sample  weight  is  measured 
by  four  strain  gages  mounted  in  a Wheatstone  bridge  (two  in 
compression  and  two  in  tension).  The  bridge  is  excited  with  a 
1.5  V DC  voltage  and  balanced  for  no  load  conditions  (no  sample 
hung).  The  DC  signal  resulting  from  the  imbalance  of  the  bridge 
in  the  measuring  circuit,  when  the  sample  is  hung,  is  a measure 
of  the  weight.  This  DC  signal,  of  the  order  of  several  ^V,  is 
amplified  and  sent  to  a strip  chart  recorder.  Known  weights  are 
used  to  calibrate  the  measuring  system.  The  recorder  gain  is 
adjusted  to  give  a reading  of  10  g/in  on  the  chart.  Surrounding 
machinery  induces  undesirable  mechanical  vibrations  in  the  beam 
which  result  in  a strong  low  frequency  noise  in  the  output  of  the 
DC  amplifier.  Noise  of  a frequency  less  than  10  Ha  could  not  be 
suppressed  by  the  built-in  filters  of  the  amplifier.  A special 
filter  was  designed  by  this  investigator  to  suppress  all  AC 
signals  in  the  range  1 to  10  Ha  and  proved  very  successful  in 
filtering  out  the  electrical  noise  caused  by  the  mechanical 
vibrations.  The  assembly  (DC  power  supply,  balancing  unit,  DC 
amplifier  and  strip  chart  recorder)  is  referred  to  in  Fig.  8.1 
as  data  acquisition  unit  1. 

The  water  vapor  partial  pressure  in  the  vacuum  chamber  is 
measured  continuously  by  a Panametrics  model  1000  Hygrometer 
(unit  2)  which  uses  an  aluminum  oxide  transducer  probe.  The 
total  pressure  is  monitored  by  a MacLeod  Pressure  gage  (unit  3). 

As  mentioned  above,  an  attempt  was  made  to  estimate  the 
average  electric  field  in  the  cavity  by  measuring  the  amount  of 
microwave  power  absorbed  in  n-propyl  alcohol  circulated  in  a glass 
loop  through  the  cavity  at  a constant  flow  rate.  The  microwave 
power  loss  in  the  alcohol  is  determined  from  the  temperature 
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diffarcnca  batwaan  tha  Inlat  and  tha  outlat  maaaurad  with  a Doric 
digital  tharnocoupla  tamparatura  indicator  (unit  4).  Howavar, 
only  a qualitative  monitoring  of  tha  avarage  flald  strength  was 
possible  due  to  the  difficulty  of  measuring  the  effective  volume 
of  alcohol  which  absorbed  microwave  power.  Moreover,  the  con- 
tamination of  the  alcohol  by  the  Tygon  tubing  connections  makes  the 
validity  of  utilizing  the  dielectric  data  of  the  alcohol  questionable. 
However,  this  system  is  used  as  an  indicator  of  melting  occurence 
in  the  frozen  sample.  In  effect,  when  melting  occurs,  as  the 
dissipation  factor  of  the  wet  meat  is  much  greater  than  that  of 
the  frozen  meat  the  power  absorption  in  the  sample  rises  suddenly. 
Thus,  at  the  same  time,  the  microwave  power  absorbed  in  the  alcohol 
drops  as  indicated  by  the  outlet  temperature  drop.  This  is  due 
to  the  fact  that  the  microwave  power  fed  into  the  cavity  is  relatively 
constant . 

The  total  microwave  power  loss  it.  the  termination  line 
(including  twist,  tuners,  window,  connectors,  leakage  sample 
and  alcohol)  is  measured  from  the  forward  and  reverse  microwave 
powers  'atected  in  the  bidirectional  coupler.  The  detected  for- 
ward and  reverse  powers  are  determined  by  a Hewlett-Packard  431 
B power  meter  wnich  uses  a thermistor  mount  (unit  5)  A manual 
switch  allows  one  to  select  either  powers.  The  power  measurement 
was  initially  intended  to  monitor  the  power  los^  in  the  sample. 
However,  a quantitative  measurement  was  impracticable  due  to  t lu‘ 
relatively  high  fraction  of  power  lost  in  leakage.  Determination 
of  this  fraction  was  not  possible  as  it  was  found  to  be  a function 
of  the  load  in  the  cavity,  anv.  the  tuning  adjustment.  Thus,  the 
results  of  the  power  me.asuring  unit  were  limited  to  comparisons 
between  runs.  The  Hewlett-Packard  power  meter  reading  (mW) 
was  calibrated  against  a factory  calibrated  power  me*er  (Bird 
Electronic,  Thru-line).  The  latter  power  meter  could  not  be  kept 
for  on-line  direct  power  measurements  as  it  cculd  only  give  re- 
liable readings  under  low  reflection  ratios. 
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2 General  experiment^.-l  procedure 

2 . 1 Tuning  of  cavity 

Tuning  of  the  cavity  is  necessary  in  order  to  obtain  a 
relatively  uniform  electric  field.  The  uniformity  of  the  field 
;vas  checked  on  a test  run  by  visual  examination  of  a partially 
dried  sample,  cut  in  several  directions.  A planar  ice-front  was 
observed  which  indicated  a relatively  uniform  heating.  The  tuning 
consists  in  performing  an  analysis  of  the  admittance  vs  frequency 
of  the  termination  line  for  a calibrated  signal  generated  by  a 
sweep  oscillator  (cold-^test ) . The  frequency  band  swept  by  the 
test  oscillator  corresponds  to  the  frequency  band  of  the  magnetron: 
2450  MHs  ^ 25  MHs.  The  tuning  criterion  is  to  adjust  the  3-Stub 
or  the  E-H  tuner  in  order  to  obtain  a flat  profile  of  the  power 
difsipated  in  the  termination  line  vs  frequency.  A flat  power 
profile  means  that  microwave  power  will  be  admitted  into  the 
termination  line  equally  for  each  frequenciis  of  the  frequency 
band  generated  by  the  magnetron.  It  was  hoped  that  this  would 
also  correspond  to  a maximum  number  of  modes  (a  mode  corresponds 
to  a given  eigen -frequency  or  resonant  frequency  of  the  cavity) 
admitted  in  the  cavity  and  thus  yield  a more  uniform  field.  With 
this  tuning  criterion  it  is  believed  that  a relatively  uniform 
field  distribution  can  be  obtained. 

The  termination  line  section  tested  consisted  of  the  line 
section  located  to  the  right  hand  side  of  the  bidirectional  coupler 
in  Fig.  8.1. 


2.2  Sample  preparation 

Two  batches  of  samples  of  different  choices  of  beef  meat 
with  different  preparation  were  used.  The  first  one,  used  in  the 
study  of  the  effects  of  the  power  input  on  th^  drying  curves 
consisted  of  an  eye  of  the  round  choice  which  offers  a rather 
homogeneous  fat  free  texture  (moisture:  71%).  The  outer  fat  was 
trimmed  and  the  meat  deep  frozen  to  approximately  -syc.  Slabs 
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of  approxlm*t«ly  1.5  cm  thlcknast  war*  cut  p*rp*ndlcul*rly  to  th* 
flb*rs  from  the  frosr  i^uacl*  with  an  electric  band  saw.  The 
samples  were  then  wrapped  separately  and  stored  in  a freezer  until 
later  use.  However,  there  were  two  drawbacks  to  this  method  of 
preparation.  The  thickness  and  the  contour  of  the  slab  were  irvegular 
and  difficult  to  reproduce.  Thus  another  preparation  method  was 
employed  in  the  study  of  the  effects  of  pressure  on  the  drying 
curves.  However,  sirce  a larger  muscle  cross  section  was  needed, 
a different  choice  of  meat  was  used.  It  consisted  of  top  of  the 
round  (average  moisture:  67X,  average  fat:  2%).  The  whole  muscle 
was  tempered  at  about  -4‘’C  and  then  pressed  and  formed  in  a cylindrical 
die,  with  the  general  orientation  of  the  fibers  parallel  to  the 
cylinder  axis.  The  meat  cylinder  was  then  chopped  automatically 
in  slabs  of  a relatively  constant  thickness  (approximately  1.5  cm). 

7.5  cm  diameter  discs  were  then  punched  out  of  the  meat  slab 
with  a metal  cutting  form  and  then  placed  on  a metal  tray  and  re- 
frozen in  a blast  freezer.  It  should  be  noted  that  the  freezing 
method  does  not  seem  adequate  since  the  fast  freezing  from  the 
bottom  of  the  sample  resting  on  the  tray  yields  a concave  top 
surface.  A different  method  might  have  given  better  results. 

The  frozen  samples  were  again  wrapped  separately  in  paper  and 
stored  in  a freezer  at  about  -37 *C. 

2.3  Experimental  procedure  for  a run 

Before  each  run,  a known  weight  was  suspended  to  the  cantilever 
beam  in  place  of  the  sample  and  the  reading  of  the  weight  on  the 
recorder  was  allowed  to  reach  a constant  value.  This  normally 
took  one  to  two  hours.  Once  the  reading  we«  constant,  the  gain 
of  the  recorder  was  ad-*jsted  for  a direct  reading  in  grams.  The 
refrigeration  system  and  '’acuum  pump  were  started  30  min  to  1 hr 
prior  to  each  run  with  the  vacuum  chamber  isolated  in  order  to 
achieve  a constant  condenser  temperature  of  approximately  -40°C. 

The  meat  sample  was  then  taken  from  the  storage  freezer, 
weighed  and  its  cross  section  area  and  thickness  were  measured. 
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Tha  sampla  was  than  attachad  to  tha  flbar  glass  string  from  the 
cantllavar  baan  walghing  systam  by  maans  of  thraa  taflon  hooks, 
its  surfaca  horizontal  and  Its  adga  facing  tha  mouth  of  tha 
mlcrowava  powar  faad  (wava  gulda).  Thus,  tha  polarisation  of  tha 
alactrlc  flald  at  the  sampla  surfaca  was  axpactad  to  ba  mainly 
that  In  the  wave  guide  (TEj^  , l.e.,  vertical  (normal  to  the 
sample  surface). 

The  doors  of  the  microwave  cavity  and  vacuum  chamber  were  then 
shut  and  the  chamber  opened  to  the  vacuum.  The  desired  vacuum 
was  found  to  establish  in  less  than  5 min.  A pressure  less  than 
200  ^Hg  was  not  possible  due  to  air  leakage. 

After  usually  30  trin  the  microwave  power  was  turned  on.  The 
desired  power  Input  was  selected  by  using  the  H-tuner  (F  In  Fig. 
6.1).  The  sample  was  then  allowed  to  dry.  It  was  considered  fully 
dried  when  the  weight  reading  on  the  recorder  reached  a constant 
value.  Readings  were  taken  at  variable  time  Intervals  for  the 
following  parameters:  partial  pre.ssure  of  water  vapor,  total 
pressure,  forward  and  reverse  microwave  power,  temperature  IN  and 
OUT  of  the  alcohol  glass  loop,  and  flow  rate  of  the  circulated 
alcohol  (a  constant  flow  rate  could  not  be  achieved). 

The  drying  curve  (sample  weight  vs  time)  was  directly  given 
by  the  reading  on  the  strip  chart  recorder. 

It  should  be  noted  that  in  certain  instances,  especially 
during  the  runs  at  high  pressure  (3.5  mmHg),  a corona  discharge 
occured  at  the  end  of  the  dehydration.  The  high  pressure  runs 
(3.5  mmHg)  were  not  relevant  due  to  local  melting  and  vaporization 
of  the  frozen  core.  An  otherwise,  corona-free  operation  was 
reached. 

The  microwave  power  measurements  were  not  quantitative,  as 
mentioned  earlier.  Thus,  the  electric  field  peak  strength  in  the 
vacuum  near  the  sample  surface  was  evaluated  from  the  drying  rate 
using  the  dielectric  data  of  Kan  and  Yeaton  (1961).  This  point 
will  be  discussed  in  the  next  chapter,  sec.  2. 
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- Chapter  IX  - 


EXPERIMENTAL  RESULTS  AND  CORRELATION  WITH  MATHEMATICAL  MODEL 


Freeze-drying  experiments  have  been  conducted  tor  beef  samples 
with  microwave  dielectric  heating  at  2450  MHz  in  order  to  compare 
experimental  drying  curves  with  those  predicted  by  the  theory. 

The  experimental  set-up  and  procedure  which  were  used  in  the 
experiments  are  described  in  the  previous  chapter.  The  beef 
samples  had  the  geometry  of  a slab  of  about  7.5  cm  in  diameter  and 
a thickness  of  approximately  1.5  cm.  Two  different  choices  of 
beef  meat  were  used  (eye  of  the  round  and  top  of  the  round). 

The  frozen  samples  were  prepared  by  two  different  methods  (see 
chapter  VIII,  sec.  2.2). 

1 . Experimental  results 

Drying  curves  have  been  obtained  at  various  microwave  power 
inputs  for  a constant  vacuum  pressure  of  approximately  .3  mmHg 
(referred  to  as  power  runs)  and  at  various  pressures  for  a fixed 
microwave  power  input  (referred  to  as  pressure  runs)  corresponding 
to  an  electric  field  peak  strength  of  about  125  V/cm. 

Two  different  batches  of  samples  were  used  for  the  power  and 
pressure  runs.  An  eye  of  the  round  choice  of  meat  was  used  for 
the  power  runs,  while  the  pressure  runs  were  performed  with  top 
of  the  round.  The  two  batches  further  differed  by  the  geometry 
of  the  slabs  and  the  preparation  of  the  samples,  as  discussed 
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in  th«  previous  chapter.  It  should  be  noted  that  the  thickness 
of  the  slabs  was  not  uniform  for  any  sample  and  varied  between 
different  samples  due  to  difficulties  in  their  preparation.  A 
variation  of  1 to  2 mm  about  the  mean  thickness  (approximately 
1.5  cm)  was  observed  for  a sample. 

Table  9.1  shows  the  operating  conditions  used  for  each  run 
and  the  average  value  used  to  calculate  the  theoretical  drying 
curves.  The  values  given  for  the  vacuum  pressure  and  partial 
pressure  of  water  vapor  correspond  to  the  starting  point  when  the 
microwave  power  is  turned  on.  They  were  observed  to  decrease 
slowly  during  most  of  the  drving  (a  drop  of  about  .05  to  .1  mmHg) 
and  relatively  faster  in  the  final  stage  during  which  an  additional 
drop  of  approximately  .05  mmHg  was  observed.  The  initial  temper- 
ature is  not  given  as  it  was  not  possible  to  measure  it  directly 
from  the  sample  when  enclosed  in  the  microwave  oven.  The  results 
from  the  theoretical  study  have  shown  (Chapter  VII)  that  the 
starting  temperature  is  not  a significant  parameter  of  the  process. 

As  discussed  in  the  previous  chapter,  the  net  microwave 
power  fed  into  the  cavity  could  not  be  used  to  determine  the  power 
input  into  the  sample  or  the  corresponding  electric  field  strength 
in  the  vacuum.  This  was  due  to  unknown  power  losses  in  the 
termination  line.  Moreover,  the  field  distribution  in  the  cavity 
and  thus  the  power  loss  in  the  load  were  observed  to  be  a function 
of  the  load  and  the  tuning  adjustment.  Therefore,  the  power 
measurement  was  used  primarily  for  comparison  between  runs  (al- 
though this  was  not  always  possible  as  the  tuning  adjustment  was 
changed)  and  for  detecting  anomalies  during  a run.  It  was 
particularly  helpful  to  detect  occurence  of  melting  (shown  by 
a drop  of  the  power  absorbed  in  the  alcohol  loop  as  more  power 
is  absorbed  in  the  sample)  and  of  a corona  discharge  in  the  cavity 
(shown  by  a jump  of  the  reflected  power  in  the  transmission  line 
as  the  admittance  of  the  cavity  drops). 

Assuming  that  the  variations  of  the  sensible  heat  of  the 
sample  and  the  heat  losses  to  or  from  the  surrounding  atmosphere 
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Table  9.1  - Experimental  conditi 


were  negligible  compared  to  the  microwave  power  generated  in  the 
sample,  the  electric  field  strength  in  the  vacuum  (n>^ar  the  sample 
surface)  was  estimated  from  the  drying  rate.  With  this  assumption 
the  power  generated  by  the  microwave  field  in  the  sample  could 
be  eq'oated  to  that  consumed  by  the  sublimation: 

<5>T  = m (9.1) 


where  the  total  microwave  power  generated  can  be  further 
related  to  the  electric  field  strength  by  Eq.  (3.5)  to  give 


4>t  = 


) vp 


(9.2) 


The  dissipation  factors  and  which  are  a function  of  the 
temperature  were  estimated  from  the  numerical  data  shown  in 
Chapter  VI  by  assuming  average  temperatures  for  the  frozen  region 
and  dried  layer.  An  average  temperature  of  5°C  was  assumed  for 
the  dried  layer.  An  average  temperature  between  -10°C  and 
-1.5*C  was  assumed  for  the  frozen  region,  as  Kp  does  not  vary 
significantly  between  -lO'C  and  0°C  (see  Chapter  VI).  However, 
a temperature  lower  than  -10“C  may  be  expected  in  the  frozen  region 
for  the  run  at  the  lowest  power  level  (100  V/cm  in  Fig.  9.1). 

Thus,  feedback  from  the  theory  was  used  to  estimate  the  temperature 
in  the  frozen  core  (-20'C).  The  volumes  Vp  and  V^  of  the  .^rozen 
and  dried  regions  were  determined  using  the  interface  position, 
derived  from  the  drying  curve  at  the  point  where  the  calculations 
were  made  (usually  half  dehydration  point  and/or  maximum  of  the 
drying  rate  curve).  An  average  16%  shrinkage  (as  observed  exper- 
imentally) was  further  allowed  for  the  determination  of  the 
volume  of  the  dried  layer.  Table  9.1  shews  the  results  of  the 
calculations  for  the  electric  field  peak  strength  in  the  vacuum. 

The  starting  fraction  of  the  initial  moisture  content, 

^1  ~ ^'"l  ~ ) (corresponds  to  the  time  when  the  micro- 

wave  power  is  turned  on),  varied  from  one  run  to  the  other  due 
to  difficulties  in  reproducing  the  operating  conditions.  Thus, 
the  microwave  drying  curves  have  been  first  normalized  to  permit 
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their  comparison.  The  normalized  values  were  then  corrected  to 
the  average  starting  moisture  content  (.94)  to  include  the  start 
up  stage  of  the  process  under  average  operating  conditions. 

Thus,  the  following  quantity: 


m - 


i»  = 


™f 


m^^-  l,aver. 


(9.3) 


has  been  plotted  vs  time  in  Figs.  9.1  and  9.2.  The  average  value 
of  observed  for  boch  the  pressure  and  the  power  runs  is  .94. 
Typical  start-up  p>ortions  of  the  drying  curve  have  been  plotted 
in  Fig.  9.1,  using  data  of  Run  16  and  Fig.  9.2,  ujing  data  of 
Run  27 . 

Fig.  9.1  s.jows  the  drying  curves  obtained  at  various  electric 
field  strengths  for  an  average  chamber  pressure  of  .3  mmHg. 

The  common  portion  during  the  start-up  st  ■*>  of  the  process 
(before  the  microwave  power  is  turned  on)  exhibits  the  same 
characteristics  as  that  predicted  by  the  theory  in  Chapter  VII. 

An  initially  high  drying  rate  exists  as  the  dried  layer  is  very 
small  and  the  mass  transfer  driving  force  high  due  to  the  low 
partial  pressure  of  vapor  in  the  vacuum.  It  decreases  rapidly 
as  the  temperature  of  the  frozen  material  drops  to  the  wet -bulb 
temperature  of  the  vacuum  and  becomes  practically  constant  until 
the  microwave  power  is  turned  on.  The  drying  rate  starts  to 
increase  as  microwave  energy  is  absorbed  in  the  sample. 

It  can  be  seen  from  Fig.  9.1  that  the  drying  rate  increases 
substantially  when  the  electric  field  strength  is  increased  as 
expected  from  the  theory.  Toward  the  end  of  the  drying  the  curves 
start  to  level  off.  This  corresponds  to  the  stage  in  which  the 
frozen  core  has  disappeared  and  the  residual  adsorbed  moisture  is 
being  removed  from  the  dried  product.  This  latter  stage  is  usually 
of  little  interest  as  only  removal  of  the  frozen  moisture  content 
is  sought  in  industrial  freeze-drying  operations.  The  residual 
moisture  content  can  be  determined  from  a plot  of  the  drying  rate 
vs  time  as  the  curve  exhibits  a sudden  decrease  of  the  slope  when 
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Fifiure  9.1  - Experimental  drying  curves  (power  runs). 


the  froMn  core  disappears  (less  microwave  power  is  absorbed  by  the 
sample).  This  is  shown  In  Figs.  9.3,  9.4,  9.5,  and  9.6.  It 
has  been  found  by  this  technique  that  the  residual  moisture 
content  is  less  than  5%  of  the  initial  sample  weight  or  7%  of  the 
initial  weight  of  water.  The  actual  residual  mcisture  content 
was  found  to  be  approximately  1.5*  of  the  ini'.ial  sample  weight 
for  all  runs  except.  Runs  No.  12  and  15  for  which  it  was  of  the 
order  of  4.3*.  This  is  reflected  on  Fig.  9.1  by  the  early  leveling 
off  of  the  corresponding  drying  curves. 

It  should  be  noted  that  a visual  examination  of  the  dried 
sample  obtained  in  Run  No.  13  showed  local  burns  of  the  material 
which  indicated  a non  uniform  field  distribution  for  this  run. 

The  tuning  of  the  cavity  was  readjusted  following  this  run  and  the 
problem  was  corrected. 

When  a higher  microwave  power  is  used,  melting  of  the  frozen 
material  starts  to  occur  (Run  No.  15),  followed  by  a local  vapor- 
ization (Runs  No.  16  and  No.  14).  The  local  vaporization  indicates 
that  the  dielectric  heating  is  not  perfectly  uniform  although  a 
relatively  planar  ice-front  has  been  observed  at  medium  and  low 
power  levels.  This  was  done  by  cutting  in  different  directions 
a partially  dried  sample.  It  seems  that  for  the  lower  power  input, 
the  drying  rate  is  small  enough  to  permit  the  temperature  non- 
uniformities in  a cross  section  to  leve"*  out.  This  is  less 
likely  to  occur  at  the  higher  power  levels.  Thus,  a non  planar 
ice-front  and  hot  spots  may  develop  in  the  frozen  region  inducing 
an  early  melting  of  the  frozen  core. 

Melting  is  believed  to  start  occuring  with  Run  No.  15  as  the 
drying  rate  exhibits  (Fig.  9.3)  a sudden  increase  at  time  t = .8  hr. 
This  was  confirmed  during  the  experiments  by  the  observed  drop 
of  the  microwave  power  loss  in  the  alcohol  loop.  Indeed,  as  a 
higher  microwave  power  is  used  (Run  No.  16  and  No.  14)  local 
vaporization  takes  place.  As  the  ice  melts  and  vaporization 
occurs  a much  higher  drying  rate  develops  locally  due  to  the 
higher  dissipation  factor  of  the  wet  material.  This  phenomenon 
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can  b«  obsorvad  visually  as  a dark  spot  and  puffing  develop  at 
the  surface  of  the  sample.  The  occurence  of  vaporization  is 
also  reflected  by  the  drying  curve  as  itc-  slope  increases  due  to 
the  higher  drying  rate  (t.ee  Fig.  9.1). 

Fig.  9.2  shows  similar  plots  for  the  drying  curves  at  various 

pressures  and  a fixed  miciowave  power  input.  An  electric  field 

peak  strength  of  125  V/cm  has  been  estimated  for  the  three  runs 
by  the  method  outlined  earlier.  It  can  be  seen  from  this  figure 
that  the  sample  dried  faster  fo.-  a pressure  of  1 mmHg  than  for 
.2  mmHg  as  expectec  from  the  theory.  This  can  be  explained  by 
the  increase  of  the  dissipation  coefficient  with  the  temperature, 
mainly  in  the  frozen  core,  due  to  the  higher  resistance  of  the 

mass  transfer  in  the  dried  layer.  In  fact,  when  the  pressure  of  the 

vacuum  chamber  is  increase^  to  3.5  mmHg,  melting  of  the  frozen  core 
occurs,  early  during  the  dehydration,  and  is  followed  by  a local 
vaporization . 


2.  Comparison  with  simulated  drying  curves 

In  order  to  obtain  theoretical  drying  curves  for  comparison 
with  those  obtained  experimentally,  the  hef.t  transfer  coefficient 
has  been  estimated  from  the  start-up  portion  of  the  drying  curves. 
This  was  done  by  assuming  a linear  temperature  profile  in  the  dried 
layer.  The  Dew-Point  temperature  corresponding  to  the  pressure 
in  the  vacuum  chamber  was  us’d  to  approximate  the  ir*»-.xont 
temperature  (essentially  equal  to  ’•he  wet -bulb  temp>erature  of  the 
vacuum)  as  the  drying  rate  was  small.  It  can  be  seen  from  an 
energy  balance  on  the  dried  layer  that. 


(Tc 


h(T, 


^s,i)  = ^ 


-Vi 


= Qs  1 = .H^W^ 


(9.4) 


h can  then  be  related  to  the  drying  rate  before  the  microv;ave 
power  is  turned  on  by 
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/iH 


^ ^ W 

^R-  ^i,l  ^ 


(9.  5) 
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wh«r«  p it  givtn  by 
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■ Biot  Numbtr  ■ 


(9.6) 


and  is  the  dried  layer  thickness. 

Values  of  h calculated  in  this  manner  are  shown  in  Table  9.1 

-4  2 

for  each  run.  The  average  value  of  1.9  x 10  cal/cm  /sec/*C 
found  for  the  first  batch  (eye  of  the  round)  was  used  in  the 
simulation.  Due  to  the  concave  shape  of  the  slab  in  the  second 
batch,  the  thickness  of  the  dried  layer,  could  not  be  determined 
accurately.  Thus,  this  might  result  in  a higher  uncertainty  on 
h. 

It  ha-^  been  found  from  preliminary  computer  runs  using  the 
average  operating  conditions  of  the  power  runs,  that  melting 
would  occur  for  an  electric  field  peak  strength  of  205  V/cm. 
However,  an  experimental  value  of  approximately  125  V/cm  was 
observed.  Although  this  discrepancy  could  be  attributed  experiment 
tally  to  a non-uniform  heating,  it  was  felt  that  the  diffusivity 
values  used  in  the  calculations  might  be  overestimated.  The  actual 
values  for  the  beef  samples  could  be  somewhat  lower  by  as  much 
as  a factor  2.  Consequently,  Eq.  (6.8)  in  Chapter  VI  was  modified 
using  Harper's  value  of  D*  for  beef  (corrected  at  265*K)  and 
data  on  (at  265*K)  for  turkey  breast  (no  experimental  value  of 
D|^  could  be  found  for  beef  in  the  literature).  It  constitutes 
a reduction  of  the  diffusivity  by  as  much  as  40%  at  high  pressures 
(several  mmHg  and  over).  The  difference  becomes  negligible  at 
low  pressures  (near  .1  mmHg).  This  modification  led  to  a new 
value  of  the  maximum  field  (without  melting)  of  170  V/cm  (peak 
strength ) . 

Theoretical  drying  curves  were  calculated  using  the  average 
opeiation  conditions  (shown  in  Table  9.1). 

In  order  to  compare  the  theoretical  curves  with  experimental 
values,  the  experimental  drying  curves  were  modified  in  order  to 
represent  the  variations  of  the  ice-front  position  with  time. 
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In  tfftct,  tht  let-front  position  can  bt  rolatod  to  tho  wolght 
loss  by, 


S(t) 


m - (m.+  tn  ) 
f r ' 


m(t)  - R>j.(t) 


in  “ in . in  “ in 
of  o ; 


(9.7) 


is  the  weight  of  residual  water  adsorbed  in  the  dried  layer, 

and  m / (m  - m.)  the  fraction  of  adsorbed  moisture, 
r o i. 

Since  m^(t)  was  not  known  it  was  estimated  from  its  final 

value  m . assuming  a linear  variation  with  time,  m _ was 
r,2  ^ r,2 

obtained  from  the  drying  rate  curves  (Figs.  9.3,  9.4,  9.5,  and 
9.6).  The  following  formula  was  used  to  determine  the  position 
uf  the  interface  at  any  time  from  the  experimental  drying  curves, 


S(t) 


m ( t ) - m^ 


"o‘“f 


(9.8) 


Figs.  9.7,  9.8,  9.9,  9.10,  and  9.11  show  olots  of  the  variations 
of  the  interface  position  determined  from  the  experimental  results 
by  the  above  method.  The  theoretical  curves  predicted  by  the 
simulation  are  also  shown  in  these  figures.  The  electric  field 
strengths  used  in  the  simulation  were  adjusted  in  order  to  give 
the  best  fit.  It  should  be  pointed  out  that  the  final  values  of 
the  electric  field  strength  used  to  calculate  the  curves  shown 
in  Figs.  9.7,  9.8,  9.9,  9.10,  and  9.11  agree  within  4%  of  th' se 
determined  experimentally.  A generally  good  agreement  is  also 
found  between  the  theoretical  and  experimental  drying  curves  as 
shown  by  Figs.  9.7,  9.8,  9.9,  9.10,  and  9.11.  The  small  differences 
observed  may  just  be  the  inherent  differences  between  the  model 
and  the  experiments.  However,  they  may  also  reflect  an  uncer- 
tainty in  the  numerical  data  used  in  the  model  (especially  diffus- 
ivity  and  dissipation  factor)  or  non-ideal  experiment.al  conditions, 
e.g. , 

1.  non-uniform  field  distribution  in  the  sample 

2.  non-constant  standing  wave  pattern  in  the  cavity 

3.  non -constant  microwave  power  fed  into  the  cavity 

4.  non-uniform  thickness  of  the  sample 
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Figure  9.7  - Experimental  (Run  No.  12)  and  simulated  drying  curves.  P = .3  mnHg 
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Figure  . H - Experimental  (Fun  No.  1 ) and  simulated  dryino  curves.  Pp  = . 3 mmHg 
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F-'igure  - Exp>er i mon ta  1 (Run  No.  11)  and  -^iniulated  drying  .urvps.  P = .1  minHg 


Figure  y.i')  - Experimental  (Run  No.  25)  and  simulation  drying  curves.  Pp  = 1 romHg 
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Figure  9.11  - Experimental  (Run  No.  2~ ) and  simulated  drying  curves.  P = .2  mmHg 


The  possibility  for  a non-uniform  field  is  '.ndicated  by  the 
local  vaporization  which  occured  in  Runs  No.  14  and  16,  and  the 
local  burns  of  the  sample  observed  in  Run  No.  13.  It  would 
explain  the  greater  discrepancy  obfierved  for  this  latter  run 
between  simulated  and  experimental  drying  curves  (Fig.  9.9). 

Indeed  the  sudden  decrease  of  slope  in  Fig.  9.9  of  the  experimental 
drying  curve,  with  respect  to  the  theoretical  one,  corresponds 
to  a peak  of  the  drying  rate  curve  (Fig.  9.6)  at  t = 1 hr.,  which 
may  have  been  caused  by  a local  melting. 

A non-constant  standing  wave  pattern  or  a non-uniform 
field  distribution  in  the  sample  may  have  caused  the  difference 
observed  in  Figs.  9.8,  9.9,  9.10,  and  9.11  toward  the  end  of  the 
dehydration.  Although  it  may  have  been  the  fact  of  the  non- 
uniform  thickness  (this  stands  particularly  for  the  pressure 
runs.  Figs.  9.10,  and  9.11).  Indeed,  it  is  not  possible  to  differ- 
entiate between  the  multiple  causes  of  non  ideality.  Given  the 
uncertainties  on  the  experimental  conditions,  the  general  agreement 
obtained  between  the  theoretical  and  experimental  drying  curves 
can  be  considered  as  good,  considering  the  small  differences  observed. 
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- Chapter  X - 


CONCLUSIONS 


1.  A one-dinensional  model  of  the  freeze-drying  process  using 
microwave  energy  has  been  derived  from  the  analysis  of  the 
transport  phenomena,  assuming  an  infinitely  sharp  sublimation 
front.  The  complete  transient  mass  and  energy  transfer  equations 
have  been  written  in  both  the  frozen  and  dried  regions.  The  trans- 
port and  dielectric  properties  are  allowed  to  vary  with  both 

time  and  location  in  the  sample  as  a function  of  temperature  and 
pressure.  The  adsorption  effect  of  the  water  vapor  in  the  dried 
layer  is  neglected.  Pick's  equation  is  used  to  describe  the 
flow  of  water  vapor  in  the  dried  layer . 

2.  The  model  enables  one  to  describe  the  sublimation  of  the 
frozen  moisture  content  before  and  after  the  microwave  power  is 
turned  on.  The  final  drying  stage,  of  little  industrial  signi- 
ficance, in  which  all  frozen  material  has  disappeared  and  the 
adsorbed  residual  moisture  (less  than  5%  of  the  initial  sample 
weight)  is  being  removed,  is  not  covered  by  this  study. 

3.  A numerical  solution  of  the  transient  energy  and  mass  transfer 
equations  in  both  the  frozen  and  dried  layer  is  proposed.  It  is 
based  on  an  approximation  of  the  p>artial  differential  equations 
by  Crank -Nicolson  method.  The  finite-difference  equations  are 
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■olvtd  by  Gausi'i  •llmlnation  technique.  An  iteration  oi  the 
interface  temperature  is  used  to  handle  the  energy  and  mass 
transfer  equations  which  are  coupled  at  the  interface.  A second 
iteration  loop  is  used  to  predict  the  interface  motion. 

4.  A new  numerical  method  is  prooosed  to  accelerate  and  control 
the  convergence  of  an  iteration  sche.r.e.  Use  of  this  method 
enables  one  to  reach  convergence  in  2 to  4 iterations. 

5.  The  results  of  the  theoretical  study  show  that  the  freeze- 
drying process  is  heat  transfer  limited.  Drying  rates  are  essen- 
tially a function  of  the  microwave  power  input.  The  total  essure 
and  the  partial  pressure  of  water  v.por  in  the  vacuum  mainly 
affect  the  temperature  within  the  sample  during  freeze-dry: ng  and 
have  little  effect  on  the  drying  time.  Low  pressures  should  be 
used  in  order  to  ensure  temperatures  as  low  as  possible  in  the 
sample  and  thus  allow  one  to  shorten  the  drying  time  by  using 

a higher  microwave  power  input.  However,  pressures  below  75  ^Hg 
are  useless  as  their  effect  becomes  insignificant  while  the  flow 
regime  in  the  pores  becomes  mainly  Knudsen  Diffusion. 

Ihe  initial  temperature  of  the  sample,  prior  to  Prying,  is 
not  an  important  parameter. 

6.  A reduction  of  the  sample  size  results  in  shorter  drying  times 
and  a greater  ou  put  of  the  freeze-dryer  as  a higher  microwave 
power  level  is  possible  without  melting  or  overheating.  An 
optimal  operation,  yielding  a maximum  freeze-dryer  output,  would 
correspond  to  an  operatnn  near  corona  and  overheating  (or  melting) 
conditions . 

7.  With  the  assumption  that  a corona  discharge  would  occur  in  the 
microwave  applicator  beyond  255  V/cm  (peak  strength)  and  tha'  the 
diffusivity  values  cf  turkey  breast  are  applicable  o beef  meat, 
drying  time  as  short  as  1^  hours  would  be  possible  for  a ^ in.  thick 
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-4  2 

6.  A value  of  the  heat  transfer  coefficient  (1.9  x 10  cal/cn  /sec/”C) 
has  been  detarnlned  from  the  start-up  portion  of  the  experinental 
drying  curves.  Use  of  this  value  in  the  prediction  of  the  drying 
curves  gave  an  excellent  agreement  between  the  theoretical  and 
experimental  start-up  portions  of  the  drying  curves. 

9.  Comparison  of  the  theoretical  (170  V/cm)  and  experimental 

(125  V/cm)  electrical  field  peak  strengths  at  which  melting  starts 
to  occur  in  the  frosen  core,  seems  to  indicate  th^t  the  actual 
diffusivity  in  the  experiments  may  be  somewhat  lower  than  that 
assumed  in  the  calculations.  However,  the  difference  may  also 
reflect  experimental  heating  non-uniformities. 

10.  Finally,  the  predicted  drying  curves  are  in  good  agreement  with 
those  obtained  experimentally.  The  small  differences  observed 

may  be  due  to  an  uncertainty  on  t.ie  value  of  the  diffusivity 
(or  that  of  the  dielectric  properties)  or  non-ideal  experimental 
conditions,  e.g.. 


1.  non-uniform  electric  field  distribution  in  th^  sample 

2.  non -constant  standing  wave  pattern  in  the  cavity 

3.  non-uniform  thickness  of  the  sample 


11.  The  mathematical  model  can  be  used  to  simulate  with  a good 
r.c7u^'acy  the  sublimation -dehydration  process  using  microwave  dielec 
■‘.iir  bea-^ing. 
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- Chapter  XI  - 


RECOMMENDATIONS 


1.  Future  improvement  of  the  mathematical  model  may  include  the 
moisture  adsorption  effect  in  the  dried  layer.  This  would  enable 
one  to  account  for  the  variations  of  the  dissipation  factor  of  the 
dried  product  with  the  sorbed  moisture  content,  as  experimental 
data  are  obtained.  It  would  also  permit  one  to  der-.ribe  the  final 
drying  stage  in  which  the  residual  sorbed  moisture  is  removed 
from  the  dried  product. 

2.  There  is  a need  for  more  experimental  data  on  the  physical 
properties  of  frozen  and  dried  beef  meat.  These  data  should 
include  especially  the  following: 

1.  the  effective  diffusivity  in  the  operating  p/essure 
range  (.05  to  1 mmHg) 

2.  the  dielectric  properties  in  the  temperature  range, 

-50°C  to  0“C  for  frozen  beef  and  -50°C  to  60°C  for 
dried  beef  as  a function  of  composition  (moisture,  fat, 
and  salt  contents) 

3.  Shorter  computer  times  in  the  mathematical  simulation  may  be 
obtained  through  a more  efficient  organization  of  the  numerical 
scheme  and  the  determination  of  an  optimal  time  increment  size 
strategy.  The  latter  would  require  an  extensive  stability  analysis 

4.  More  extensive  and  comprehensive  experimental  data  on  the 
microwave  freeze-drying  process  are  still  needed.  There  is  a need 
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lor  a systematic  experimental  study  of  the  effects  of  the  process 

constraints  and  variables  on  the  drying  time,  in  order  to  support 
the  theory. 


S.  Special  care  should  be  used  in  the  design  of  an  experimental 
freeze-dryer  in  order  to  ensure  conditions  of: 

1.  corona -free  operation 

2.  uniform  electric  field  distribution  in  the  sample 

3.  constant  standing  wave  pattern  in  the  microwave 
applicator 

4.  in-proc  .ss  determination  of  the  electric  field  strength 
in  the  vacuum  near  the  sample  and/or  the  microwave 
power  loss  in  the  sample 


appendix  a 


transfer  equation  in  the  dried  layer: 


And: 


p<--  w-/ 

ftCx,  k'f(TA-To) 

('/>./.  /o) 

r^./.  n) 
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A. 2 ENERGY  BALANCE  AT  THE  INTERFACE 


An  tntrgy  b«l«nc«  it  ntdt  on  t slab  ont  space  increment  thick,  centered 
at  the  interfacei 

/ ^jcx3ri\/ecto/\^^  j r l^c(UitritUatioh^^l^o^u^t'^i^onsum^f>pyil 

where: 


Icohdi/doyf  iiijs.  - 4 /tiiS 

lx 

(ccnvtdiofn^^l^l/fi  Cf>^  ^ ^ 

( 3 ^ 

lconS,Uyt1pi9>)/m  yfc 

ll>rociuc^^j  ^ 

ItucwmM^  Is  (^Q>  ^ Ip  ^ 

Substitution  of  the  transformation  and  normalitation  equations 

(see 

Chapter  IV,  Sec.  4)  into  Equations  (A. 2.1),  (A. 2. 2),  (A. 2. 3), 
(A. 2. 5),  and  (A, 2. 6)  gives 

(A. 2. 4), 

CO^/D  s.  ^ Oii  J~  iSf . ^ ^^0 

s /-s 

J') 

cnuft  - TV;  i.  -aw, 
k^CTi^To) 

(&.^.^) 

fROD  r f Cu^ChS)  ] 

{^A)0) 

-teuffs  /|.«I  i)»>  r ide  TVZ  ' 

^ ^ htr  ^ p<TOi'p(jP^  ^ 

4j  ^ 

with 


1/:  ^ fdr\ 

' " u/’S) 
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The  energy  balance  becomes 


CONO  e ^CtOH  ~ PAcC  -f  CONS  (fi.-O  .l2>) 

where  the  terms  COND,  CONV,  ACCUM,  PkOD,  CONS  correspond  respectively 
to  the  terms  of  conduct'  convection^jj^,  accumulation,  production 

(microwave),  consumpti  'sublimation)  in  the  energy  balance  athe  inter 
face. 


The  corresponding  finite-difference  equations  are: 
COND^  z.  ^ <Xp  i /}i^\  __  I 


COhjV 


_ 'IM'  ^ 

?1  - - ~o— 


(f1 


CG'VC..  _^J= 


’ ” , o / ■ 


/'  v; 


_ L-  j£o£o_ 

“ I Pf.  Cp 


w i t h 


.V, 


(" ^rp-.  ^ J /d^ 


L(  f-  ) 


Typical ly , 


iMf 


(^.QQl) 
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Calculation  of  the  interface  temperature 

The  terapcrature  at  the  interface  at  which  would  have  satis 

fied  the  energy  balance  (as  upposod  to  that  guessed)  is  given  by^ 


ft  -*•  -t  PAOOmtl/^  - CqajS^^^‘^&C2  mtj/4 


Calculation  of  the  interface  position: 
From  Eq . ^ 5 . 14 ) : 


^ ^ * »S 


ifi.Ji.  -J?3 ) 


a-  ACC2  in  Eq.(A.2.22)  corresponds  to  the  sum  of  the  two  latter  terms 
in  the  right  hand  side  in  Eq.  (A. 2. 18).  AINT  in  Eq.  (A. 2. 22)  represents 
the  coefficient  of  (^®^)p+3^  i"  Eq.  (A. 2. 18). 
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A. 3 BOUNDARY  CONDITION  AT  THE  OUTER  SURFACE 


An  energy  balance  at  the  outer  surface  of  the  slab  (dried  layer  side) 

Is  used  to  derive  the  boundary  condition  for  the  temperature  as  of  Eq.(5.10). 
The  energy  balance  is  made  on  a slab  half  a space  increment  thick: 

iz  ^ It  ^ 


Substitution  of  the  transformation  and  normalization  equations  in 


Eq.  (A  3.1)  and  use  of  the  finite  difference  approximation  at 'C' 


will  lead  to  the  desired  form  for  the  boundary  condition. 
Equation  (A.3.1)  becomes 


n+?5 


■//.//  A?  1-L  cJ*  / cis  ABt  __  /,,s) 


Let : 


isM)  f..r  a%L‘>0-&)  <ii 
pfCf.  a « 


II 

(A,^.  a ) 

^ ^ ^ Sn^i  - Sn  j / PdC^ 

b Pf.Cp 

(P.3M) 

Off  _ /•/vj'f  cop  L^C/~Sn)  6"h 

Then,  the  finite  difference  equation  corresponding  to  Eq.  (A, 3. 2)  is 

a (A2>.s) 
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Equation  (A. 3. 8)  corresponds  to  an  energy  balance,  In  the  fixed 

boundary  domain,  over  the  slab  (2,..  , 2.)  between  times  tT  and  C 

1 n '^n+1 

It  Is  rearranged  to  give  the  final  form  of  the  boundary  condition  as 
required  by  the  numerical  scheme  of  solution  (Eq.  (5.10)),  One  has 


a 


q) 


where 

-f  ^ J ] / ■+  _ iS /Kn+f^  ) ( H ) 
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APPENDIX  B 


B.l  SOLUTION  OP  THB  SYSTEM  OF  FINITS-DIFPBRBNCE  EQUATIONS 


The  finite-difference  equations  of  Appendix  A art  i'uitten  under 
the  general  form 


^i  " ^ Vl  * Pi  « 


i+2  i 


i = ] 


1 ) 


while  the  boundary  conditions  are  written  as; 
at  the  outer  edge  (center  line  or  outer  surface) 

\ * ®i  “2  * \ 

and  at  the  interface  (ice-front) 

D U + U ■ q 
P ^I  “i+i  ^ 


(B.2) 


(B.3) 


They  form  a set  of  I+l  linear  equations  which  are  to  be  solved  with 
I+l  unknowns,  the  U^,'s  i*l,...,I+l.  One  has 


U1.81U2 


Ui.aiU2*PlU3 


* K. 


= t. 


^i.i"“i.iVPi-i^'in 

Ui*  - Pi 


--i 


i-l 


ti 


U +a  U+B  U 
I-l  I-l  I Pl-1  I+l  I-l 


p Uj  + 


^I+l  = ^ 


A forward  elimination  of  the  lowest  subscript  unknown  U^,  starting 
i*l  yields  the  8^*s  and  the  *s  defined  by 


“i  * «i  “i*i  ■ 


i— 1 , . . ■ , I 


(B.5) 
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it  <:on®l^ 


for  i*2 

A backward 
ii;  Eq.  (B.S 


a , ,*  '(!>.  , 
i“l  1-1 

(B.6) 

>1. 

X 

y - ■ 

, '■  i-1  '.i-i 

(B.7) 

i-1  ,\-J 

r 

, A 4 

and  are  given  o>  the  L 

ic.  undr-ry  condition,  Eq. 

substitution  of  the  nuaierical  , iljes;  c.t  the  8's  ancle's 


),  starting  i=I,  yields  the  u-iknov.nj 


q - p 


I+J  1 - p Cj 


The  corresponding  PORTRAN  IV  compute  preoraip 
(subroutine  EMMA). 


(B.8) 

I (B,9) 

i:-  g\  Cl’  is  App.?ndix 


F 
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B.2  FORCING  AND  ACCELERATING  THE  CONVERGENCE  OF  AN  ITERATION  SCHEME 

ThtoraticAlly,  Itarttiva  procadurti  aneountarad  In  tha  solution  of 
nathanttlcal  problami  corratpond  to  the  solution  of  tha  aquation 

X •^{x)  (B.IO) 

by  an  organized  method  of  trial  and  error.  The  variable  x In  Equation 
(B.IO)  may  not  be  a scalar  variable,  while  the  functional  ^ Is  not 
necessarily  known  explicitly.  may  for  Instance  result  of  the  solution 
of  a more  complex  problem  for  a given  value  of  the  "parameter"  x. 

In  ord^r  to  solve  the  above  equation  a value  x*  of  the  variable  x 
Is  "tried".  The  corresponding  value  ^/(x*)  of  the  functional  Is  then 
obtained  and  compared  to  that  of  x according  to  a predetermined  con- 
vergence criterion. 

The  usual  Iteration  scheme  Is  such  that  , If  the  convergence 
criterion  Is  not  satisfied,  the  current  value  of  the  functional 
Is  fed  back  as  the  lext  trial  for  the  value  of  x.  This  procedure  Is 
repeated  (Iterated)  until  the  convergence  criterion  Is  satisfied,  l.e. 
until  the  difference  |x  - '|^(x)[  Is  less  than  a given  number.  Mathemat- 
ically, the  Iteration  scheme  Is  given  by 

(B.ll) 

where  x^  represents  the  Jth  trial.  Geometrically,  this  Iteration 
scheme  correspond  to  the  search  of  the  Intersection  of  the  curve 
y ® 'l^(x)  with  the  first  bisector  y * x In  a coordinate  frame  of 
abscissa  x and  ordinate  y (see  Fig.  B.2.1).  Unfortunately,  this 
scheme  is  not  always  convergent.  Indeed,  In  many  Instances  It  may 
converge  very  slowly,  resulting  In  long  and  therefore  costly  machine 
calculations.  In  effect,  theoretically  this  scheme  converges  only  If 
the  slope  of  the  curve  y = •>p(x),  at  the  solution  point.  Is  less  than  1 
In  absolute  value: 

(B.12) 

An  alternate  scheme  Is  proposed  here  to  force  and  accelerate  the 
convergence  of  an  Iteration.  This  scheme  Is  such  that  the  next  trial 
Is  computed  from  the  previous  trial  and  the  previous  value  of  the 


Figure  B.2.1-  Geometrical  representation  of  the  iterat 


ion  schemes. 
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(B.13) 


functional  ty  according  to 

= (1-^)  X*  + 

where  the  parameter  ^ , which  will  be  referred  to  as  acceleration 
coefficient  may  be  adjusted  to  force  and  accelerate  the  convergence 
of  tiic  iteration  scheme.  Geometrically  (see  Fig.  B.2.1),  this  scheme 
correspond  to  the  search  of  the  intersection  of  the  curve 

y = (l-15)x  + ^t(x)  (B.14) 

with  the  first  bisector  y = x.  The  change  is  possible  since  the  two 
functionals  ^(x)  and  (l-5)x  '^(x)  have  the  same  value  at  the  searched 

solution  (c  nvergence)  point  where  x ='|^(x). 

In  theory,  practically  instant  convergence  is  reached  with  this  new 
scheme  if  the  slope  of  the  curve  of  equation  (B.14)  is  equal  to  zero  at 
the  solution  point  (see  Fig.  B.2.1).  The  main  advantage  of  the  method 
is  that  the  slope  of  the  latter  curve  (B.14)  can  be  adjusted  and  forced 
to  a near  zero  value  by  varying  the  parameter  , thus  yielding  extremely 
fast  convergence. 

The  optimal  value  of  15  is  given  by 

0 = (1-Ux  +1§q>(x)  = l-'5  +5  ^ (B.15) 

Thus 

^op,.  = 

Use  of  this  method  in  the  calculation  of  the  ice-front  temperature 
by  iteration  in  the  simulation  permitted  to  reache  convergence  in  2 to 
4 iterations  on  the  average.  The  acceleration  coefficient  was  cons- 
tantly updated  during  a given  iteration.  It  should  be  noted  that  the 
iteration  by  the  first  scheme  (direct  feed  back  of  the  computed  value) 
was  highly  divergent. 
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- appendix  c • 


SAMPLE  CALCULATIONS  (RUN  No.  15) 


C . 1 Calculation  of  h 

For  drying  on  both  sides  (symmetric  slab),  Eq.  (9.5)  gives 


AH  ^ 


(C.l) 


with 


( C2) 


From  Fig.  9.3  at  t * t^ 

'6,2  g/hr  * 1.7  x io‘^  g/sec 
Approximate  A^  by  A^.  Thus  from  Appendix  E 

Ag  “ 61.5  cm^ 


Assume  a,  room  (vacuum  chamber)  temperature  of  20'C.  Approximate 

by  the  frost  temperature  of  the  vacuum  chamber  at  = .3„n,Hg 
(cf.  Table  9.1).  Then,  from  the  ice-water  vapor  equilibrium  curve 

^i,l  * -29.5'C 

rrom  m^ , m^,  m^ , andL^  (cf . Appendix  E)  determine 
e m^ 

■ m - m^  ) 4 ' -936)(.78)=  .05  cm 

O X 

where  L is  approximated  by  L . 

o 

From  Fig.  *.4,  at  P^^  = .3mmHg  one  has 

•4 

ko  = 10  cal/sec/cm/®C 
From  Table  6.1 

= 675  cal/g 
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Substitution 
h ■ 2 


of  the  numerical  values  in  Eq. 
1 X 10  ^ cal/sec/cm^/®C 


(C.l)  gives 


C . 2 Calculation  of  E 

From  Eqs.  (9.1)  and  (9.2)  one  has 


From  Table  6.1 

=675  cal/g 

Assume  an  average  = - s*C.  Then  from  Fig.  6.5 
Kp  * 4.  X lo"^  cal/sec/cm^/(V/cm)^ 

Assume  an  average  T^^  = lO'C.  Then  from  Fig.  6.6 
Kjj  * .9  X 10‘^  cal/sec/cm^/(V/cm)^ 

The  volumes  and  V^^are  determined  by  (symmetric  slab) 

V =SV°«2SL  A 
^ F O o 


(C.3) 


(C.4) 


% * \ (C.5) 

L,  A are  given  in  Appendix  E The  overall  drying  rate  m (drying) 
on  both  sides)  is  obtained  from  Fig.  9.3.  The  electric  field  is  then 
calculated  for  two  locations  of  the  ice-front  (interface). 

At  S = .5  : 


Thus 


m = 20.6  g/hr  = 5.7  x lo’^  g/sec 
Lq  = .78  cm;  = .7  cm 

^o  " ~ 

V°  = 2 X ,78  X 61.5  = 95.7  cm^ 

= 2 X .70  X 51.0  = 71.3  cm^ 
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X .5  X 71.3 


■ 673  X 3.7  X 10'^ 

4 X i0“®  X .3  X 93.7  + .9  x lo"® 

■ 1.8  X 10^  (V/cm)^ 

E ■ 134  V/cm 
Similarly  at  S ■ .3  ; 

n ■ 11.8  g/hr  * 3.28  x lo‘^  g/sec 

and 

■ mf, 675  X 3.28  x lo'^  ^ 

4 X 10  X .3  X 96.7  + .9  x lO'^T 

* 1.38  X 10^  (V/cm)^ 

E - 117  V/cm 

Average  field  peak  strength, 

E * ij(134  + 117)  * 125  V/cm 


.7  X 71.3 
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APPENDIX  D - 


COMPUTER  OUTPUT 


The  computer  output  pertaining  to  this  investio-’ition  are 
in  the  possession  of  Professor  Yi  Hua  Ma  of  the  Department  of 
Chemical  Engineering,  Worcester  Polytechnic  Institute,  Worcester, 
Mass.  01609.  They  are  kept  in  Room  225. 


- APPENDIX  B - 


EXPERIMENTAL  RESULTS 

The  experimental  results  are  known  with  the  following  ac- 
curacies (on  an  absolute  basis): 

Am  = . 2g 
AL  = . 2cm 
AA  = . 5cm^ 

AFp=  .03-.09mmHg 

AP’^=  .OlmmHg 
H 

The  power  measurements,  which  could  not  provide  quantitative 
information  due  to  calibration  problems,  and  the  start-up  portion 
of  the  drying  curves  are  not  reproduced  here.  They  are  in  the 
possession  of  Professor  Y.H.  Ma  of  the  Department  of  Chemical 
Engineering,  Worcester  Polytechnic  Institute,  Worcester  MA  01609. 
Room  225. 
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RUN  No  25 


n 

o 

» 73.0  g 

= 26.2  g 

L 

0 

= .80  cm 

~ .70  cm 

A 

0 

* 49.0  cm^ 

= 45.3  cm^ 

"r 

m 

hr 

mmHg 

mmHg 

9 

0 

.24 

1.0 

70.3 

.5 

.45 

I .0 

63.7 

1.0 

.45 

1.0 

55.6 

1.5 

.41 

1.0 

47.3 

2.0 

.31 

1.0 

40.9 

2.5 

.24 

■ .0 

36.3 

3.0 

.20 

.9 

33.2 

3.5 

.16 

.9 

30.9 

4.0 

.13 

.9 

29.0 

4.5 

.11 

.9 

27.7 

5.0 

.09 

.9 

26.8 

5.5 

.07 

.9 

26.2 
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RUN  No  27 


m = 
0 

78 . 0 g 

- 25.4  g 

L = 
0 

.80  cm 

.75  cm 

A = 
0 

48 . 7 cni^ 

= 45.0  cm^ 

■“r 

Pr 

m 

hr 

DiDiHg 

mmHg 

9 

0 

.21 

. .21 

75.3 

.5 

.23 

.23 

68,7 

1.0 

.21 

.21 

60.5 

1.5 

.20 

.20 

51.6 

2.0 

.18 

.18 

43.9 

2.5 

.16 

.16 

38.3 

3.0 

.15 

.15 

34.9 

3.5 

.14 

.14 

32.3 

4.0 

.13 

.17 

30.2 

4.5 

.12 

.16 

28.7 

5.0 

.11 

.12 

27,3 

5.5 

.11 

.12 

26.4 

6.0 

-- 

— 

25.7 

6.5 

-- 

- - 

25.4 
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RUN  No  29 


■ 74,2  g * 26,3  g 

Lq  ■ .75  CO  * .70  cm 

■ 49,6  co^  * 39.1  cm^ 


Pr 

m 

hr 

omHg 

mmHg 

9 

0 

. _ 

3.4 

73.3 

.5 

.55 

3.5 

65.0 

1.0 

.52 

3,7 

52.4 

1.5 

.42 

3.7 

41.1 

2.0 

.39 

3.7 

31.3 

2.5 

m m 

m m 

27.6 

3.0 

m m 

m m 

26.6 

3.5 

-- 

m m 

26.3 

r'"  ^ 

RUN  No  12 


rn 

0 

- 74.5  g 

= 22.1  g 

L 

0 

= .71  cm 

= - 

A 

0 

= 58.6  cm^ 

Af  = -- 

P" 

R 

'^R 

m 

hr 

mmHg 

mmHg 

9 

0 

.23 

.32 

72.1 

.5 

.29 

.29 

66.7 

1.0 

.26 

.26 

60.4 

l.S 

.23 

.23 

54.9 

2.0 

.21 

.22 

50.3 

2.5 

.20 

.22 

46.5 

3.0 

.19 

.20 

43.0 

3.5 

.18 

.20 

39.8 

4.0 

.17 

.20 

36.9 

4.5 

.16 

.19 

34.1 

5.0 

.16 

.18 

31.7 

5.5 

.15 

.18 

29.3 

6.0 

.15 

.17 

27.2 

6.5 

.14 

.17 

25.3 

7.0 

.12 

.17 

23.9 

7.5 

.12 

.17 

22.9 

8.0 

.10 

.17 

22.3 

8.5 

— 

- _ 

22.1 
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RUN  No  13 


m = 
o 

88,0  g 

"Jf  = 26.2  g 

L = 
o 

• 74  cm 

= - 

A = 
o 

— 

Af  = - 

m 

hr 

mmHg 

mmHg 

g 

0 

.27 

.29 

84.9 

.5 

.20 

.26 

77,2 

1.0 

.25 

.21 

69.3 

1.5 

.22 

.21 

61.1 

2.0 

.21 

. 2U 

54.3 

2.5 

.18 

.20 

49.4 

3.0 

.17 

.18 

45.2 

3.5 

.16 

.17 

41 . 7 

4.0 

.15 

.17 

38.4 

4.5 

.15 

.17 

35.3 

5.0 

.14 

.;6 

32.7 

5.5 

.14 

.16 

30.3 

6.0 

-- 

-- 

28.5 

6.5 

-- 

-- 

27.1 
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RUN  No  14 


in 

0 

■ 85.6  g 

in^  * 26.0  g 

L 

0 

• .71  cm 

1 

1 

N 

A 

0 

■ 60.5  cm^ 

> 

H 

1 

1 

Pr 

m 

hr 

lamHg 

mmHg 

9 

0 

.23 

.24 

79.9 

.5 

.26 

.26 

65.5 

1.0 

.23 

.23 

48.5 

1.5 

.19 

.19 

39.3 

2.0 

.17 

.17 

34.2 

2.5 

.16 

.16 

30.9 

3.0 

.15 

.15 

28.3 

3.5 

.14 

.14 

26.6 
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RUN  No  IS 


in 

0 

= 87.5  g 

= 25.3 

9 

L 

o 

= . 78  cm 

o 

II 

cm 

A 

o 

~ 61.5  cm^ 

= 51.0 

2 

cm 

’‘I 

Pr 

m 

hr 

mmHg 

mmHg 

9 

0 

.30 

.30 

^ 83.2 

.5 

.30 

.30 

73.1 

1.0 

.28 

.28 

61,8 

1.5 

.24 

.24 

51.3 

2.0 

.21 

.21 

44.1 

2.5 

.19 

.19 

39.1 

3.0 

.18 

.18 

35.2 

3.5 

.16 

.16 

32.1 

4.0 

.15 

.15 

29.5 

4.5 

.14 

.14 

27.5 

5.0 

.13 

.13 

26.1 

5.5 

.12 

.12 

25.4 

-l4l- 
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RUN  No  16 


ID 

0 

■ 87.5  g 

* 27.0  g 

L 

0 

■ .76  cm 

S - -- 

A 

0 

.c,  « 2 

■ 61.0  cm 

II 

1 

1 

Pr 

m 

hr 

mmHg 

mmHg 

9 

0 

.30 

.31 

83.7 

.5 

.32 

.32 

72.4 

1.0 

.30 

.30 

58.5 

1.5 

.25 

.25 

46.8 

2.0 

.21 

.21 

40.5 

2.5 

o 

CM 

.20 

35.9 

3.0 

.19 

.19 

32.5 

3.5 

.16 

.16 

30.0 

4.0 

.15 

.15 

28.0 

4.5 

.14 

.14 

27.1 
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I 


RUN  No  17 


m 

0 

■ 94.0  g 

lUj  ■ 28.3  g 

L 

o 

“ .71  cm 

“ .62  cm 

A 

o 

“ 60.8  cm^ 

= 56.5  cm^ 

-"r 

"r 

m 

hr 

mmHg 

mmHg 

9 

0 

.27 

.26 

86.3 

.5 

.22 

.22 

77.5 

1.0 

.20 

.20 

68.4 

1.5 

. 19 

.19 

59.6 

2.0 

.18 

.18 

51.8 

2.5 

.17 

.17 

45.4 

3.0 

.16 

.16 

44.4 

3.5 

.16 

. 16 

36.5 

4.0 

.15 

.15 

33.2 

4.5 

.14 

.14 

30.6 

5.0 


.13 


.13 


28.7 


oooooooooooooooooooo 


- APPENDIX  F - 


SUBf^OUTlNE:  EMMA 


THIS  SUBROUTINE  SOLVES  THE  FINITE  DIFFERENCE  EQUATIONS 
FORMING  A SET  OF  M LINEAR  EQUATIONS,  THE  METHOD 
USED  IS  similar  to  THAT  KNOWN  AS  GAUSS'S  ELIMINATION 
TECHNIQUE,  SEE  APPENDIX  C,1  DF  PHD  THESIS  FOR  SIGNIFI- 
CANCE OF  arguments  of  Subroutine, 
note  that  m in  the  subr.  arg,  has  to  be  lt  2e0 

subroutine  EMMA(M,DE1|XLA1|P#OiU) 

COMMON  N,TImE,AL(199) i BE ( 199 ) , G A ( 199 ) 

DIMENSION  HE ( 2Ca ) ,XLA (200)10(201) 

DE(1)*DE1 
XLA(  DsXLAI 
DO  1 Ii2,M 
DENiAL( I-1)-0E<  I-l) 

IF(DEN,EQ,P, )G0  TO  3 
DE(  I )«RE(  l-D/OEN 

1 XLA(I  )>(GAM-1)-XLA(  I«l)  )/OEN 

IF(P*DE(M) ,E0,1. ) GO  To  3 

U(M^l)t(OeP#XLA(M) ) / ( 1 , -P*DE ( M ) ) 

CO  2 I«l,  -I 
J«M+1«I 

2 U( J) sXLA( J),QE( J)#U( J*l) 

RETURN 

3 WRITE(5i  200)  n.MF 

200  FORMATdX, 'DIVISION  PY  ZERO  IN  SUBRl  AT  TIME  I '|F7,2) 
PAUSE 
END 
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OO.  OOOC,  0000000000000000000000000000000000000000000000000000 


subroutine 


MWTD 


THIS  SUbROUTiNE  SIMULATES  THE  E HEE2E-DEHYDP AT  I ON 
PROCESS  USING  HlCROwAVE  OlELECTHlC  HEATING,  IT 
IS  DIHECTly  applicable  to  the  EREE?E-DRYlNr,  OF 
beef  meat  (A  OIEEEPEnT  type  of  material  to  re  dried 
HILL  necessitate  A CHANGE  OF  THE  PHYSICAL  DATA 
OF  BEEF  MEAT,  ASSUMED  HERE)! 

ENTER  OPERATING  CONDITIONS  THROUGH  ARGUMENTS  OF  SUB. 
ROUTINE  ACCOROINU  TO  NOMENCLATURE  AND  UNITS  BEl.OW, 

FIELD!  ELECTRIC  FIELD  PEAK  STRENGTH  IN  VACUUM  INSIOE 
THE  CAVITY,  V/CK, 

PK:  total  PRESSURE  JN  VACUUM  CHAMBER,  MM-HC, 

PHRi  Partial  pressure  of  water  vapoR  in  vacuum  chamber 
mm-hg  , 

TR;  AMRIanT  temperature  IN  CAVITY,  DEG  K, 

TB!  initial  temperature  OF  FROZEN  PRODUCT,  DEG  K, 
thick:  product  slab  ThICKI.ESS,  cm, 

WATER!  initial  PRODUCT  MOISTURE  COnTEnT-wET  BASIS  (NOT 
IN  )(),-. 

FATI  fat  content-wet  basis,  S-, 

MElTi  returns  TO  MAIN  PROGRA'-'  1 IF  MElTIN'G  OCCURS,  OTHERe 
WISE  EQUAL  TO  ZERO, 

general  NOMENCLATURE! 

OK(P)l  Thermal  conduct! VlTYaORIED  PRODUCT,  CAL/CM/SEC/ 

DEG  k;  ARGUMENT!  LOCAL  "RESSURFi  MM-hC, 

FK(T)!  thermal  CONOUCTl VITYbFRO^EN  PRODUCT,  CAL^CM/SEC 
/DEG  Ki  ARGUMENT!  LO"AL  TEMPERATURE !DEC  K, 

D(P)I  EFFECTIVE  DIFFUSIVITY,  CmZ/SEC, 

DL(T)!  dissipation  COEEF, -DRIED  PRODUCT,  CAL/Cm3/SEC/ 
(V/CM)i  ARGUMENT!  LOCAL  TEMPERATU»E  IN  DEG  K, 
EL(T)i  DISSIPATION  CCEEE , sEROjIEN  PRODUCT,  DAL/CM3/ScC/ 
(V/CM)>  argument  I temperature,  DfG  K, 

EE(T)!  ICE-VAPOR  equilibrium  RELATIONSHIP;  ENTER  EOU. 
temperature  in  deg  K,  returns  EQU,  CDNCEf TRATION  OE 
WATER  VAPO°  IN  GR/CM3, 

wcapi  specific  heat  of  water  vapor,  cal/gr/oeg  k, 

OCaPi  specific  HEaT  or  DRIED  product,  CAL/GR/DFG  k, 

FCaPI  specific  heat  OE  EROFEN  PRODUCT,  CfL/GR/pEG  k, 

DEnSI  DEviSITV  OE  ICE,  GR/Cm3, 

EDENS!  bulk  density  OE  EROFEN  PRODUCT,  i,h/CM3, 

ODENS!  bulk  density  OF  DRIED  PRODUCT,  GR/CM3, 

htc  I HEAT  transfer  coefficient  at  ext.  DUREACFi 

CAL/CM2/SEC/DEG  K, 

HSi  HEAT  OF  Sublimation  of  ice,  cal/gr, 

CTPi  water  vapor  concentration  at  triple  point,  GR/CM3 
CR!  AMBliNT  water  VAPqR  CONCENTRATION,  GR/CM3, 

PGCI  PERFECT  GASES  CONSTANT,  MMaHG/DEG  K/(GR/CM3), 

MM!  number  of  space  INCREMENT-EROEEN  CORE!  MUST  BE 

c less  (OR  euual'  than  2/a, 

c KKl  NUMBER  OF  SPACE  I NCREMENTS-UR 1 ED  L^YERi  MUST  bE 

-lUS- 


LESS  (OR  EQUAL)  THAN  2l22. 

Euri  microwave  enie«cy  usage  factorj  ratio  or 
total  energy  nEed  roR  sublimation  / 
total  microwave  energy  CONbUMfO, 


subroutine  Mwroir ieldiPRiPwh;TRi Tifl, Thick  I water, faTimelt) 

COMMON  N.TIML, Al(199) iBE(199) ,GA(199) 

DIMENSION  CN1(21) ,TNnlt21) ,TNri(21) 
dimension  rN2(2l) ,TND2(21) ,TNF2(21) 

DIMENSION  FC(2iai) . DC  (2?l)  ,011(201)  ,012(201)  ,FP1(2k)1)  ,FP2 
1(20:),DP1(2.-^1)  ,np2(2k)l ),  or  1(201  ),riF2(  2 Jl)  ,001(201)  ,002  ( 
2201) 

INTEGER  EPSI 
real  LAI 

FUNCTIONS  plOCK: 

DK(P)fl,F.4*(i,i26*ALOGtP)*(B,13‘>l4ALOG(P)*(O,01P65»ALOG 
1(P)*(-3,137Fo3^ALOG(P)*(»5,403E-A) ) ) ) ) 
FK(T)«AK*8k*S3RT(271.65-T) 

D(P)a23,l/(l,344P) 

rL(T)«l,ro6*(294,7^T*(33.62a*T*(l,464E-2-1.922E-S*T)  ) ) 
OL(T)  *l,E-6*(14,39*T*(-0,l577*T#(5,5,77r»4-c  ,924E-7*T)  ) ) 
FE(T) ■EXP(-63.23*0,2969*T-4,03eE-4*T*T ) 

ENTER  CONSTANTS! 

CALL  0FIlE(2J, 'PRINTI' ) 

CALL  CFILE(21,  'hIST* ) 

CALL  0FIlE(22, 'PROF*  ) 

CALL  0FIlE(23,  'PRINT?'  ) 

CALL  0FIlE(19,  ’CHECK' ) 

DATA  WCAP,FCAP,0CAP,0£NS,rDENS,DDENS/,5, , 4 ^ , 36 , , 9? , , 96 , 
1,32/ 

DATA  HTC,HS.CTP,PGC/1,9E-4,675, ,4,84E*6,3,46E3/ 
data  KK,hh/20,23/ 

AKsl,E-3*(2.9i5a2?.#FAT) 

fcK*l,Ee3»(C,i'2*i,8*(wATER«r  ,65)  ) 

POR*wAT£R#FUEr,S/PEf'S 

CR»PwK/Pr,c/TR 

FC0sFK( TC) 

FU0*FCi?/FOENS/FCAP 

FP0*FL(T(')*FIELO**2. 

DELi!al,/FLOAT(KK) 

OELY«1,/FlOAT(HM) 

DELT«0,072 

CETAU»0ELT*F03/THICK*«2. 

cphtc*thick/fc? 

KleKK*! 

Al«FDENS*rCAP« ( TRaT0  > /DENS/POR/HS 
A2*DDENS*0CAP/FDE’!5/FCap 
A3*FP0*ThICK*»2,/FC0/ (TR-T0) 
A4«0ENS*WCAP«P0R/DDENS/0CAP/0EL2 
A5*0. bawCAPaTHlCKPDELi/FCB 
A6»HS*THICK/FC0/(TP-T0) 

A7«ThICK/OENS/POR/FD0 
Atj«ro(fl«(rTpBC-n/THiCK 
AV«0ENS«WCA**P0R/F0ENS/FCAP* ,9075*A2 

PREPRINT  ASSUMED  EXP.  CONDITIONS! 

WRITE (20,200)  FIELD,PK,PWR,TR,T0,THlCK,wATrR,FAT 
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WRITE (23,2110)  riELD#PR|PWR|TR,T0,THICK,WATERiFAT 

ENTER  ALL  INITIAL  VALUES* 

TIMEt0, 

Ni0  ' 

S0*1, 

Sl«, 999*30 

TEMKERaTURE  PROriLES! 

TNS«1,/(1,*OK(PR)/HTC/ThICK/(1,sS1) ) 

DO  1 

TNri(MJs0, 

00  2 Kal.Kl 

TNDKK)  sTNS*(l.srLOAT(Ksl)/rLOAT(KK)  ) 

TN02(K)aTNDl(‘<) 

GRTOl*  ( ’*,01(KK*1)«TnD1(KK)  )/DEL2 

GRTris( TMFl(MM*l)-TNrl(MM))/DeLY 

DRYING  Rate  from  energy  PALANCE  AT  INTERFACE! 

C0EFl*A6-A5*GRTDl3.25*A7*(A2*(2,-DEL?i*DEL?*GRTDl"(2,8DE 

lLY)*UELY*r,RTFl) 

C0N01»bDk (PR)/FCP/ ' 1 , "Si )*GRT01 
DR1«C0NQ1/C0EF1 
concentration  PROFILE! 

GRC1*0R1*(1 ,-Sl)*THICK/(CTP-CR)/D(PR) 

00  3 K*1,K1 

CN1(K)»GRC1*FL0AT(K-1) /FLOAT (KK) 

CN2(K)=CN1(K) 

heat  consumption  at  interfaces 

0Il*C0M0l8(C0EFl-A6)*0Rl 

DELQI«5', 

miscell:anuous» 

Tiisa, 

OELTI10, 

XI0,75 

SF1*S1*S1 

SUls(l.-Sl)*(l,-Sl) 

FIELO1»0, 

MELT*0 
TMWE*0 , 


simulation  loo^-head: 

CONTINUE 

TIMEoTIME*OELT/3600. 

MN*1 


UPDATE  all  PROPERTIES! 
NTCs1,E2*(T0»TNF1(D*(TRsTP')  ) 
IFCNTO.LE, 27165)00  TO  A0 
MELT*1 

WRITE(5,101)T;MEiN 

return 

CONTINUE 

FIEL02*FIEL0 

IF(N,Lr,71O)FIELD2«0. 

AF1«a3*UETAi'/FP0*FIELD1*FIELD1 

AF2aA3*DETAU/FP0*FIELU2*FJELD2 

A01»AF1/A2 

A02»A'’2/A2 

DO  6 Mai, Ml 


TFaT0*TNFl(M)*(TR3T0) 
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rc(M)«FK(Tn/fC0 

FLTFaFLC^n 

FPKM)  »FlTF*'IElD1«FIELD1/FP0 
FP2(M)3FlTF*‘-  IElD2*FIELDJ?/FP0 
0F1(M)«FlTF*AF1 

6 QF2(M)=FLTF*Ar2 
DO  7 K»1,KI 
TD»T0*TND1 (K) *( THaTa) 
PKE«PR«PwR*PGC»(C"*CN1(K)*(CTP-CR) )«TD 
DC{K)iDK<PRE)/rca 
DI1(K)sD(PRE)/FO0 

DLTO«DL<TO) 

DPl{K)iULTO*FIELDl*n6LOl/FP0 
D?2(K)aOL  TD*F I ELD2*F I ELO2/FP0 
QDl{K)«DLTn»A01 

7 0D2(K)aDLTD»A02 
C 

C INTERFACE  POSITIONI 

OIGaaii^DELQI 

ITERal 

70  CONTINUE 

QI«0,9*<tJll4QIG) 

S2«S1-DETAU*A1»QI 

C 

C EXT,  SURFACE  BOUNDARY  CONDITIONI 

CDl-(DC<l)4DC(2n/2,/(l.eSl)/DEL£ 
CD28(DC(l)*0C(2n/2./(l.-S2)/DELZ 
CV«,5*A9#(R2rSl)/0ETAU 
ACa»25»A2*nEL2/0ETAU*(2.-S2-Sl> 
PRlr,5#A3#nELzi*(l.aSl)*(  ,25*DPl(2)*,75*DPia>  ) 
PR2a,5»A3#DEL2*(l.-S2)*( ,25*DP2(2)4 ,75*DP2(1) ) 
DEla(CD2-CV-,5*AC)/(-C-C02*CVel,5*AC> 
LAla.(TNDl(l)*(eCaC0l*CV+i,P*AC)*TN01(?)*(r01-CV*,5*AC 
l)^PR2*PRl42»CJ/<»CeC02+CVel,5*AC) 

C 

C iteration  LOOPoHEaD: 

ICDUNTal 

EPSIalo 

TIGaill^OELTI 

5 CONTINUE 
C 

C temperature  PRCPlLE^nRlED  LAYERI 

R01s,5»(l,eSl)»(S2aSl) 

RD2a,!>»(l,,S2)»(S2-Sl> 

SD2B(l,eS2)*(l.sS2) 

SD«SD1*SD2 
DO  b Ka2,KK 
PKaA4*FL0AT(Kel) 

RD«DC(  K )/A2#D'’TaU/DELZ/DEL? 

DRDa,25*(DC(K*l )-OC(KBl) )/A2*DETAU/DFL2/DEl2 
AL(K.l) •-(2.*RD*S0)/(RnB0RD^PK*RD2j 
8E(K-l)a(R0*OROoPK*RO2)/(RDaDRD*PK#RO2) 

6 GA(K-l)a*(  (RDi05D*PK*R01)*TN01<K-l)M2,*RO-Sn)* 
1TND1(K)*{RD^DRD-PK#R01)»TN01(K*1)*QD1{K)*S01»0D?(K)# 
2S02)/(R0o0RC*PK*RD2) 

CALL  ERMA(kK,OE1,LA1,0, ,TIG,TND2> 

C 

C temperature  PROFiLFeFROZEN  COREI 

RF2a,;>«S?*(S2oSl) 

RFla,5*al*(S2oSl) 

SFc-b2»S2 
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SF»3F1*SF2 
00  9 h«2|MM 

rf«fc(m)*oftau/dely/de:ly 

DRFi,25*(FC(M*l).FC (M«l) )«0ETAU/DELH/0ELH 

AL<M*1)«-(2,*NF*SF)/CRF-0RF‘*FL0AT(M-1)*RF2) 

BE(M-l)»(RF*0RF*FL0AT(M8l)«*r2>/(RF*DRF-FL0AT(M-l)*RF2) 

GA(M-l)«-((RF»ORFa''LOAT(M-l)*RFl)*TNFl<M*l)-vii*RF-SF) 

l#TNFl(M)*(RF*0RF+FLOAT(Mel)*RFl)*TNFl(M^l>*QFl(M)*SFl* 

2UF2(M)*SF?!/(RFsDRF-FL0aT(M-1>*RF2) 

CALL  EmMA(HMis1,  .01  ii"i  »T1C|TNF2) 

THERMAL  EQUILIPRIU^  AT  INTERFACEj 
TE«T0^TIG  (TRsT?) 

CE*FE(TE) 

CNl*(CE»CR)/(CTPaCR) 

approximate  diffusivityj 

DO  90  K*1,K1 
TO«T0*TND2(K)*(THtT0) 

PREspR-PwR.PGC*(C,>CN'2(K)»(CTPeCR)  )«TD 
DI?(K)aO(PRE)/FD0 


concentration  profile 

SO«SO*POR 
00  10  K82.KK 

RClaOIKKI^OETAU/OELZ/OELi 

RC2iOI2(K)*DETAU/DEL2/OELE 

0RCl*,25*(0Ii(K*l)s0Il(Kel) > «DET AU/DEL2/DElH 
DRC2s,25>*(0I2!K4l)-0i2(K-l)>«DETAU/DEL2/0ELZ 
AL(Kel) a-(2,*RC2*Sn)/<RC2-OHC2*FLOAT(K«l>*RD2) 

BE ( K-l ) s (RC2*ORC2«FlOAT (Ksl >*R02>/(RC2-DRC?*FL0AT(K«l)*Rn2> 
GA(K-1)5*{(RC1s0R"1*FL0AT(K-1)*R01>*CN1(K-1>-(2,*RC1-SD>» 
ICNKK  , ♦ RC1*Q’'C1»Fl0AT{K«1)*R01  )*CN1(K  + 1>>/(RC2-0RC2*FL0AT( 
2Kel)*R02) 

call  EMMa{KKi0,.0.i0.iCnIiCN2) 


ENERGY  BALANCE  A THE  InTERFACEI 
GRT02«(TN0?(KK>1):  TND2(KK) )/OELZ 
GHTF2»(TNF2(:-lM*l)-tTNF2{MM))/DELY 

GRC2sICN2(KK*1)-CN2(kK) )/0EL2 

C0EF2«-'A5*GRT02*A6-.25*a7«(A2«(2,-DELZ>*DELE*GRTO2»(2,- 

10ELY)*DF.LY*GRTF2) 

CONOl^o.StllFCIMM  + D+  FCIMMjj/SlaGRTFl+IOClKK^D+OCIKK))/ 

l'a,-Sl)*GRTDl)  ^ 

C0NO2a-,5*( (FC(MM*1)*FC(MM))/S2*^GRTF2*(0C(KK*1)+DC(KK) >/ 
1(1,-S2)*GRTD2) 

PH0018,5*A3*( { .75*0P1<KK^1)*,25*DP1(KK>  >*(l.-Sl)*DEL2*< , 
175*FP1(MM  + 1)*.25*>FP1{MM)  )«S1*0ElY> 

PR0D2».5*A3*{(,75*0P2(KK»-1)^.25*DP2(KK>>*(1.-S2)*DEl2*( 
1,7&*FP2(MM*1)* ,25*FP2(mM) )*S2*DELY> 

ANEARb.0625*( a2*(2.»SI=S2)«CELZ*(TND2(KK>-TND1(KK) > v(S2* 
ISi; •DELY*(TNF2(MM) "TNFI(MM) > )/oetau 

AlNTER* .1875»( A2*(2.-S1*S2>*DELZ*(S2+S1>*DFLY'/DFTAU 

OH1*.I>»A0*(DI1(KK*1)*DI1{KK))/(1,"S1>*GRC1 
0R2«.!>*Ae*(0I2(KK*l)*0l2{KK)  >/(  1,-52)  *GRC2 

OELT  'Cp .5* (C0N01  + C0NG2eC0EFl*DRl-C0EF2*0R2*PR001*PR002o 

12,*ANEAF )/AlNTER 
TI»TI1^DELTIC 


ITERATION  LOOPoTAILI 
IF(  ICOUNT.EQ.DGO  TO  12 


1F( IC0UNT,lT.3>G0  to  ll 
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t. 


t 


11 


12 


C 


C 

14 

15 

16 

C 

r* 

V 


150 

C 

C 


160 


C 

C 

161 


TC1»TC2 

TC1»TC2 

TC2»TI 

TC2*TIG 

IF(  IC0UNT.lt, 3)G0  TO  12 
NY1»1,E7*aL  o( TG2-TG1) 

IF^NYI.lt ,15)G0  to  12 

SL0»(TC2-TC1)/(TG2^TG1» 

NR«l.E7*A8'!(bL0-l) 

IF(NR.LT,1?)G0  to  14 
X»1,/(1.-Sl0) 

CONTINUE 

NY«l.E7«ABS(TIGaTI  ) 

IF(NY.LE.EPSI)G0  to  i5 

TIG»(1.-X)*TIi,*X*TI 

IC0UNT=IC0UNT*1 

control  iteration  convergence: 

IF( lCOUNT.FO.10)EPSls50 

IF(lCOUNT.EO.30)EPSl=la0 

ir(lCOUNT.EQ,;j?)EPSI*500 

irdCOUNT.EO.  100)  EPS  1*1000 

IF( ICOUNT,NE.1000)GO  TO  5 

I«RITE(5»100)  n 

RETUPN 

ITERATION  loop-out: 

TI2i(TIG*TI  )/2, 

GO  To  16 

TI2i(1,-X)*TIG*X*TI 

C/ELTIaTI2-TIl 

AVER2*AVER2*Fl0AT( I COUNT ) /FLOAT ( N ) /2 . 

HEAT  CONSUMPTION  AT  THE  INTERFACE** 
0Sl«C0NDl-(  C0EF1-A6)  #0Rl  i-PROUl 

ClS?sc0N02-(C0Er2“A6)*0R2*PR0D2 

ACCUM«aINTER*DElTIC*ANEaR 

OI2«OS2>OS1-2.*ACCUM-OI3 

IF ( ITEH.GE.2)G0  TO  150 

0lGaQI2 

ITERalTER^i 

GO  TO  70 

CONTINUE 

0ELQIaQl2-aii 


N1»(N^509)/§S^'’  ‘^^'^^f’^^ION/SUBLlMATION  POWER  CONSUMPTION 

FMWPa(FP2(i)«EP2(Ml) )/2* 

DO  160  Ka2,KK 
0MWPe0MWP*0P2  < K ) 

FMWPaFMWP*FP2(K) 

DMWPeDMWP#FP0#THlCK*(l,-S2)*DEL2*4.1«4 

FMWPaFMWP.FP0#THlCK«s2*DELY*4,184 

TMWPaOMWP^FMWP 

TMWEaTMWt*TMWP*OELT*50, 

C0NSaHS*DR2*4, 184 


[Ir  characteristic  VARIABLES  HISTORY) 

IF<N.NE.1)G0  to  161 
WRITE(20,201) 

Kal 
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f 


162  K«K>1 

I IF(K.CT.K1)G0  to  163 

i IF(TnD2(K) ,GT ,TN02(K-1) )GO  TO  162 

163  TOMiTNO2(K-l'*(TR-T0)*T0-273,l5 
EM«(K-2)*DE:l2 

TC«TnF2(1)*(TR-T2)+T0-273.15 

TT»TI2*(TR-T0)^T0-273.15 

i WRITE(2kJ,202)  N , T TME  , A VER2 , 2M , S2 , Tc , TT  , TDM  , FMWP , JMWP , CONS 

i RRITE(21,2P3)  time, 2.M.S2.TC.TT, tom, FMWP,TMWP, cons 

I 164  CONTINUE 

I C 

C FILE  ON  OISK  TEMPERATURE  AND  VAPOR  PRESSURE  PROFILES! 

N2« ( N-959 ) /500 

lF(N.NE.l,AND,N.NE.5iiJ9.ANO.N.NE,  709.AND,N,NE.7D9, 
lAND. (N-959)  .NE .N2*500)GO  TO  166 
1F(N.EQ.459)G0  to  166 

; WR1TE(23,2P4)TImE 

DO  165  Lsl,Kl 
E»OElE*(ls1) 

TF«TnF2(L)*(TR-T0)*T0=273.15 
T0«TnU2(l  )*(  TR<-T0)^T0-273.15 
PWaPGC*(CN2(L  )*(CTP-CR)-»'CR)»(T0^273,15) 
C0NCsCn2(L)*(CTp-CR)*CR 
U'R1TE(23,206)  E,TE,TD»Pw.CONC 
I 165  WRITE(22,205)  2,rF.TO»CONC 

f 166  -ONTINUE 

' C 

c check: 

HRITE(19,APi0)AVER2iTNF2(l),TNF2(10),TNE2(2l),TNO2(lb) 

; 1.TNO2(1),CN2(10) , Cn2 ( 2l ) , GRC2 i GR TD2 , CRTF2 , X , DELS 

400  FORMAT(^6. j,11F10.5,e14.5) 

C 

C UPOATE  ALL  PARAMETERS? 

. 00  17  Ksi.Kl 

TND1(K)«TND2(K) 

17  CNl(K)*CN2(K) 

DO  la 

18  TNFl(M)aTNE2(M) 

GRTDltGRTD2 

GRTFle&RTF2 

GRC1sGRC2 

COEFiaCOEF2 

SD1SSD2 

SF1#SF2 

DELS*S2-S1 

S1«S2 

TIHTI2 

011*012 

FIEL01=FIEl02 

AVER2*AVER2#FL0AT(N)/Fl0AT(N*1) 

C 

c simulation  loop-tail: 

IF(N.GT.10000)RETuRN 
IF(S1.GT.0,001)GO  to  19 
TENeOENS*PDR*TH I CK*HS 
EUF«TEN/TMWE/4 . Ib4 
HRITE(20,207)TIME,EUF 
RETURN 

19  CONTINUE 
IF(N.EO.100)GO  TO  ?1 
IF(N,EO.500)GO  to  2l 


GO  TO  4 
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21  DETAUPDETAUPie. 
rELTpDELTtia, 

X*.2 
GO  TO  4 
C 

c 

100  formatc  iteration  stuck  at  time  Increment  n»  m5,'j  co 

IMMANO  IS  RETURNED  TO  THE  MAIN  PROGRaMM') 

101  FOPMaTC  melting  occurs  at  Si  A9  CENTER  a T 1 1- £ » ' , F6 , 3 i ' 
IHRS,  i-  ',17,/'  Cum.-iAnd  is  RETUtiNED  TO  flAI.^i  PROGRAM') 

200  FORMATC  experimental  CONDITIONSt '/3X'FIELDs  ',EH.4/3X' 
IPR*  ' . - ■*.4/3x'PwR=  ',E11.4/3X'  r<«  • ,E11,4/3X'T0*  ' 

2iEll. a/3X'TH1C<NESSs  ', Ell . 4/3X ■ WATER  CONTENTs  ' 

3.E11. 4/3X'FAT  contents  ',E11,4/> 

201  FORMAT!//'  time  H 1 STOR Y : ' /2X • N * i 4X ' T 1 ME ' , 1 v ' A VER  . ' » 2X ' ?M 
lAX ' , 6X'S2' ,8X' TFMAX* .9X'T1 ' ,8X' TDmAX ' I0X 'FMWP ' ,PX' TKWP' 
2i9X'QI ' ) 

202  FORMAT! 15,3F6.2,7(1XE11.4)) 

203  F0RMaT(2F6.2i7E11.4) 

204  FORMAT!//'  TIME*  '.F6.2,'  HRS ' /IX ' ? , Y ' , 7X ' TF ' 1 10X ' TO ' 1 10X  ' 
IPW  ,11X'C  ) 

205  F0RMaT!F5.2,3!1XE11,4)  ) 

206  F0RM/.T(F5,?,4!1XE11.4)) 

207  FORMAT!//'  TOTal  DRYING  TIME:  *iF5,2,'  HRS'/ 

1'  MICROWAVE  ENERGY  USAGE  FACTOR:  ',F5,3) 

end 
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- NOMENCLATURE 


Letters : 

2 

A sample  cross  section  area,  cm 

3 

C concentration  oi  water  vapor  In  the  pores,  g/cm 

specific  heat  of  the  dried  material,  cal/g/*K 
Cp  specific  heat  of  the  frozen  material,  cal/g/*K 

C specific  heat  of  the  water  vapor,  cal/g/’K 

^ 2 
D effective  dlffuslvlty,  cm  /sec 

D°  effective  dlffuslvlty  at  1 atm.,  cm^/sec 

2 

effective  Knudsen  dlffuslvlty,  cm  /sec 
E electric  field  peak  strength  In  the  vacuum,  V/cm 

uiwctric  field  peak  strength  at  breakdown,  V/cm 
E^  ..'lectric  field  peak  strength  In  the  dielectric,  V/cm 
f microwave  frequency,  Hz 

fg  Ice-vapor  equilibrium  relationship.  Relates  the  concen- 

tration to  the  temperature 

f^  Ice-vapor  equilibrium  relationship.  Relates  the  normalized 

er'-.icentration  to  the  normalized  temperature 

2 

h heat  transfer  coefficient,  cal/sec/cm  /*K 

AHs  enthalpy  of  sublimation,  cal/g 

k effective  thermal  conductivity,  cal/sec/cm/*K 

3 

K dissipation  coefficient  defined  by  Eq.  (3.6),  cal/sec/cm  / 

(V/cm)^ 

m weight  of  the  sample,  g 

A overall  drying  rate,  g/sec 

P total  pressure,  mmHg 

partial  pressure  of  water  vapor,  mmHg 

2 

Q heat  flux,  cal/sec/cm 

S normalized  Interface  position 

t process  time,  sec 

T temperature,  ”K 

V volume,  cm^ 
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2 

W •ffactlv*  mass  flux  of  watar  vapor,  g/sac/cm 

X apaca-coordinata,  cm 

X Intarfaca  position,  cm 

x^  initial  moisture  content  of  beef  sample  (wet  basis) 

x^  initial  fat  content  of  beef  sample  (wet  basis) 

y normalised  sfxxce -coordinate  in  the  frozen  region 

s normali;:ed  space -coordinate  in  the  dried  layer 

Greek  letters: 

2 

a effective  thermal  diffusivity,  cm  /sec 

p Biot  Number 

g thickness  of  the  dried  layer,  cm 

€o  dielectric  constant  in  vacuum  (CGS) 

€'  relative  dielectric  constant 

€"  relative  loss  factor 

0 normalized  temperature 

^ fraction  of  the  initial  mois;ure  content 

corrected  fraction  of  the  initial  moisture  content  defined 
by  Eq.  (9.3) 

^ acceleration  factor  in  the  iteration 

fr  number  Pi 

3 

p density  of  ice,  g/cm 

3 

P-  density  of  the  dried  material,  g/cm 

^ 3 

Pp  density  of  the  frozen  material,  g/cm 

a porosity  (void  fraction) 

Cli  electric  conductivity,  mho/cm 

*e;  normalised  time 

^ total  microwave  power  absorbed,  cal/sec 

3 

<ji  density  of  microwave  power  absorbed,  cal/sec/cm 

2 

microwave  power  absorbed  per  unit  area,  cal/sec/cm 

Indices  (superscripts  or  subscripts): 

0 initial  value;  before  the  microwave  power  is  turned  on 

1 starting  value;  immediately  before  the  microwave  power 
is  turned  on 
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2 final  valut;  as  tha  frotan  cora  dlsappaars 

i final  valua;  at  tha  matarial  is  fully  driad 

in  tha  driad  layar 
F in  tha  froaan  ragion 

i at  tha  intarfaca 

R in  tha  vacuum  (room) 

S at  tha  outar  surfaca 

max.  maximum  valua 

* normaliaad  valua 
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